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ABSTRACT , r . 

This is the secc^'vi of a two-volume minicourse xn 
astrophysics.' It con-t;ains chapteiis on the following topics: stellar 
nuclear energy sources and nucleosynthesis; stellar evolution; 
stellar structure and its determination; 'and pulsars. Each chapter 
gives much technical discussion, mathematical treatment," diagrams, 
and examples. References are -included with each chapter, , (BB) 
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S!^LLAR IfUCLEAR EHEROY SOURCES • 
lUCLEOSYHTHSaia 

It Introduction . ^ . . 

' — ^ /. ■ ■ 

. ■ / . • . • . 

In this odnltext the subject of nuclear reactlone In stars ylll be 

surveyed vlth three alms In mind: to shoy what nuclear processes play 

Inportant rolm In generating the energy radiated by stars during the 

principal stages of their evolution, to see vhat physical factors govern 

the rates of energy production by these processes, and finally to demon-r 

stx^ate the roles these processes play in building up, heavy elements froa 

the main primordial material of the Uhi verse-, hydrogen and helium. 

Two quite dlff extent sources of energy supply the energy radiated hy 

' • (>■ . 

a star: gravity and nuclear fusion. As a rule these sources alternate 
in Importance In successive steps in the evolution of a star from birth 
to Its final state. Initially gravitational contraction condenses and 
hearts interstellar gas and dust to form a star, and when this jpracess has 
raised interior temperatures into the range of «10^^K, fusion of hydrogen 
into helium by various processes takes over the Job of providing the star' 
energy. At this point gravltationq|Pc6ntraction ceases and gravity no 
longer plays ah active role in energy generation. However, when the hy- 
drogen in the central region of the star has been converted into helium, 
the fuel ifor the initial forms of nuclear energy generation is exhausted.' 
These fusion processes then cease and gravity takes over, "catising the star 
to contract and heat further until tenqperatures of the order of 10*^K are 
reached. At these temperatures helium nuclei can fuse into carbon nuclei 



with peleaae of energy, and a nev. source of nuclear energy e«n^i:;wi« i;he pic- 
ture. ' 

Eventually the fuel for this process is exhausted. Qft^tevV i. > t v^ i in 
once more end heats the star until anothex" exothermic nuclear rG>i?o!vl^: r 
hecoihes possible, and the process repeats in t>vls fashion until it in^ 
terrupted by one of ^several possible terminal steps which vill be ^iUscusse 
later* 

The key factor in production of energy by nuclear fusion is tVie fact 
that heavier nuclei are bound more tightly than lifter nuclei, and when 
they form from fusion of lighter nuclei, tiit^ difference in binding energie 
of the initial and final nvi^ltti is given off in. the form of energetic 
electrons, positrons* gansa rays, and neutrinos « The connon mahifestation 
af this fisustor can be found in the familiar chemical periodic table of the 
elements where the mus of a helium atom is less than the total mass of 
fo\ir hydrogen atoms that can fuse to form one helium. This mass defect, s 
it Is called, is evidence of the binding energy increase, and the energy 
glvea off in the fusion process is E « Anc^ by the Einstein relation. 

To give an explicit example, consider the fusion reaction . 

hn^ He* ~" ~ ^ ; — 

1* x 1.008 U.003 + .0029 gm 

mass of k mflUBS of 1 mass equivalent 
^ moles^ of mole of of the excess 
hydrogen helium energy 

The excess energy equivalent .0029 gm/wole of He formed is equal by the 
, Elnstedn relation to 2.6l x 10^' erg8»or 6.3 x 10^^ ergs/gm of He formed. 
This is quite a lot of energy; in fact it is approximately 10^ times the 
energy released in a typical chemical reac:t ion, even so energetic a reac- 
tion as the explosive burning of hydrogen in air. Such an energetic proc« 



iB necessary, however, to supply the obsenred energy production rfttaa of 
stars. 

The Sun. for example, radiates k x lo" ergs/Bee from a isBS of 
2 X 10" gms.. or about S^ergs/sec per graa. It can continue at this 
rate for some 10 billion years or «3 x 10»' see. yielding a total of 
6 X 10»' ergs/gm. in good agreeaent with the value calculated above for 
energy production from hydrogen fusion. 

Gravitational contraction is an efficient energy source, but It oper- 
Hites on a much shorter time scale than'^nuclear energy production. For 
example, gravitational energy released by contraction of the Sun from an 
initial gas cloud to Its present siie would fuel the Sun at its present rate 
of radiation for only 20 million years. But geological evidence on Earth 
ihows that the Sun has been radiating at its preaent rate for at leMt 
3.5 billion years, and present theory of stellar evolution predicts that 
a star like the Sun will radiate at nearly its present rate for approxi-. 
mately 10 billion years before undergoing any substantial changes. There- ^ 
fore, nuclear burning stages in the life of a star tend to be relatively 

p. 

long and stable while gravitational energy production stages are relatively 
Shorty aad often dynamically unstable. 

The general processes -of-fusionaeadjto prod^^^^^ heavier elements 

from lighter ones. The building up of heavier elements proceedi^wi^^ 
in the sequence H He c N > 0 y Me Mg Si ..... Pe. Many 
steps and details have been omitted in this outline, and some of these will 
be filled in later. There are several features of this bare outline of 
nucleosynthesis thaWeserve comment at this point. First, the sequence 



St ' ' ii ■ u . 

stops at Fb . Thl8 nuelous is the most stabije one in nature, and the 
conBtxnietlon of any heavier nucleua than Pe"Vequlre» a net input of 
energy rather than giving off energy, ao energy production hy fuaion stops 
at Pe". It' is possible to Bake heavier eleaontu than iron in stars but 
the process becomes Bore difficult and requires rather special conditions, 
as ve shall see . 

The second point to notice i.n the building up scheme above is that 
there are lots of gcqae in the periodic table of elements made by simple 
fusion. For Instance , "there is no Na«^ or Bo"'or Li', etc. , in the sequence 
above. These ' eleioents are actually produced in small amounts in some of 
the nuclear raaetlons that fil^ihe main products listed above, and ve 
shall show hov these minor constituents are created. Toward the heavier 
end of the- sequence there are a number of nuclei that are somewhat neutron- 
rich in that they tiave more neutrons then protons. These elements cannot 
be create^ without the presence of free neutrons , and some of tl^ steps in 
production of intermediate elements Just beyond atomic weight «20 do occa- 
sionally yield neutrons which can be captured, to give the isotopes missing 
from the preceding list. 

II. Nuclear Reaction Rates 

It is clear that while stars last a long time by our reckoning, they : 
ar^. nevertheless dynamic, evolving objects. Consequently ^ the rates at 
which various processes that go on in a^rt in 
understanding its. behavior. In particular, we will want to know the rates : 
of poasible nuclear reactions both from a standpoint of the. stars' energy 
production and in order to understand how elements are synthesized by nuclear 
reactions in stars. ■ 



ThB rate of a partlciilar nuclear reaction la the product of several 
tors ^hat are not difficult to understand In broad tennB, although bom 
of these factors may he quite eoniplex In detail. In Its slioplest form* a 
nuclear reaetlcm aay he considered as the possible result of a dolllslon 
betveen tvd particles » These particles nay be of various types Including 
nueleons (i.e. » protons » helium nuclei (alpha particles) and other nuclei 
such as deuterium^ He\ C^^» etc.), and lighter particles such as electrons » 
positrons, neutrinos 9 and ganma ray photons. 

It la. easy to see that the reaction rate (nuaber of nuclear reactions 
per second per unit vo].uae) nnist be proportional to the collision rate, 
the number of collisions per second per oi'. This nay be calculated in a 
straightforward way. Let the effective "siae" of a particle of one kind, 
say Type 1, be described by radius r^ so it may be considered to collide ^ 
with any other particle it approaches closer than r.. If this particle " 
moves at velbclty Vj^, in one second it sweeps out a cylinder in space of 
height and balse area wr^ which we designate a^. This "collision 
cylinder" of area o, , length v, , is* such that the oype 1 particle' will 
Collide with anything within it. The area la called the collision crosa 
section of a Type 1 particle. The cylinder swept out in one second has a 
volume ^2y^» (See Plguare 1) 

■» 

Pig% !• Colliaion cylindtir sw«pt out in on« second by a Typ# 1 particle. 



If therm are Ng partlelea of Type 2 per cm' in the gas, a Type 1 
particle vill coUide with aiVj^Hg 'PyP* 2 particles each second, and if 
there are particles of O^pe 1 per cm* the total collision rate will 
be OiVj^Nj^Ng. This axaaiple was phased on the aosuaiption that every Type 1 
particle moved at velocity v^ and that all Type 2 particles are at rest. 
Inreality^both types of particles have a wide range of velocities, and 
a suitable average of the product ((Jv)niuBt be made over the distribution 
of particle velocities characterised by the gas temperature T. This vel- 
ocity distribution is the MaxwellrBoltsman distribution lAich gives t^he 
fraction of particles, n(v), having speeds between v and v + dv 

n(v) dv - Uir /JL_\'^' ^V^'/2kt 4^ (1) 
\2iTlrty 

sfollisions between nuclear particles take place on so small a sise 
scale («10-*' cm) that quantum effects are very important, and the nimpln 
classical picture above must be modified. The principal change is that 
nuclear dimensions such as must be expressed in terms of the correspond- 
ing^ deBroglie wavelength, h/mv. This means that the collision cross section 
o will be proportional to the square of the deBroglie wavelength, i.e., to 
v"', and the product ov in the collision rate must be replaced by v~*. 
This makes the collision rate Rj^ proportional to v~* , or 

Rl2 ■ ANj^ll2 v~* collisions per cm' per second (2) 
vheire A is a constant that includes Planck'.s constant and the. masses of the 
colliding particles. It would be premature to iaiverage (2) over the velocity 
distribution (l) at this point because there are. more velocity dependent 
factors yet to be included in- the overall reaction rate. 



We must txmi'noxt to sone baalc facta of nuclear phyBica to aaa hov 
cloae nuclei nuat approach each other for a nuclear reaction' to take place. 
The central feature of nuclear forcea, for puriMiBea of this dlacuaaion, 
ia that they are Tary short range forcea and quite strong. He are accus- 
toaed to forces like gravity and electrostatic fbrces vhose potentials 
faU off slowly with distance in r-» fashion and to dipole forcea that 
fall off as or eren vender Haals forces between noleculea that gp as 



The forces between nuclei fall bff much more rapidly than any of 



these, and it is a fairly good approximation to think of the nuclear force 
as being derived from a "square-well" potential like that shown in Wgure at 



r. 



ui o 
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Pig* 2e Nuclt4r "aquart-wtll" potential 
The range, r^, ;of this nudlear force is of the order of lO"*' cm, the sise 
of a nucleon. and the force is attractive, meaning that once two nuclei 
touch they coalesce much as do two small droplets "of liquid that ccblnel 
under the influence of surface tension forces once their surfaces come ' 
into contact. This^alogy with a liquid is widely used in nuclear ihy.io.,. 
and it is a useful way to. think about collisions leading to nuclear reactions. 

-Most of the nuclear processes- of interest in stellar Energy production 
involve, mtei-act ions between charged nuclei such as protons . alpha particlw 



(doubly charged H<» nuclei ) , and so forth. These particles , beln^ posi- ' 
tively chared, repel each other strongly at close range, and since they 
must approach to distances - 10-»» cm to react, they must ha^ high ini- 
tial velocities (i.e.. high temperatures)' to overcome the strong electro- > 
static Coulomb repulsion. The Coulomb repulsion ot iyo protons separated, 
by 2 X 10-»» cm is « .6 Mev (l Mev - 1.6 x 10-» ergs).' This correspond, 
to two protons colliding head-on each with initial velocity « 1.2 x 10» 
cm/sec which is .Ol» times the speed of light. In terms of the average 
energy. | M. of a particle in a gas at temperature T (k is the Bol^izmann 
constant I.38 x 10-»« ergs/degree), the Coulomb repulsion correspond^ with 
the average energy of hydrogen gas at 5.5 x 10 •'degrees K. Such tempera- 
tures are far above those foumi in the interiors of normal stai-s whose 
central temperatures range mostly around 10-30 million °K. It is true 
that a few veiy high speed particles are present in a gas at a given temper- 
ature. even pwiiiclea whose kinetic energy are as much as UOO times the 

average. particie energy as required for the proton-proton collision in 
the example. But the fraction pf such energetic particles is vexV small ^ ' 
indeed (« e"-" « lO"*"). much to8 small to allow nuclear reactions to 
proceed at a rate sufficient to account for the energy output of a typioai " 
star.. There must be some way by which the Coulomb barrier -a be foiled, 
and the method was first discovered by Oamow in connection with a related 
problem", natural radioactive decay of heavy elements W alpha parti«a,e " 
emission. G«««>ow showed ^that^ charged particles which do not have enough^ ^^^^ 
energy to go over the barrier maj^. instead' go through it by the process of 
quantum mechanical tunnelling as shown schematically >q Figure 3 (next page). 



Fig •3* Schema tc4.c^ r«pr«aontatlon of tunnallng of a chargad 
partlcla through a CoaloMb barrlar* 

• ", • <r V • ' . ." 

Quantum nechanlca deacrlbes the poaltlon of the incident particle 
in the fonJ^pf deBroglie vavee of vavelength «« (e - O^r) )^ vhich becoMa 
Imaginary within the barrier whe^e E - Uts-) is negative. Elementary theory 
of wave propagation shows that such a wave will propagate with eicponen- 
tially decreaaing amplitude Into a region where its wavelength is imaginary. 
This is true of electromagnetic waves (e.g., penetration into metals to the 
"skin depth") and acoustical waves (penetration into dissipative Mdium), 
to cite two familiar, examples, and it is also true of matter waves (deBroglie 
waves) . Therefore » the matter waves Associated with an incident charged \ 
particle penetrate into the barrier between r^ and r^ irtiere the wavelength 
X, is imaginary, but the sinplitude decrea^ exponentially with distance into 
the barrier rather than going to zero abruptly at E » U(^) irtiere claasl- 
cally the incident particle would stop. This means that there is a small 
but finite probability for the incident partite to tunnel through the 
barrier and bo f6und inside r^. This probability for tunnelling ciin be 



calculated quaxitum laechanlcally and tha result it that the barrier trana- 

miaalon (fraction of incident particles that get through) is 

T « exp (-Uir'z^ZgeVhv) * (3) 

0 ■ ■ , ■ . ■ • • 

Once the incident particle has penetrated the barrier, there is no 

certainty that a nuclear reaction vill take places the particles might 
separate again almost at once , or some other nuclear reaction than the one 
of interest might take place. A factor, F, is therefore included in the 
reaction rate expreaaion. It is called the "nuclear factor" and ia the 
probability that the desired reaction vill occur once the Coulomb barrier 
has been penetrated. This factor is very difficult to calculate vith any 
accuracy at the present state of knowledge about lev energy. nuclear reac- 
tions > and it is usually determined experimentally « 

Putting all the factors together, the overall reaction rate can be 
expressed in the form 

R . J (c V3 V-') . (v»T-/» . (.- ::^^) r dv! U) 

° \ 1 V \ . 

Collision Maxvell-Boltzmann Gamov factor Nuclear 
factor velocity dlstri- (barrier factor 

bution ^ transmission) 

It is evident that integration of the above expression gives a some- 
vhat complicated result. For particle energies irtiere T does not vary rapidly 
with energy (i.o., where there are no "resonances"'' or energy regions of 
enhanced nuciear reaction probability) the rate of energy production, e, 
which is proportional to (U) takes the form 



vfaere B and b are constants, p Is the gas density, X^^ and are the frats- 
tlons Xsy weight of the Interacting nuclear species (lypes 1 and 2 In the 
earlier discussion), and T Is the absolute tengperature . Even this expres- 
sion Is Inconyenlent for aany uses and can be. expressed In still simpler 
ofom over a limited temperature range, ^ ^ . 



e - ApXj^XgT^ 



(6) 



vhere n Is a positive number that varies soMvfaat vith energy but vhich 
can be taken as constant over a reasonable teBqE>erature range. The exponent 
n Is large ranging from about k for the pp cycle described below at testper- 
atures around 2 x lO^^K to about 20 for the CNQi cycle at temperatures arouod 

1.5 X Ip^OK. ^5 

' .. ' ■ - ' " - • 

HI . Nuclear Energy Production In Stars s 

The simplest and most efficient energy producing nuclear reactions In- 
volve the building up of He** nuclei from HV nuclei (protons) by two princi- 
ple ittodes. One of these, the p-p chain, involves only light nuclei up to 
.He** in its principal form while the other involves carbon,- nitrogen, and 
oxygen nuclei and is comc^nly referred to as the CN chain, or CNO bi-cycle 
since it consists of two Interrelated reaction cycles. 

The p-p chain in its principal form (p-p I) consists of three steps 
that take place in series as shown below; 

(1) + + e^ + V 

(2) . + He^V+ Y ) P-P I chain ^ (7) 

(3) He^ + He^ He** + 2H^ 

In these reactions is a deuterium nucleus, e* iti /a positron, y is a g 
rsy, and v is a neutrino. Other possible reactions such as He', 
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.vid«.t th.t th. pv.r.u *ff.ct Of tM. ««tlo. cum 1. to f... four 
protons into an alpha particle* . ^ - 

_ - loj. 
I>» 2 p..Uronr.ooo find .Lctron.. aanlhllat.. „4 n,I4 a y.. 
ror . total 6 ^ pi„. 2 

(He*) foraed. ®' 

Th.« .r, t» addltlo^l 8„«p. Of r,^,,„ ti»t t^«rti, 

fo» nucl.1 Of .o^ut M^o, ^. to,, B... „i a«) W fouov- 
ing chains: o 



He'c + He* ■*■ Be' + y 
Be' + e- +^ 
Li' + H* Be* + Y 
Be" 2 He* 



P-p II chain (^j 



-n un.t..l. .1th a half-Uf. of only 2.6 V„-.„o. Io.t..d .of. 

captju-ing a frise electron, the Be' in th« * ^ , \ - 

, in the first^step of (9)'can coliidp^vlttf 

a proton to form boron:' / " 



' Be' + H* -»• B* + Y 
B* -^'Be».+, 
Be' -»• 2 He* 



B-VBe-:.e**v , , ^^^^ 



The sloveit step in the p-p chain Ji- the "initial one. 2 H» D» .* + v. 
and this rate limiting step controls the whole chain. 

It is possilile to set up a group of linear differential equation, for 
thV time-variation of the ahundances of each species in the p-p chain, aho^. 



These equations resemble those of radioactive decay equations » and their 
solutions are In the/fcrm of aseries of eaqponential buildup and decay 
terns for the individual species. These solutions of the abundance equa-. 
tions are useful in dealing vlth rapid phases of stellar evolution that 
take place on tine scaleis too short to allov equilibrium to be established 
among nuclear speciies. Or the equations can be used to predict equilibrium 
abundances among the various nucleax: species during sloy stages of evolu*-^.;, 
tlon. These equilibrium abtmdances are inversely proportional to reaction 
rates involving the species, ijfuclei that have slow-reaction rates will be 
abundant ybile those that react rapidly vlll have cozl'espondlngly low.coa> 
centratlons, and .these relative abundances are loGPortant in establishing 
the buildup of heavier elements by the various nodes of nuclecfsynthesls . 

A star containing only H and He (ie.g. , extreme Population II stars)-, 
must operate on the p-p chains, but there is another option ayailable, for ^ 
jtiurs^wlth an admixture of heavier elements, specifically C and H. In 
tils case C*'' and^H^*~serve_as_^atd3^^ a double cycle of reactions 1 
by which k ♦He'* as befoue. The basic CaKr~cycle-t8^ follows: 

V.^C*^ + H* -*.'h*' + y . ] 
N*' -»■ + e+ + V • I 



r 

• I 

Y ♦ 

I' 



•i ('♦•N^" + ^H* ol' * 

I 0*' N*" +'6+ + V \ ' 
' N»» + H» - He" + C»2— ' / 

• or.P*- '+ Y . (11) 

I 0*' + H* P*' + Y • . . 



• I Ol' + e+ + V ? 

•1..' ■ - ,■ r 

. V ' I + H» He" t N^" • 

' - . '» V or I 

O .. ■ .. - \ . . 

ERJC . F*'+ Y ° 13 



The rate limiting, react ion is N** (p, y) 0*' so N*** eventually becones ' 
thie most abundant species in the bi-cyole. The initial C** and C*' nuclei 
become converted in time mostly into N***. • 

Be8^4ea.^^h^^]fimary fusion reactions contained In the p-p chain and " 
CHO bi-cycle there^r^other reactions that fuse other nuclei than to 
produce significant contributions to stellar energy generation under cer- 
tain conditions. The first of these to be considered is the fusion of 
three alpha particles (He** nuclei) to form a C*^ nucleus. The.stral^t- 
forward reaction 

3He** C^^ . \ (12) 

is a very •;3Elikely reaction bee aiiso it would involve a a-body collision » 
moreover one in irtiich two Coulomb barrier tuxmeillngs must occur simul- 
taneously. The probability for such an event even under conditions at 
the center of a star is so small as to rule out this reaction as a viable 

energy source. I 

■ ■ ■ . ■'■ ' . . . / 

There is^ however, a way for the triple alpha (3ot) reaction to go |n 
two steps. The first m^es a Be* nucleus via the inverse of the decay 
reaction for this nucleus ^ namely 

. He" + He" t Be' ('13) 

As mentioned earlier, the Be' lasts only a very short time since it a. decay 
lifetlipe is 2.6 X 10"*' sec. In a, hot, dense He* gaa deep in a stellar 
interior this short lifetime is sufficient to give a reasonable, chance ' 
this,t the Be' will be struck by another He* before it decays, yielding 

C** by the reaction . ^ i 

I ■ ' ' ')■■•■ - " ■ : "■ . ; 1 . • ■; 

Be' + He" $ C** +: 2y i (ll») ' 



lU 



Even this gtep would not produce appreciable carbon were it not for the 
fact that there is an excited state of the whose energy level is near 
the combined nMs-energy of Be* + He**. Even eo, the ^ excited C*' state 
dec^B back into Be' + He* most of the tirie and Into a stable C** state 
bnijr occasionally. The large Coulomb barrier for both (a,« ) and (Be',«C ) 
collisions leads to a very strong temperature dependence of the overall 
reaction rate for the triple alpha process. The energy production rate 
for this process can be represented approximately by a TP^ dependence where 
n varies with temperature but ranges between 30 and 1*0 for temperatures 
between 10* and lp*OK where this process is an important energy generator. 

Prom the standpoint of nucleosynthesis the 3oi reaction is very impor- 
tant for it i« the bridge from very light element synthesis to the pro- 
duction of heavier elements. It is a hard bridge to cross, but one that ' 
V can generate substantial amounts of C** in iatef^evolutionar;f stages of 
passive stars. C As the fifth most abundant nuclcsus in the solar system 
and presumeahly also in the Universe, testimony to the operation of the 3a 
process despite the ^rat her unlikely prbspects for its success. 

The third most abundant element is 0*', formed in yart by the CWO 
cycle and in part by the direct capture of an alpha particle by C** via 
the reaction 

C'* + He" 0*' + Y (15) 
In a similar way we can account for Ne** synthesis, 

0*' + He" t Ne^" + y ' (l6) 

Under conditions favorable to this reaction however , the Ne^° rapidly picks 
up an alpha particle to form Mg^" leaving a low abundance of Ne^". This low 
abundance of Ne^" predicted on the basis of alpha capture does not agree ':_ 



vith the observed natural abundance of Ne"; which puts it foxirth in order 
behind 0 and ahead of C, In Population I stara Mg'" la eighth in order of 
abundance, well belov Ne" in contrast with predictions from alpha pro- 
cesses. The order of observed abundances in Pop. I stars is 

• H, He, 0, Ne, C. N, Si, Mg, S, Ar, Pe, Na, CI, Al, Ca, F, Ni (if) 

Clearly the observed Ne abundance must be accounted for by processes ' 
other than (16) above. 

The He burning does not Involve temperatures high enough to build ele- 

nenta much teyond Mg. The Coiaomb barrier for alpha capture becomes too 

great and the alpha capture rate decreases rapidly so relatively, little 

synthesis of lieavler elements can take place before all the He is exhausted 

in the core of a star. The star's central temperature will continue to rise 

as a result of gravitational contraction of its inert c'* core formed as 

He > C, and if the star is massive enough, the core temperature may reach 

» 6 X lo'OK Whereupon C'^ begins to combine with itself in a carbon-burning 
'phase' '/ . . •■ 

C** + C" -»• Mg*'' + Y ^ " ' 

■■^ ■-»• Na*' + H^: -■ ,v 

Ne*"' + He"- 

Mg** + n 



.23 . _ (18) 



-►'vO** + 2 He** 



One of the reactions above yleida free neutrons , and there is another 
Important source of these particles through the reaction chain 
^ C* + N»? ^ C'' + e+ + V 

C' V + He" - 0'«' + n * ' ^^^^ 



T^e proton that is captured by C*« to form H*' in the first step coaea 
from the second line of equation (l8)^ 

The free neutrons liberated in these reactions are available for cap- 
ture by virtually everything present in the gas, and this capture is rela- 
tively easy because there is no Coulomb barrier between a neutron' which !• 
a neutral particle and an ion. The neutrons are therefore available to 
synthesize a variety of heavier elements. . \ 

The alpha particles liberated in two of the C-buming steps qan^ 
captured by C>% 0>«. He", and Mg*\ and the temperature is now high enough 
to make these captures go at -i^ensible rate even with Hg*** to form 81* •. 

Once all the carbon has been consuMd, the central ten?^ 
star will rise agaln^due to gravitational contf action, and if the star is 
massive enough for this; temperature to reach » IO'ok, oxygen burning wlH 



I ,^ pJl ^ jjl . 

1 I ■ ■ 

. * I .. : . • . - 

I S'* + n 

' • ■ Si" 4 He" 



(20) 



Mg*** + He'* 

. ■' •■ ■ . 
This reaction series Is very much like that for C-burning. and again there . 

Is one direct neutron producing reaction. 

At temperatures, between 6 x lo" and 10'°K where carbon and oxygen . 
burning take place, the thermal photons 'hkve enough energy to cause dlsln- '. 
tegration of many of the nuclei that have been built up by the reactions 
described above. This photodlaintegration la analogous to thermal ionlxatlon 
of gases such as H which takes place at »10"ok'. Because nuclei are bound 

■ IT ■ . .' 



together much more tightly (ty a factor of « 10 •) than electrons are bound 
to atoms, it follows that nuclear photodislntegratlon should take ■place ftt 
much higher temperatures, of the order of 10*% For example, Ne" vhose 
abundai^e we have noted already poses a problem is further broken down by 
photodislntegratlon at temperatures above ID'oR: 

. Y + Ne" 0»« + He* ... (21) 

This Is Just the Inverse of the alpha capture by 0»%' and In general photo- 
dlslntegratlons are the Inverse of earlier bulldlng-up reactions that yield 
a gamma ray. I.e., the reactions A + a * (A + U) + y goes both ways at high 
\ enough temperatures. 

At still higher temperatures, "beyond about I.3 x lo' .°K where oxygen 
burning is complete, the nuclear reactions are mostly of a rearrangement 
type In which the existing nuclei are photodlslntegrated, and the' products 
of this disintegration are quickly captured- by another nucleus. In such 
reactions there is a trend toward equilibrium in which. mainly the. most 
stable nuclei; survive . There 13 a broad maximum in the nuclear stability 
curve, the curve of binding ener^, in the vicinity of Fe", so equilibrium, 
favors the formation of nuclei in the vicinity, of the iron peak. 

The rearrangement reactions begin to. be important , . 1 . e . , to occur ^ 
rapidly, at about 3;x,10' °K, and at these temperature/ Si M. ip the' pre- 
dominant nucleus remaining since,, it is more tightl/bound than the other ' 
products of oxygen burning which have largely been phptqdisintegrated by 
this time. Accordingly, the rearrangement reactions leading to iron peak 
nuclei are often considered to begin with a predpminantly Si" gas. A . 
great variety of reactions proceeds simultaneously with the result that the 
details become very difficult to follow. However,, the overall trend toward 



th« iron peak elements !• clear. These reactions tending tovard equlijt- 
briuai concentrations of the eleoents ar« referred to as e- proeessss in 
*the parlance of nucleosynthesis. / ' 

in order for equllibriua to be approached, It is necessai^for the 
rarloua rearrangement reactions to proceed faster than the eleoents so . 
fpr»d are able to undergo beta decay, and this is vhy the e- broceas re- 
quire, tenperaturea above 3 x 10 ' degrees. At slightly ItSwSr^ten^wratures 
the 0-decay processes are as rapid as the photodlsintegrations ind lead 
to other than an equilibrium mixture of nuclei. This non-equilibrium mix- 
tp-e is very difficult to calculate in contrast with the equilibrium abun- 
dmices which are easily calculated by .aethods of. statistical mechimlcs from 
relative binding energies of the nuclei .• These ; binding' energies are ac- 
curately known from laboratory measurements. ^ 

The calculation of equilibrium abundances followa: the same procedure 
as calculation of parent-daughter ratios in the secular equilibrium of . a 
natural radioactive decay chaiii. The method can be illustrated by a simple, 
three element chain: 

Jl) ■ ' - ' 

■ " .(2);, 
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Each line in this figure represents one element designated by number^ 1, 
2, 3. The differential, equation for the decay of element (1)' is 'simply: 



dN 



1' 



dt 



0 • 



where is its decay rate. For element (2) the equation la a little more 
conqpllcated since N2(t) la Increasing as a result of decftys of N^^ vtalch 
convert (l) Into (2), but Is also being depleted by decay Into (3). Thus 



dNg 



Finally, Is built up entirely by decay of (2), so 



dHg 

dT" ■ ^2''2 

Thade three coupled differential equatlona can be solved \sy rather alinpie 
methods but In a much more complex scheme such as occurs In a star In yhlch 
there az^ very many elements and several Interconnected buildup snd At0iy 
chadns the equations become formildably complex. In the equilibrium caae In 
vhlch all time variation^ dN^ have become zero, the equations reduce to a 

set of simultaneous algebraic equations rather than differential equations*/ 
In the above example, these reduce simply to Xj^N^- - XgNg ■ 0 or Ng/Nj^ « 'X£yx2* 
Thufl the equilibrium ratio of N2 to Nj^ is given simply by thfe ratio Xj^/;^^ 
of decay constants. Even with complex stellar buildup schemes the equilibrium 
ratios can be obtained €U3 ratios of reaction rates without solving the dif- f 
ferentlal equations. - . . . w 

The production of neutrons becomes copious in the temperature ran^e 
above 10 by processes already mentiojOLed, emd these are available to 
build the neutron-rich isotopes of elements lighter than the iron peak; 
They can also build elements heavier than the iron group clear up to Bj^^^*. 
The next heavier element; P^^^**, is unstable and decays by a-decay as • 
rapidly as it forms in the relatively slow neutron addition processes that 



go on in evolirtlonaxy stages of^^riflial stars. These* neutron capture pro-* ' 
cesses that take place In late' red giant stars have been temed s -processes 
(**slov" neutron processes). ^Detailed calculations based iq>on neutron cap- 
ture cross sections measured In the laboratory (e.g. ObJl Rldge> with 25 ke^i . 
neutrons haye. been rather successful in predicting observed relative abun*- 
dances of elements heavier than Fe*\ up to Tl^^^ and even beyond to Bl*** 
with somevhat less accuracy. 

At still higher temperatures (> 5 x 10* ^k) neutron fluxes are greatly 
Increased, and neutron capture rates greatly exceed 3-decay and a-decay ' 
rates permitting the buildup of unstable elements heavier than Bl,^®* all 
the way to Cf^'* and perhajps beyond, though such high atomic weight elements . 
are very unstable and decay quickly. Not only Is high temperature but also 
high density required to build elements above atomic weight 209. The key 
to production of these elements Is rapidity of formation to overcome the 
cbnpetitlbn of rapid decay rates, and we have "seen that the reaction rates 
are proportional to density so high densities are required. In practice, 
densities above 10* gms/cm^'are general^ needed for the rapid neutron 
capture processes to be successful. These are generally calle'^ r-^rocesses. 



and the ten5>erature8 and densities that make them possible as element builders 
are thought to occur only in 'the explosion of a supernova. Details of these 
catacysmic events will be treated in ^^other\lnitext of this series en- 
titled "Explosive Astrophysics". It will be 'sufficient here to, say that 
within the order of 15 minutes, during which 'a typical supernova is believed 
to collapse and explode, the conditions are right for the r- processes to 
take place, and in this short period the synthesis of all elements beyond 

2 0 9' 

Bi occurs. In addition significant amounts of many elements between Ge ^ 



and Bi arm also produced r-proceises. For axamplo, from 30 to B0% of 
the flolar system gold abundanbe is believed to have been produced by r- 
processes. This aeans that most of the jsold in j)ne*s 

dental fillings,^ or in the rings on one's finger was created in the interior 
of a supernova* Such explosions also solve the problem of hov, once a 

heavy element is created deep Vithin a star, it manages to escape into the 

• . •■ • ' ■ ^ ' 'f . ■ .. . 

interstellar regions of space. Therr it may later be reassembled Vy gravity 

into» another star to give spectroscopic evidence of its presence, this time 

in the outer layers of the star, or in the case of the Sun, in its surrotmd- 

ing planets that vere formed from the same reprocessed intetstellar matter. 
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STELLAR EVOLUTION - 

I . Introduction • 

A star loses energy, so it must change vith time. The equilibrium 
structure- of a star is uniquely determined by its mass and composition. 
As it radlaties avaiy energy, it must lose mass (E « mc*), and this effect 
alone might be expected to o^^use a change in the" star. " In very luminous 
stars mass loss both by. radiation and by stellar wind can have soAe notice- . 
able effects, but in stars like the sun the mass change is very wnall. - - 
The sun, for example, loses less ttfah 10-'^ of its mass in a billion years' 
throufjh radiaibion,^ hardly enough io produce an appreciable change in strucr 
ture. , , . 

In the energy producing core of a star a more subtle but more Influ- 
ential mass loss occurs.. There hydrogen, is being converted into heavier 

elements, He^ C, N^' 0, and this causes 'the molec\xlar veight of the core to 

... _ ■'(*•.,, ' . ■ . . 

increase. Thj.s increase vis substantial . Suppose a star begins at "zero 
age" on the main .sequence vlth a typical cbinposition for solar type stars:! 
X » .72, Y = .26, Z a .02 giving y w 0.608. After ^m^ient time has 

(1) This statement is known, as the Vogt-Russell theorem^^icf holds for nearly 
all types of stars. \ . - •% 

^2) The notation here is conventional in astrophysics. X and Y ar'e respec- 
tively the fractional compositions by weight of hydrogen and helium. 
Z is the weight fraction of all elements heavier than helium. Since . ' 
. these cure fractional' compositions , X + Y + Z » 1. 

The "molecular weight" of the gas, y, is the mean weight per particle. 
Since the gas in the interior of a star is fully ionized, the electrons 
freed by ionization must be counted among the peurticles in determining 
y. The light weight electrons dilute the "gas" inside the star, and its 
effective molecular weight is considerably smaller than the usual, values . 
we are. accustomed to in chemistry. Ionized hydrogen, for example, hais 
y ■ 0.5 Insteaid of its un-ionized value of 1 familiar in, chemistry. De- 
. . tails of the derivation of the formula for y of an. ionized gas can be 
^rtound in the minitext on Stellar Structure or in any standard text such 
as^M. Schwarzschild, "Structure and Evolution of the Stars", Dover, •(1966;)« 
The restili is . ' I . ' : 
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p&ssed (alO*® years for a at ar like the sxm) moat of the hydrogen in the 
core has b^en changed Into helium, and the core composition might typically 
he X » .02, Y » .95, Z ■ .b3%iying y » 1.303, more than double its initial " 
Talue. Such a change in must cause a large change in structure, and ob- 
▼ioualy the evolved structure is quite different from the initially homo-- ^ 
geneous one. The "old" star has an almost pure He core whereas it began 
as mostly H. Other consequences of th^ hydrogen burning follow j too, as " 
will- be- seen. v 

We shall trace the evolution of stars of soveral different masses and 
two widely different compositionn from the initial stages of formation from , 
Interstellar gas and dust, through the main stage of life spent on the Main 
Sequence' and on to the various fates that can eventually befall stars as. they . 
grow old and exhaust their central 'supplies of nuclear jruel, ultimately col-^ 
lapsing sunder the grip of gravity. ; 

A star^ travels through life by a succession of alternating energy pro- 
duction mechanisms. First , gravity causes contraction of a gassy,' dusty 
interstellar cloud, raising its centred density and temperature in the process 
until both are high enough to begin hydrogen fusion into helium. Then nuclear 
energy production takes over the Job of supplying the eriergy the star radiates; 
away. Since this energy source is both efficient and well supplied witl^.fuel, 
it operated for a relatively long time, but not forever. The time eventually 
^. comes when hydrogen in the core is used up, and the star turns again to gravity^ 
to supply its energy requirements. The core contracts, and in so doing it ^ 
heats up sufficiently to ignite fusion in a shell of H around It. This shell- 
burning source then' takes over energy production, and the star once again re- 
lies upon nuclear energy. Gravity continues to squeeze the constantly ffrxTWlng 
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helium core, and eventually, if the star is aaaBiye enough, this core bejsina 
to %um" He to form C** hy the triple alpha process: 
.3 He"*^ C^* 

In this my a new source of nuclear energy Is added. Eventually the He in 
thei core is all converted into carl>on, and this nuciltear energy source ceases. 
Gravity once more takes over, squeezing and heating the core further until. 
If the star Is massive enough, even the carbon may begin to fuse into several 
products such as Ne^^ + He*', Ma*^ ♦ H, and occasionally Mg*' + h, an Impox^ 
tant iBOurce of neutrons for building heavier elements. A helium burning 
shell will alio form before the core begins bumlhg carbon. 

If the star is not heavy enough to proceed to the next possible core 
burning stage, it will still have one or more sliell bimiing sotiroes, of energy* 
and gravity will continue to operate on the core until it becomes degenerate, 
whereupon contraiption virtually ceases . Degeneracy Is a consequence of the 
Pauli Exclusion Principle which governs "the occupation of available energy 
"states in a dense collection of interacting spin particles such as electrons. 
Ordinary metals exhibit the features of a degenerate "electron jgas" that are 
of Importance in degenerate iBtellar interiors, namely high thermal conductivity. 
' and low compressibility. These metallic characteristics of a degenerate core 
will b€l of importance later in understanding certain phases of stellar evolu- 
tion. ^ ' ■ ; ' *. o- ■. : 

Our knowledge of stellar evolution rests upon two foundations. First, 
there sre qbservations of stars from which such properties as mass surface 
temperature, surface compositioxi, ^radius, and rate of radiation of energy 
can be determined.; A great deal of effort and ingenuity by observational 
. astronomers has gone into the gathering and interpretation of the data oni'* 
yihieh our knowledge about the mass, Ixnninosity, tenperatuire, radius, and bom-* 
position (M, L, T, H, y) of the stars depends. Much of this obseryatiemal 



Inforaatlon is sunimarlzed graphically in the two chief diagrams of astronooyt ; 
the Hertzsprung-RusBell (H-R) diagram and the Masfl-Luminosity (M-L) diagram, 
which ahow what values of M, L, T, U are exhibited by real stars. 

The second foundation of stellar evolution rests upon the branch of 
astrophysics, that deals with stellar structure. It is the task of this icienee'' 
to reproduce by^ calculation from first principles stellar models that corres- 
pond in observable properties (M, L, T, R, p) with real stars. Sinbe real ; 
stars f'all'^ into well defined and often quite natrrow regions of the H-R and 
N-L diagrsma , it is evident that natxire has placed tight restrictions upon 
the range of possible structures of a. star « To put it another way, we might 
say that the laws of physics select only those stellar structures that con- 
form to observed stars, and it is the task of the science of stellar structiire 
to leaiarn how properly to apply the laws of physics to yield realistic models 
of stars:. \ [ 

\Among the great numbers of observable stars are to be found a wide gamut 
of the properties M, L, T, R, y plus one additional variable: x, or "age". 
That Is, we are surely observing stars with a great range of ages from very o / 
young, 'newly formed stars, perhaps less than 10^ yeariT old, to those more than 
10^^ years old, compeurable to the age of' the Universe. All of these stars lie 
in well defined regions of the H~R diagram, and it is one of ; the tasks of the 
study of stellar "evolution to discover how the observable properties of a star, 
change with time and how a star moves in the H^R dia^am as it ages. 

There is clearly a very close connection between the study of stellar 
structure and the study of stellar evolution. The procedure to be followed 
in discovering how s^tars evolve Is straightforward. One starts with a model 
of a real star lying somewhere on the H-R diagram, usucLLly on the main sequfnoe 
' as a beginning. Thel internal structure of this star together with a knowledge . 



of Its Wnoalty or rat« of en.rgjr loss p«n»ite caleulitlon. of the rat* at- 
vhlch the coBpoaltlon chMgea Inside the star where" nuclear trMsmutatlonB. 
convert one kind of riucleus into another (e.g. H He) and the rat^ at which 
gravity squeezes the interior during evolutionary phases in which gravita- 
tional energy production is significant. A new model may then be calculated 
based upon the altered con5)08ition, u, of the slightly, evolved star.^ This . 
process, is repeated many times to obtain a whole series, of models eyolvl.ng 
with time. The observabie parameters, L and T,* forj'this series will be re- 
presented by a track on the H-R diagram. These tri^cks are the principle 
means by which the results of evolutionary studies are presented, and they 
: form the pbint of contact between theory and observation: the tracka must 
lea* through regions in which real stars are found whose characteristics 
agree with those of the theoretical models.. The time variations., of internal 
structures of a series of evolved models give further -detailed information ; 
about the changes that take place within a star as- 11^ ages, and^a full appre- : 
ciation of stellar evolution can only be gained by a study of such models. 
We will therefore supplement discussion of evolutionary tracks on the H-R ^ 
diagram by description olf the corresponding internal changes that take plkce 
in aging stars. The descriptions of evolution will necessarily be qualitative, 
but this should not obscure the fact that they are based upon carefully cori- 
structed, quantitative models in which such properties as density, temperature,, 
etc., have been yorked out as functions of radius. . 

ir. Earliest 'Stai^es : Star Formation . 

The.beginninf?s and endings of .the stars remain the least understood 
phases of stellar evolution; middle life is the best known phase. Therefore, 
In discussing -the initial' stages in the formation of a star we are dialing 
with an area that remains sketchy and rather speculative. Most treatments of 



stellar evolution are content to begin with the obvious fact that stars 
do form and' to start dealing with them in the stages Just preceding arrival 
at the main sequence period where they remain for most of their life burning 
hydrogen into helium. i ''\ . ' 

We shall look back further than this, however ^ in an attempt to outline 
some of the main featxires' in the birth of a star from the gravitational con- 



traction of a large cloud of interstellar gas and dust. 



Astronomers have convincing evidence that star, formation was not a one- 
time phase in the evolution of the Universe but is a more or less continuous 
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process which is going on even today. They also observe that the new, young 
stars form only in regions of the Qalaxy that have high concentrations of . 
gas and dust. The gas L chiefly hydrogen, can be detected by radio "astronooor 

by means of a spectral' line at a frequency of 1U20 MH^ which lies in. the 

■ . - ■■(:■■■ ■ 
r^eadlly accessible mi(|rowave region at a wavelength of about 21 pm. The dust 

consists of microscopjlc solid particles generally thought to be in the size 

range of a few thousand &igstroms and to consist mainly of graphite and slli- 

I ■ . " . .■ > . 

cates, perhaps with l^hin coatings of ice. { ' ^ 

Interstellar gas makes up 98-99^ of *the matter from which stairs condense^ 
and its concentratio]|i is roughly 1 atom per cm'. Thus the matter needed to 
form a star of aboutj one solar mass {^2 x lO'^ gm) must come from an inter- ^ 
.stellar cloud whose ^initial radius' is of order 10^^ cm (10 light yelrs). The 
tismperature of this |gas has been estimated to be ^lOO^K (from its infrared 



emission, among other things), and its stability may be tested by inquiring 



iuiri 



whether its gravitational attraction is strong enotjgh to prevent its outer 



atoms, from evaporat 



ng away. The gravitational potential energy of an atom 



of mass m at the. edge of a cloud of mass M and radius R is where p is the 



gravitational constant (6.67 x 10* in cgs. units). If this J,s equal to 
the thermal energy of the peripheral atom, i kT (where k is the Bolt zmann 
constant), then the cloud will Just be stable. It will collapse if 

which is the bondition that the pull of gravity be stronger than the dissi-. 
I>ating tendency 01' thermal motions.'^' It is more convenient to put this in- 
equality in;to a form involving the cloud density > P> eliminating? R by the 
expression ^ ^ a jgj^ Thia gives 



3 1* X 10-^*t' gma/cm-' (l) 



or 



p > U X IO^^^T' fM«l gms/cm* for contraction (la) 



It follows from this relatiob that a cloud of 1 solar mass (M ■ M^) 

"1 

at 100^ k will be stable against thermal dissipation if 4ts initial density 
is greater than about U x 10"^* gm cm"*. How does this comi>are with the 
density of the in^.ersteilar cloud containing approximately one atom per cm^t 
If the atom is hydrogen, the cloud density will be (6.03 x 10^^)""^ or 
1.67 X 10""^^ gm cn*^, about a factor .10* too low to condense. Such a cloud 
will dissipate unless coBq>res8ed by some external force, and the compression^ 

if it occurs, must raise the density by a factor «2 x lO^^jLlthout increasing 

. . ■ . \ . h 

Its temperature. If T increases during compression', Eqh (l) tells us that 
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the compression factor must be even greater than 2 x 10^ but thiS' is likely 
to produce an even greater rise in temperature and is likely to be a losing 

... J . . ■ ' ^ 

proposition. It might, be argued that if the compression wex*e slow enough, 
the cloud might, be able to radiate away most of its excess thermal energy 



and remain near 100° k. : This procesB, hoverer, would require a long time, 
of the order of 10 million years according to calculations, and it is diffi- 
cult to laagine an effective compressive mechcuiism (remember the' compreaa ion 
factor must l>e i.lO*) that operates so slowly. \ 

We are forced to the conclusion that a star such as the Sun probably did 

. .. } ♦ 

not cohdanse neatly and efficiently from a cloud of about one solar mass. In- 
atoa<fu It noBt likely condensed from a much larger cloud together with a 
Shower of pther stars, probably with quite low efficiency leaving behind most 
of its pairent gas cloud ^uncondensed into stars. Observations on the sites of 
interstellar clouds that are fairly well defined in our part of the Galaxy 
Indicate a range of several tens to perhaps a hundred light years in si«e 

10*" cm) with hydrogen particle densities already up by a factor 10 over 
the average Interstellar value of 1 per cm' . Such clouds may contain several 
hundred solar masses and are often observed to be associated with groups of 
young stars in accordance with the suggestion that stars rarely if ever form 
singly. 

A 100 cloud at 100^ needs a density greater than U x 10"^^ gm cm"^* 
to be stable against theirmal dissipation, according to Eqn. (l),aad a cloud 
of 10^ would need p > x lO"^^. Compared with p « 1.6 x 10"^^ gm cm"' 
for a particle density of 10 atoms/cm', even clouds this large are unstable, 
but not ]>y so large a factor as before.' More modest compression by external 
forces (radiation pressxire from nearly stars, collisions with other gas clouds 

magnetic fields) riow becomes conceivable -as a means of raising the density to 
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the critical point that will initiate gravitational contraction and subsequent 
star formation. . 

Once gravitational contrMtton begins, all is not simple because the 
temperature tends to rise,/ and this rapidly increases the tendency of the 
cloud to- dissipate again (note the t' factor in Eqn. (1)). As long as the 



cloud remains transparent to its ovn thenaal radiation, there is a chance for^ 
it to radiate away energy fast enough to keep Itself cool enough to prevent 
the contraction from stopping short of star formation. The presence of dust 
grains is especicdly valuable in converting gas atom thermal energy into 
thermal radiation: the gas atoms collide yith the grains, heating them, and 
the grains radiate at long wavelengths into space like microscopic blackbody 
radiators/ The gas is quite transparent to the far infrared radiation of 
the grains as l9ng as temperatures remain low, and the whole mixture of gas 
and 4ust appears able to radiate well enoxigh 'to let contraction continue once 
it begins. Many details are still not clear and cannot be pursued here. 

As the cloud contracts it heats up, and it also spins faster since its 
angular momentum is conserved. Let us first consider the consequences of. 
angular momentum conservation becauise it poses a problem in star formation 
every bit as formidable as the contraction vs. thermal disslpaticm problem 
outlined above. Consider a gas cloud with a radius 10^® cm and a mass, of 
100 Me a 2.x 10^*, grams. If it is a uniform spherical cloud, its initial 
angular momentum can be calculated by elementary principles to be |. MR^ (i*^. 
moment of inertia) times its initial angular velocity w^. 

The cloud will have, on the average, an initial angular velocity imparted 
by the shear forces that arise from differential rotation about the center of 
the GcLLaxy.^ ^Is rotation of the cloud arises from the fact that .parts of the\ 
* cloud closer to the center of the Galaxy haVe higher orbital speeds than 
those fcurbher away (Kepler's 3rd Law siaiys that the orbital period varies as 
the 3/2 power of the mean radius of the orbit). The figure below shows how 
differential rotation arises from the fact that to > u)''. o 



»-> w-. Thii phenomenon 18 called rf«r<rMHfl/ rofflf ton. pcn,.ncrr. 

Figure 1 ■ 

In. oiir region of the Galaxy the differential angular velocity, 
(Oq « 3 X ip""** sec""*; therefore a cloud of the size under consideration would 
, have angular momentum — MR^ (Oq « 2 x 10*^ gm Cn^ b«c** • Compare this with 
the angular momentum of the solar system, approxiioately 2 x 10^^ gm cm^ aec**^/ 
The .tvo values differ^ a factor 10^ when allowance is .made for the fact that 
the cloud contains 100 solar masses while the Sun contains only one. In short 
the 100 cloud must lose all but io*^ of its Initial angidar momentum if It ■ 
is to give birth to 100 well behaved stars. Put another ?raar, conservation of 
angulflor momentum requires that R^u) remain constant as the cloud contracts. A 
factor Ru) of this quantity is the equatoral,, surface velocity of the rotating 
cloud; so Rv. must be constant. Initially v„ 3 x 10^ cm/sec when R « 10*® ci 

By the time R has decreased by a factor of 10*, to 10*** cm, v has iOLrlsen to 

s ^ „ 

the speed of light! And 10*** cm is still quite large, about 1000 times the 

• ■ ■> ■ ■ ■ .■ ' ■ ■ . ' ■ ■ . ■ ■ ■ 

radius of a star like the Sun. 
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Clearly the contracting cloud must reduce its angular momentum somehow if 
it is to siicceed in forming stars. ^ At present ve don't know hov it accom- 
plishes this. Although there have been many speculative suggestions ov^r 

the years » noncj has been supported by fully ctJnvlncing arguments. One of 

■ ■ ■ * . "\ ' * • • •• 

the more plausible suggestions is that the cloud forms binary stars prefer- 

entially, and these are able to contain much more angiilar momentum than 
single stars. Observations confirm that a majority of stars (« 60%) do indeed 
belong to binary or mxaltiple star systems , so formation of binaries does seem'*, 
to be the norm. ' 

In order to- discuss the relative angularnnomenta of different configur- 
ations, it is convenient to deal with emgular momentum per unit mass which 
will be designated, L. The primordial clpud we have been discussing has 

; L a 10^* cm^ iaec*"* while the solar rstem haslT^lO^^ cm^ seC"* . A binary^' 
system consisting to Two solar mass stars in circular orbits at a distance 
r from each other has *L = i 7rr^/2P where P is the orbital period of the 
binary pair. Kepler's 3rd Law provides another relation, Mr? ■ Kr', between 
P and r,and this leads to L « P^'* or L « r^ '^ for a binary. \ Putt inj^^ in the 
constant K and the other numerical factors and inserting a sliati stickily 
weightel^ average period for observed binaries, ve find that L « 2 x'lO** cm* sec 
for typical binary stGurs^ a value about 200 timed' greater-than |ihat for an ' r ' 
average single star-. From this kind of reasoning we see that binary/ formation 

.vis a much more efficient mechanism for storing angular momentum thantsingle . 
star formation. Even so, the primordial cloud xinder consideration has^enough 

T to supply « lO' binaries, far more than we believe actually form from such - 

' ■ ■ , ^, ■ ' ' \ ' ■ 

a cloud Fon the basis of observational evidence. If only a small fraction 

of the cloud mass actually condenses into stars we^rban account for observed 

* ■ • . j- ■ ■ • ■ ^ ^ . " " 

rates of new star formation, euid' the cloud retains most of its initial angular 

11 



momentum. Details of the fragmentation process and .what happens to halt 
conti»actlon of the whole cloud once star formation begins remain obscure. 

There Is a third Ingredient to be taken Into account^ In the early stages 
of star formation, the interstellar or galactic magnetic field. v^ls field 
la very small (of the order of a few times 10- ^ gauss) in the spiral arms of 
the Oality -where stars are formed. Compare thia galactic field with the 
Earth^s magnetic field of « 0.6 gauss to appreciate its small size. However, 
even so small a field ,^ acting as it does over vast distances, can have a 
very strong Influence upon certain stages of star condensati^on. As long as 
the jLpterBtellar gas is cool enough to remain neutral , the fl^^does not 
interest directly with it. Once the gas becomes ionized, however i the fields 
and the charged particles interact strongly, and the motions of the ionized 
gas are greatiy influenced by the prVsence oif^t he field. Likewise the field 
lines are Influenced by charged particle motionia and the whole M«plex-of 

field plus ionized gas becomes an interdependent entity, a magnetohydrpdynwiic 

...■0' --^ 

plasma. The properties of such plasmas cannot be discussed in any detail 
here , but "we can illustrate in a slnple way one serious problem that comes 
from the presence of the magnetic field. 

ConsMer a protostar of about 10*^ cm radius with sufficient density to 
be^able to contract to form a real star of radius about 10|^ cm, a realistic , 
value (the Sun's radius is « 7 x 10^* cm). If the protqptar of 10*^ cm radius 
has an appreciable amount of ionized hydrogen, as calciilatll'bna' indicate it will, 
these ions are >ouad to the magnetic. field lines. Sufficient ions will be 
fo»ed in the cloud l^ .-osmic rays and by ultraviolet radiation from nearby 
hot ^ blue "stars . Througlix the mechanism of collisions between ions and neutral 
atoms, the latter are indirectly connected with the field. Unless the frac- 



tlon of gas that is ionized is extremely small, this indirect connection 



between gravity and the magnetic field through the intermediary of neutral 
particles colliding with ions is strong enough to effectively tie gravita- 
tional contraction to the magnetic field. As the gas contracts it therefore 
co«presses the field in the plane perpendicular to the fields as illustrated 
in the^ sketch belov: . ' 




(a) Initial gas cloud 



Figure 2. ,Ck»preaBion of magnetic field lines in a contracting partially 

ionized gas cloud. / 

... . . . v ^ . / - 

Conservation of magnetic flux in such a contraction leads to the con- / 

■ ' • • \ ^ -. K • 

elusion that the product BR* is conserved vhere R is thexradlus of the cloud. 

Thus B ^.R^^pSnd in going from R « l^J^cm to R 10^/ cm, B increases fncm 

initial values of a few times IQ-^ gauss to well over 10* gauss. ^ Such /high . 

\ . . .' • 

fields have never been found in norWl stars. The strongest fields obsexnred 

have, been fo\md in. certain i>eculiar AUtype stars with values of B to «30,p00 

gauss. Most stars like the Sun have rields ranging from a few gauss to a few 

hundred gauss. Clearly the star finds some way to rid itself /o^^^^ 

flux lines that pass through its parent gas. cloud, but the means by. wh|ich it 

accomplishes this are still not understood. 

We have gone to some lengths to illustx^ate serious difficulties that 

hamper our present knowledge of the earliest stages of star formation. The 
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fact that these difficulties can be presented in terms of elementary phy»ic«J- - 
principles makes it appropriate to discuss them in Aa introduction to stellar 
evolution such as this. But whatever means stars employ to overcoBse the 
barriers to their formation, they somehow manage to do so, and once they 
reach a certain stage of condensation, (radius « 10^' cm, density p » lO^^gm 
cm-', temperature a lO'^K) we can folloy the evolutionary history of the star 
with much greater confidence except during a few short lived stages late in 

;life.' ■ ' . , • 

' Let us first consider paths on the H-R diagram of stars from the state tj^*^ 
described above until they reach the Main Sequence. As long as the iprotostar 
obtains its radiant energy from gravitational contraction (before nuclear 
energy sources belgin to supply its energy when it reaches the Main Sequence), 
there is a very useful general physical theorem, the Vi rial Theorem, that Is 
very helpful in xinderstanding the sta;r'8 development. This theorem comes 
from considerations of the overall effect of converting gravitational potentia?. 
energy, -V, into kinetic or thermal energy, K.E. The theorem says that in such 
conversions 
\ A K.E. » ilAVl 

\ • • ^ . ■ ' • ■ ■ ' • ' 

That is, half the potential energy released by contraction goes into kinetic 

energy of thermal motion (i.e. into raising the temperature of the gas) ;- the ^ 

other half leaves the star as radiation. This simple relation permits us J 

to follow the average temperature and luminosity of the contracting star as 

long as it obeys the conditions for validity of the Virifljl Theorem. These 

conditions apply as long as the star is contracting slowly ^^ough to maintain 

approximate hydrostatic equilibrium in which gas pressure balances the force 

■ ' ..... 

of gravity • o 



There are phases^ of contraction during which the Virial Theorem does , 
not hold. The first of these is encountered when the mean temperature of 
the gets reaches approximately 2000OK and the molecular hydrogen present 
dissociates into neutral atomic hydrogen. This is an endothermic reaction 
that absorbs thermal energy, convi^rtlng it into potential energy of dissocla- 

tion. That is, 211^-:; energy, vnien disBociation occurs at a rapid rate; 

' ' ■ ■ ' ''^ ' ■ ■ . . • . " . •■ 

the star undergoes rapid contraction to supply the required* dissociation 

* \ ' , ' ' " 

energy. Its temperature remainn nearly "constant , but its radius decreases 
and its luminosity must drop, (too, ain^e L ^UttR^oT**. This phase, therefore, 

' ' • ■ ■ . . ■ \ ' • ■ • ' ■ 

appeeurS'-on the H~H diagram as a nearly vertical drop in L at nearly constant 
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■ . ^ ' ' ■ ■ ■'^^'^w , ■■ ■ . /■ \-V. 

Figure . 3 • H-K p:.a£ram £ho>vinr;; r.:<> ^\ gravitational .contraction iptagte of \ ■ \ 
star formation v 



This diasociation phase of the collapse begins at R =* 10^* cm and is com- 
plete when tlie radius has decreased to about 1000 solar radii (« 7 x 1,0*' cm) 
Figure 3 does not attempt to trdce the collapse betveen the initial cloud 
(R « 10*' cm, T =»>100°K) dovn to I^oint A. This is a dynamic f>^^}e-fall pro- 

■ •■ . A ■ : ■ ■ . 

eeas« and t^e cloud is transparent to the long vavelength infrared radiation _ 
(mainly at 58 microns )\ by vhich it radiates its excess energjr. Thio phwe; 
is characterized as free-fall because, the cloud Is so tenuous that there la 
insufficient pressure, to halt or even aloy appreciably the gravitational pull 
on individual gas moleculea and atoms. These fall initially inward toward 
the center of the cloud altnoVt as if there were no other particles present to 
• Impede their fall. Th^ time rfl^quired for a particle to fall to the center 
.xuider such idealized conditions is csLlculable by elementary methods; It is 
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The free-fall time for the initial gas\ cloud turns out to be T x 10* years, 

■■. . .- . . ■■■ ■ r ■ ■ -x. " ; ■ ■ ■ . 

assuming its ma£5s is initially - 100 ^before It fragments Into s61ar mass * 
size stars sometin^ during this phase • ' , . - 

After dissociation is complete at C« the gas temperature can rise' again 
as gravity contracts the star farther; ^e rise in temperature is aided oon- 
siderably by the fact that it is nov opaque; in fact calculations show that 
it Just becomes opaque fi)r. radiation originally at the center at A Just prior 
to dissociation. This helps the star trap its radiation Inside temporarily 
and speeds up the dissociaiibn process. The Interior temperature continuee . 
to rise rapidly to « goOO^K where le^izatlbn of atomic hydrogen takes place* 
•This too is an endothermic reaction, and lAiile it operates « the star la far. , 
from being in eqalltbrlum. , It continues to shrink rapidly along the dTftflhed 
line (dashed because of uncertainty abput details in this phajEie) from C to D. 



-The radius decreases agaiCi by about a- factqr 20, but this time the surface 
temperature of the star (the abscissa of the Figure) does not remain con- 
stant as it did before • In going from'A to B the star vas not yet fully 
opaque. At first all radiation except that originating at the. center could 
escape vithout being absorbed at least once (Point A) ; later between A and B> 
as the density increased, radiation originating farther and farther from the 
center vas Just able to escape without absorption. Only at t^he end , at 
had the opacity risen to the point that all radiation except that; from a 
^ thin layer at the sixrface was absorbed before escaping. .Actually, at C and 
beyond most radiation in the interior is absorbed and reemittcd many times' \ 
liy atoms in the Btellci- interior before it finally reaches the surface. This 
takes some time ^ and the siirface radiation no longer reflects the instantaneous 
state of the interior as it did earlier when the star was still, someirtiat trans- 
parent, v The net result is that the interior heats up very rapidly at B to the 
gOOO^K or so required to ionize hydrogen, but the surface teniperattire plotted 
on the H-R diagram lags in time, rising until it finally reaches - gOOO^K at 
Point D. The whole process of dissociation, ionization, and final buildup of 
surface temperature, the trip from C to D, only takes about 20 years for a 
solar mass, star , and this must be reckoned as virtually instantaneous on the 
Udual time scale for stellar processes. 

The little' U-shaped excursion at the end of the dashed p^th just before 
Point D is the result of a minor implosion of sorts. The star collapses so 
rapidly that am inward radial shock wave forms, and when it reaches i^the now 
rather dense core near the center of the star, it is reflected back out, a 
process vividly doscribed as ,,a "core-bounce". The shock wave travelling back 
out reaches the surface where its higher temperatiure becomes evident for a 
short time as a transient temperature and luminosity spike as shown adjacent 



to 0. This transient flelre up laisTts about three months In a- solar oaas 

■ - . ' ■ ■ ... ■■ • • . • •;. 

star. • ; , . ■ * . ^ 

At JPolnt b the Interior of the star is turbiilehi as a result of the " 
rapid heating 9 the shock vaTe8> and the generally non-*equilibrlum nature of 
Its prepedirig collapse sttiges. Under these conditions It transports energy 
from the Interior to tlief' sixrface mainly by convection. This is, a relatively ' 
efficient form of heat transport » and It serves to stabilize the surface 
teinperature of the steur as it continues to contract toward its eventual else 
that it irill reach 'on tbfi Main Sequence* Curing Its contraction along D to 
e and f at nearly constant the interior temperature continues to rise 
rapidly, especially in the central regions. The convective^ turbiilent condi- 
tions in the interior insure ef-fective mixing of the material of the star so 
vhen^he central temperature eventually reaches about 2 x 10* ^ at E> any 
deuterium » lithium > berylliusn, boron, C^*, M*', etc are capable of under- 
going nuclear reactions that, destroy thea« Because of the cohvectlve mixing 
bf. gas into the hot core in this phase, virtually all of these esoteric light 

nuclei in the star are "burned- before the star reaches the Main Sequence , and 

* . " ■ ' * • . 

» « . * . ifrw • ■ * ' ' ' , 

dbservations conf ir6 abnormallyAConcentrations of these nuclei in normal stars. 

As the star contracts from 0 to E along the liearly vertical track (the 

"Hayashi trafck" as it is called after its discoverer) it remains in quasi-* 

hydro8l;atic equilibrium, a condition a)nenable» to study by standard methods of 

. calCTilatlon of stellar interiors. The. star has finally reached a vell-behaved 

condition th^t permits ^^some Insight into its internal structure, and a \. 

number of researchers have calculated detailed series of evolutionary models 

*of stara in this nre*Main Sequence, phase-. These calculations show, for example, • 

that* at Point D the star begins to depart from a condition of complete turbulance 



and a "radiative core'' develops ^ so-called because energy transport vlthin 
it is by means of radiative transfer. That is, energy moves outward in the 
form of radiation that is successively absorbed and reemitted by the gas 
particl^ss, nov almost entirely ions and electrons at the* temperatures^ of the 

core. As time goes on the radiative core grows at the expense of the con- 

.' ' . ■ . . ■ • * • . ' • ^ ' ■ • 

v^ctivis region, and the boundary between the .two gradually receeds toward 
the surface. Radiative transfer is a quieter but less efficient means of 
energy transport than convection, and it drastically changes the interncLL 
character of the star. 

Figure U, gives the results of a number of pre-MS evolutlonairy tracks 
calculated by Iben, and Table I following it lists the times in yeeurs required 
to reaeh the corresponding points on Figure U beginning at Point D in Figure 

3, t^e Btai^ of the Hayashl track-. 

/ 
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Fia. 4 Pre-MS evohitioiiary tracks on the H-R diagram for models 
of maL M (solar wits) ^ 0.5. 1.0. 1.25/1.5. 2.25. 3.0. 5.0. 9.0. and 110 
(indicated alongside the appropriate tracks). Numbered pomts alon^.the 
tracks are the point? reftrred to in Table I , 
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1 



kKjiptmuniu) 



Point 


15.0 


9.0 5.0 


3.0 


12S 


1 


6.740(2) 


1.443(3) 2.936(4) 


3.420(4) 


7.862(4) 


.2 


3.766(3) 


1.473(4) 1.069(5) 


1078(S) 


5.940(5) 


3 


9.3S0(3) 


3.645(4) 2.001(5) 


7.633(^. 


1.883(6) 


4 


2.203(4) 


6.987(4) 2.860(5) 


1.135(6) 


1503(6) 


5. 


^657(4) 


7.922(4) .3.137(3) 


1.250(6) 


2.818(6) 


6 


'3.984(4) 


.1.019(5) ' 3.880(5) 


1.465(0 


3.319(6): 


7 


4.585(4) ' 


1.195(5) 4.559(^ 


1.741(fi) 


3.993(6) 


8 


6.170(4) 


1.505(5) 5.759(5) 




5.855(6) 






. ■ ' Af 






Point 


1.5 


1.1S 


1.0 ~ 


0.5 


1 


2.347(5) 


4.508(5) 


M89(5) 


3.195(5) 


2 


2.363(6) 


3.957(6) 


1.0S8(Q 


1.786(6) 


3 


5.801(6) 


8.800(6) 


8.910(6) 


8.711(6) 


4 


7.584(6) 


1.155(7) 


.1.821(7) 


3.092(7) 


5 


8.620(6) 


1.404(7) 


1529(7) 


1.550(8} 


6 


1.043(7) 


1.755(7) 


3.418(7) 




7 


1.339(7) 


2.796(7) 


5.016(7) 




8 


" 1.821(7) 


1954(7) 







t Numbers in parentheses are tbo powers; of 10 by which fbe corresponding entries ara to 
be multiplied. 
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Table I 



The general time scales for Hayashl tracks ..(D* B) of stars of various 
masses should be noted. The Sun, for example, requires about io** years to 
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reach Point in FJ^gur^ 3 (Point 2 in Figure k) where development of a 
radiative core begins and anpther 8 x 10^- yeara to reach E; where this develop*^ 
sent* is complete and the star turns to the left (points 3 to 5 in Figxire h) 
.ailffloat parallel to the M.S. The same sequence in a 3 star, on the other . 
hand, takes only « 8 x 10 years, a factor 10 less'« This illustrates clearly 
how rapidly the evolutionary time scale speeds up as one considers more mas- 
sive stars. Br contrast a 0*5 star continues all the way to the M.S. on 
its Hayashi track taking - 10^ years to do so. A star more massive than 
a 3 hiairdly has emy Hayashi track at all, and massive stars scarcely appear 
to experience any break in the leftward trend of their tracks beyond the onset 
of opacity. The rapid pre-MS evolution of heavy stars is noteworthy.'' For 
example a 13 runs through the entire gamut, 1 through 8, of Figure U in 
only 60,000 years, a factor lo' less than a solar mass star. 

The mass-luminosity relation gives a good approxisiate idea about relative 
time scales . in stellar evolution. Both theory and observation reveal that 
luminosity L is proportional to some power of the mass of ' a star. This can 
be written 

■ i ■ ■ ■' 

where n varies somewhat with a star's composition and position on the H-R 
diagram, ranging from ~ 2*3 for very low mass stars to over 5 foi* massive 
stars. The L-M relation can then be written in terms of solar, values as 

■ • . ■ ■ ■ ■ ■ . • ■ \ 

k ^ (VI f \ ' 

- , . ■ . \ ■ ■' ■ 

The luminosity of a star is its rate of energy loss and the time siient by a 

star dvtrlng a certain phase of its evolution should be inversely proportional 

to the rate at which it loses energy (i.e., T « L~M. This, taken together 

with the mass-luminosity relation leads to the expectation that the time 



22 



Boales for Btars to run through corresponding stages of evolution should 
follow the relation 

■• • / ■ -.v.: 

fiacamlnaiion of Table I reveals that n remges f rom « 1.5 to ■ 2.5 depending* 
somavhat' upon mass but laore upon the stage of evolution considered, n « 2,5 
lava fair average for the trip from 1 to the Main Sequence in Figure U, 

The ialightly rising, leftward movement of a star in Figure k representa 

• ' . • * *• ■ ■ ■ ' • *. 

a stisige before the onset of nuclear energy production in which the star still. 

obtains its energy by gravitational contraction with a hot^ compact;, Tadia:tlve 

core. When thia core becomes hot enough, at point 5, nuclear reactions iura 

ignited Bind begin converting hydrogen into helium plus considerable energy. 

The detailed reactions are fairly complex and are diacussed elaewhere in the 

mini-texts on Stellar Structure and Nucleosynthesis. The overall result may . 

be summarized 

k H^-** He* > energy. 
The energy comes from the fact that a He** nucleus is slightly less massive 
than U nuclei. The mass difference appears as enr^gy according, to the 
Einatein relation € « Amc^ where Am ia the castt difference. This energy 
appears initisLlly mainly in the form of ijehima rays a^r' ^constitutes about 0*7% 
of the total mass-energy of the reactir^g nuclei T 6 the process is leas 
than 1% efficient in converting mass into energy, hut that low efficiency is 
adeguite to keep the star going, for a ver/ long rtime, over 10 billion years 
for a solar type star. ; o 

As soon as nuclear fusion taken over th6 Job of. furnlching energy, 1;he 
starts interior readjusts itself, contracts slightly ^ and its luminoaity 



deereaaes a bit as it aettles down onto tha.Main Sequence to spend the major 
phase of its life quietly converting hydrogen into helium in its central re- 
gions, The last points shOTO for each modal star in Figure U (point 8 in 
noat cases) represent arrd>^ral on the MS. By custom, arrival on the MS is ' 
reckoned as the true blfkh^of the star and at this point it is considered a 
Zero Age Main Sequenc* star. The star' is of. course observable throughout 
much of its pre-MS gestation period where . its lumihoaity in actually higher 
than the- MS star it will become. Observationally the T Tauri and Herbig 
emission classes of stars are believed to consist of sta'^in the pre-M8 
region, and the Herbig-Haro objects may well be still earlier pre-MS stars 

In the prp-dlBBOciation stage of Figure 3. Thus atars become visible well 

• ■ ■ ■ . . .' t ■ 

before their conventional birth, but at ZAMS they become much easier to. obr 

serve because they spend so much time on the MS that there are a great many 

there at any given time for us to see. By contrast, stars move rapidly 

through many phases of their pre-MS period arid consequently at any given 

time there are few present In these phases to be obsejrved. 

• * * . ■ 

III. Middle Sta g es; Main Sequence 

Stars on the MS obtain their energy by two principle series of fusion 
reactions. Both have the same net result, namely fusion of 1* protons into 
one helium nucleur*. Biit one, series of reactions, the proton-proton chain 
Involves auly light eieiAants: He**, He*, D^; h\ and it operates effectively 
at relatively ''low" tCTiperatures (low for stellar interiors) biatween «5 to 
25 million defJT^eej?. rj-ie energy generation rate for this procooa Js roughly 
proportional to the fourth power of the temperature near the center of the 
star. The other Droceos involves the normal isotope of carbon, C*^, in a 

■ ■ V . "214 ^ ^ " 
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complex series of reactions that eventually regenerates the C^* again, so 
carbon acts as a catalyst in the process. Various-isotopes of nitrogen 
and oxygon are also involved and the proceisa is known variously as the carbon 
chain, the CH chain, or the CHO bi-cycle. It begins operating at about 
12 million degrees and continues up to a.50 x lO^^K. This cycle is therefore 
effective over a somewhat higher temperature range than the p-p chain, and 
it is characterized by a much stronger teiaperature dependence, =« T*^ at the 
low end of the range dropping to » T^** near the upper end.. 

The large temperature dependence of the carbon cycle has an Important 
consequence. It leads to a turbulent, convective core in stars in which 
this mode of energy production dominates, i.e. in stars a little more massive 
than the Sun with central temperatures above « 15 million degrees. These 
more massive stars that operate mainly on the ^carbon cycle with convective 
cores and radiative envelopes are designated Upper Main Sequence (UMS) s;tar8. 
The less massive stars which operate on the proton-proton chain with radia- 
tive cores and convective outer layers are designated Lower Main Sequence 
(LMS) stars. The difference in structure between these two broad types of 
MS stars has important .consequences in the manner in which they evolve, as 
we shall see. 

Once a star reaches the Zero Age Main Sequence (ZAMS) , it begins in 
emest to convert Its hydrogen core into helium. Only the central 15-20^ 
of the star is hot enough and dense enough to sustain hydrogen fusion, and 
it is only in this core that the chemical composition changes while the star 
is on the M.S. If the star is less massive than « 1.1 M^, the core is quiet, 
non-oonvective, arid non-turbulent because of the relatively low T**) 
dependence of the energy generation rate. The temperature rises to a ioaxloum 



at the center of such a core, and so does the rate of energy, generation. 

■ ■ . - ■ ■ . ■ ' ■ \ 

Since there is ho convective turbulence, th«re is no mixing going on in the 

core and the fusion product , heliiun, stays where it is formed. It^orms 

BK>8t rapidly at the center vhere the rate of fusion is highest, and In time 

there develops a amail but growing subcore of pure He.. - Nuclear energy pro- 

duction ceases in this aubcore since it is not hot enough to fuse He into\^. 

It remains inert, the "ashes" of hydrogen burning, while hydrogen fusion 

continues in a thick shell outside the He subcore. This He core continues 

to grow, of course, as the hydrogen burning shell. makes more helium which is 

added to the ever-groWing He subcore. This pure He core remains at constant 

temperature rather like a loaf of bread in an "oven" of hydrogen burning 

shell. The shell continues to bum outward as a grass fire burns in a thin, 

* - ■ 

ever widening circle leaving a growing, charred central area. In this faaion 
the He core grows until its mass becomes i^about 10% of the total mass of the 
star. Gravity begins to squeeze the cor4 raising its temperature signifi- 
cantly so it is no longer isothermal. Eventually its central temperature 
will reach « 100 million degrees which\is li'^gh enough to initiate helium 
fusion by the "triple-alpha" process , 3^He" y C^^ This process generates 
nuclear energy, and the star now has two nuclear energy sources, a He burning 
core and a hydrogen burning shell. All of thia^ is still deep within the star 
but the effect is felt at the surface of the star where the additonal energy 
generated by He fusion increases both the luminosity and the size of the 
star while at the same time lowering its surface temperature. The star at 
this stage in its evolution has therefore left the main sequence and become 

a red giant . . / 

Figure 5 shows i>ost-M3 evolutionary tracks of a series of stars ranging 

from 1 to 13 solar masses and Table II shows the time intervals between nun*- 

bered points on the tracks. Point (1) is the ZAMS point for each star. The 

stage (l) (2) corresponds with exhaustion ot H in the core. The stage 
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lation I" composition (X — 0.708, Z ■= O.Cf2 for all track* except the 
A# - 30 track; for this track, 0.70, Z = 0.03). 
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Table II. EvoLunoNARY LiFEnMES (Years)! 

(Inisiai composition: A'- 0.7M/Z- 0.02 /or 1.0 :S Af ^|^15.b;^A'« 6.70, 
Z - 0.03 /or Af- 30.0) 

Af (solartmits) 



Point 


1.0 


1.5 


3.0 


5.0 


9.0 


15.0 


30.0 


1 


5.016(7) 


1.821(7) 


'2.510(6) 


5.760(.5) 


1.511(5) 


6.160(4) 


2 (4) 


2 


8.060(9) 


1:567(9) 


2.273(8) 


6.549(7) 


2.129(7) 


1.023(7) 


4.82(6) 


3 


9.705(9) 


1.652(9) 


2.394(8) 


6.823(7) 


2.190(7) 


1.048(7) 


4.91(6) 


4 


1.0236(10) 


2.036(9) 


2.478(8) 


7.019(7) 


2.208(7) 


1.050(7) 


4.92(6) 


5 


1.0446(10) 


2.105(9) 


2.488(8) 


7.035(7) 


2.213(7) 


1.149(7) 


4.93(6) 


x6 


1.0875(10) 


2.263(9) 


2.531(8) 


7.084(7) 


2.214(7) 


1.196(7) 


5.45(6) 


7 






2.887(8) 


7.844(7) 


2.273(7) 


1.210(7) 


5.46(6) 


8\ 






3.095(8) 


8.524(7) 


2.315(7) 


1.213(7) 




9 \ 






3.262(8) 


8.782(7) 


2.574(7) 


1.214(7) 




10 










2.623(7) 







* Numben In porontham iri the powon often by which the lorroipondlng enirlon nra to 
b« multiplied, . c . 
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(2) (3) corresponds with development of a hydrogen burning shell source 
outside of the isothenuaLL He core » and at Point (3) the core has' reached 
the limiting mass (« 10^ of the mass of the star) beyond vhioh the core 
begins to heat rapidly, under gravitational contraction, developing a hot 
center that iiltlaately ignites He fusion. 

It is evident from Figure 5 that there is a qualitative difference " 
bttveen the stapes of the tracks (ij (2) (3) in a low mass star (i 1.1 M^) 

and tracks of more laassive stars. The low mass stars have non~convectlve 

» . « ■ •<i ■ 

cores and \develop H-humlng shells graiiually and without fuss. This produces 
a i^mpoth track without breaks or kinks as the limiting isothermal^ core mass 
is reached at (3). It takes a 1 star about 9 billion years to traverse . 
the track from (l) to (3)> as Table II shows. 

Jt is easy to estimate the time a normal star spends on the Medn Sequence. 
Let its total mass be M gra and its luminosity be L ergs/sec. The efficiency . 
of ^hydrogen bunilng is .007 meaning that 1 graa of hydrogen . converted into 
He gives •OO? c^ ergs of energy . If the star stays on -the Main Sequence 
ixntil 10)( of its mass has been convert^dd into He» the total energy released 

is 7 X lO-** Mc^ » 6.3 X 10*^ M ergs./ If the ateyr radiates energy at the rate 

■'. V - " • , 

L ergs/nec for a time T on the M.S. .'then 

' ■ • ' J ' .' ' • ' ■ ■ . 

"^as ^^-S * ^0*^(^J"cond8 

Applied to the Sun (L =« U x 10*' ergs/sec, M » 2 x lo'^ gm) one gets X ^ 

3 X 10^^ sec » 10^° years which is iust slightly more than the 9 billion years 

calculated for the Sun to evolve froia points (l) to (3) in Figure 5. For 

other stars the equation above can be put into the form of a ratio to solar 

. • ■■ . f ■ . '■ ' . ^ . . * 

i^ues:/ 



. If la posBible to go one step further by using an • .^r-Xi -iate mass- 
luminosity relation, L « vhich is' a reasonaMy good approximation over 
a large part of the Main Sequence. With this relation, the lifetima of a 



star on the M.S. becomes 

years. 



We .should, expeci a 9 M, star to have « i.U x lo' years . -which is ' not far from- 

\ ■■ - 9 ' ■ . . ■ ■ - • • 

value of-2.1 X lO' years for a star' of this mass from Table IIv It, appears, 
therefore, that this very simple theory 'gives a good account of the time ■> 
spent by a star on tlie Main Sequence. . ' 

Heavier stars than about 1.5 M^ have turb]ulent, cohvective cores initially 
since they operate more or less, on the CNO cycle with a yery- sensitive temper- 
, ature dependence = T*". The core turbulence keeps the core well mixed and ^ 
prevents the smooth growth of an inert He subcore beginning at the center; 
Instead,, ao H-buming goes on, the entire core (15-20^ of the .stellar, radius) ; 
Tjniformly changes composition, becoming m6re ^djnore rich in He as the- H con- 
tent declines. This process must eventually result in a core too poor in H 
to burn any more; the ashes have diluted the fuel till the fire goes out, or 
nearly so. At Point (2) in Figure 5 such a situation has occurred. Ae' the 
nuclear energy source begins to fail, gravity steps in to take over energy 
production by contracting the core. This raises the temperature of the core 
and manifests itself by an Increase in both luminosity and surface temperature 
of the star, causing it to move from (2) to (3) as-sfiown in Figure 5. During 
this stage of evolution a massive star contracts slightly causing the track 
(2) (3) to lie In the direction shown (tqward higher T, L). The sharp kink', 
at (2) reflects the relative suddenness with wliich the entire core runs out - 
■ ' ■ . ' • . ■ ■ \ ' ■ ■ - • ■ ' ' \ 



of hydrogen. The lov mass stars have no such sudden fuel exhaustion and 

1 ■ ■ ■ ■ " • • . ■ . ' ' 

their •cores remain isothermal between (2) and (3) vhile the H-burning shelle . 

. •• \ ■• ' > . ' ' . ■ •'- 

oontinu^e to supply nuclear energy As a consequence, the radius grows \ . 

. ' • • ■ ■ ' I" ■ ' . ' • 

smoothljT in a less massive star, its luminosity increases, and its surface 

temperature &tays nearly constant, as the Figure shows for the case of a 

. \ y\ \ ■ ■ , • . . ' • / ^ . '-rv-H ■ • ■■ 

star . like the Sun. 

C, . . \ ' ' ' - . ■ ■ ■' ■ • • ' 

in^^the UMS stars the core contraction phase (2) (3) is relatively 

rapid^ requiring on the averagf^ only 2-5Jf of the time from (l) (2). In 

these heavy stars, the core becomes hot enough at. about (3) to ignite H- 

burning in the surroimdlng shell. The exhausted core may or may not become 

isothermal for a time after the H-burnifiig shell has been Ignited. It all 

depends upon whether 'jpr not the He core is more massive than 10$ of the - 

star's mass. For stars more massive that '6 an isothermal cofe never 

develops, and the gravitational contribution to energy production continues 

to be effective, raising the core central temperature to the triple alpha 

ignition point, straightaway without a pause for the core to gain more He 

mass from the H-burn±ng shell. 

/ ■ '\ ' : ' 

There is one more .f^fiture ofvLMS stars to be mentioned. For staris of - 

< 1,3 the isothermal core .is dense , enough and cool , enough to become degea- 

erate even though its temperature may be 15 million degrees I As described 

earlier, degeneracy is a dense, almost incompressible metal -like state of 

^igh thermal conductivity. Incompressibility arises f rom tiie Pa^li Exclusion 

Principle applied to the densely packed electrons , preventing them from , ' 

.crowding closer together. A degenerate core is not bound by the 10% mass . . t 

limit discussed ubove; it can support the entire mass of a 1.3 stsr while 

'Still remaining degenerate and isothermal because of its high ..thermal con-o . " 

ductlvity. ;Now tlie degenerate core of such a star is not absolutely incorapressil 



only relatively 80, as Is a natal compared with a gas.. Such a coVe can 

contract, gravltatlonally and become- hotter aa .a 'consequence when more matter 

(He.) is added to, it from an H-burnlng shell. It r^ins isothermal but 

growB hotter as the shell burning continues until eventually the 3« 

Ignition point is reached. * • 

Once He-burning begins, the core temperature rises very rapidly. 'iTie 

h^sat generated by these reactions does not, have time to escape ^so the temper- 
ature also rises very quickly, causing the 3eC procesd to go even faster. 
The rate of core heating consequisntly accelerates -until the, pressure due to 
thermal motions of the nuclei exceeds that due to the degenerate electrons, ' 
and. the degeneracy is suddenly lifted. When this happerfs the core ""goes, J 
normal" again and becomes a conventiona] hot gaa once more. It is aa if the- 
"metallic" core were Suddenly vaporized. This thermal runaway with removal 
qt degeneracy is termed a "helium flash", and it results in a sudden reorgani- 
sation of the star's interior structure that produces a corresponding rapid 
'increase in luminosity and radius of the star. Such a sudden, extensive read- 
luStment of the "structure of a star ought to be observable , but two factors 
act to make such observations difficult. Most of the heat of the expanding 
core at the helium flash is absorbed in heating the outer layers of the star 
beyond the core where most of the star's mass resides. Thus the helium flash 
is -heavily blanketed by the overlying portions of the star, and the rise in 
surface luminosity is not very large, certainly far less than the transforma- 
tion of the core might suggest. The second point that affects observability 
is that the flash is' short-lived, and at any one short period of time, say ^ 
the past century during which photometric observations have been made, the 
chances of seeing such an event are rare. .While the Hielium flash peak itself 
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may not tver have been observed, the slower phases leading up to and beyond,, 
the flash appear fairly clftTvrly in thi H-R diagrams-: of star clusters . There 
1b one more aeiective effect: only stars in a narrow mass range between 
a 1 - 1.3 ftt>pear capable of iindergoing a helium flash. -Stars more maaaive 
than i» 1,3'M^ don't develop degenerate corer-, those leas maasive than « 1,1^ 
cannot develop high ehc?ugh central temperatures to igaite the 3 « proceaa and 
therefore have no central nuclear energy eoiirce to produce a Ijeli^ fl^sb. 
' ^ Returning to consideration of evolutionary tracks in Figure 5, the path 
between Points 3 ahd 6 represents a phase in which the core has ceased to bum 
H, having nin out of fuel, and the energy of "the star is supplied mainly by a 
hydrogen i)urning shell. This shell eats its way outward in mass but doean ' t 
actually move outward ^radially. ' It is a kind of incineifator into which gravity 
keeps pushing more fuel from outer regions of tlie star. In fact, as time goea 
on, the H-buming shell actually contracts in radius making the inert core 
plus sliell more and more compact while the envelope grows radially. The star 
with an enlarged envelope has a much larger surface area and therefore. doesn't 
require as high a surface temperature to radiate away its energy, so the H-R 
track moves to the right to lower teiqierature as the star grows in size. 

Stars less massive then 9 Me develop deep convective envelopes at Point 3 
and behave very much as pre-main sequence stars do as they contract along the 
"Hayaahi track% a phase; also characterized by. extended convection. The red 
giant star that is deeply convective evolves along a Hayashi track, but. in 

* & ' ■ . ' 

reverse ''from Point 5 to Point 6. At Point 6 helium burning begins "if the star 
is massive enough (t I'H^) and a helium flash may follow %riilch would. add a, 
small nearly, vertical spike at Point 6. 



stars less massive than « .5. cannot become hot enough to ignite the 
triple alpha He-burning process, and their cores remain degenerate until i 
eventually the hydrogen burning shell runs out of fuel and extinguishes, j 
Thereafter the star has only limited reserves of gravitational, energy, and ^ 
when these are soon pxhauated,the star simply cools slowly as a He white • 
dwarf. In the process the star contracts to a very small radius of the orler 
of a few Earth radii. We have seen the trend toward compaction of the corf 
in red giants in the manner in which a hydrogen burning shell is fed from the 
outside by material pulled inward by gravity, while at the same time the 
shell actually contracts in radius.* This characteristic trend in the giant 
phase toward concentration of most of the mass of the star into a dense, com-r 
pact core surrounded by a distended, tenuous atmosphere reaches .extremes in 
white dwarfs. In these stars the' cores are all that remain excerpt for Just 
a trace of atmosphere bo\md closely to the surface by the strong surface 
gravity of so much mass in such a compact ball of matter. Virtually the whole 
star is degenerate, therefore more or less isothermal and rigid - a massive, 
very hot, very small cinder of a star left to radiate away its thermal energy 
until it^ grows cold. ^ 

The cooling time for a white dwarf (W.D.) may be calculated simply. A 
particle (nucleus, ion, electron, atom) in a star at tempieratiire T has thennal 
energy l kT (k is Boltzmann*s constant), and a star of mass Mgjj^^: consisting of 
pure He, fully ionized, has total thermal energy aTM where a is a constat t)|at 
* includes k and the mean mass per particle of stellar material (He '•^ + 2e") . 
The luminosity L is f^imished by radiative loss of the thermal energy so 

L = d_ . aM dT 

dt dt 



The tl>6ory of the structure of white dwarfs shows that L « KT^** where K 
la another constant; Eliminating L between these two equations gives 

. dt aM ;* 

. i 

This has a ■inple solution, t « AMT"*^*^ « BMT/L and A, B are constants. j 
The time, T, calculated in this way is a cooling time, roughly the time it 
takes the star to lose a significant fraction of its initial thermal energy. 
The table below gives cooling times for He white dwarfs. ^ 
Table III. Copling Tlmejs of He White Dwarfs 

T (years) 
9 X lOV 
U.8 X 10*° 
2.U X 10** 
5.5 X 10* 
2.9 X 10*° 
1.5 X 10** 

All of these except the two 'with L/L^ « 10**^ have cooling times that exceed 
the age of the Universe so we can expect the W.D. phase of a star's lifetime 
to be comparable with its M.S. phase or much longer. 

Observed T and L values of white dwarfs place them well below the M.S. 
and somewhat to the left at high temperatures. As they cool, they move along 
H«-R tracks of nearly constant radius due to ihe relative Incompressibility 
of degenerate matter. Such tracks roughly parallel the main sequence emd lie 
below it by factors 10*^ - 10-^ in luminosity. 



M/M^ L/Lo 

•5 lp~* 

" ,5 10-V 

' .5 10-** 
°.25 " 10-^ 

.25 10-' 

.25 lO"" 



Returning again to the He-core burning stage of . the giants massive > / 
enough to have ignited the triple alpha phase, once this new energy source 
comes into operation the star settles down again into a fairly normal stSruc- " 
tore rather like an UMS star. In fact the evolutionary track from Points 
6-8 Is very much like a replay of the initial approach of the pre-MS star , 
to the MS along a Hayashi track followed by its long ^stable residencp on the 
Main Sequence. The inner structure of the star at this phase is of course 
more complex than a MS rtar, for it has not only a turbulent convective 
He-burning core but also a H-burning shell surrounding the core. When He 
begins to fuse in the core, the latter expands and core expansion is accom- 
panied by a contraction of the distended red giant envelope « The resultant 
slow rise in surfau:e temperature moves the star leftward, in Figure 5 to 
Point 8, a process requiring a few million years for stars in the 5-15 M^ 
mass range. At Point 8 He exhaustion occurs in the core. This happens all 
at once throughout the conyective core much as H exhaustion in the UMS stars 
occurs. As the star evolved to the right in Figure 5 in the latter case, so 
it does in the present stage when He exhaustion takes place in the core. The 
path on the H-R diagram is from Point 8 to Point 9 in Figure 5. An inert 
carbon core with appreciable amounts of oxygen is left after He-burning. This 
inert core will generally be surrounded a He-burning shell and, farther 
out, a H-buming shell. If the star is massive enough, it may attain central 
temperatures high enough to ignite carbon fusion to yield a variety of heavier 
products. If t'he core is degenerate when this happens a "ceurbon flash" will 
occur similar in some respects to the He flash. 

As stars between about 5 and 10 soleur masses swing back and forth on 
the H-R diagram, in their varioua giant phases, they cross a region of inatabllJ 



marked by dashed lines in Figure 5 and designated as a region of classical 
cepheidB. These are radially pul3ating stars whose iumlnositles vary period- 
ically with periods that are related to the mean luminosity (the higher the 
luminosity, the longer the period). These stars have plsyed ah' important 
role in observational astronomy as distance indicators. They are \bright and ' 
readily recognizable^ Since their measured periods of brightness variation 
serve, to determine their luminosities , measurement of their apparent magni-* 
tudeSj^r brightness as seen from the Earth, serves to determine their 
dlstcmces . . ' ■ " . " ' / 

The scenario of further evolution of massive stars can be inferred from 



what has been said above. If a star with a, carbon core is not massive enough - 
to ignite carbon burning, it may eventueuLly become a pure carbon white dwarf 
consisting eilmost entirely of carbon nuclei and^ degenerate electrons much like 



the helium W.D. discussed above, or it may have a more spectaculeir fate In 

. ■: ' ^ ■. ^ ' ■ " - ■ .. ■ 7 • • ■ . . •• 

store as discuased below. / 

/ 

If C-buming is ignited, the star then makes another left-right, looping 
excursion on the H^K diagram similar to those following eai^der core igni- 
tions. Less energy is released from these carbon and subsjequent burnings so 
they take less time; evolution speeds up rapidly, and' each' new burning phase 
leads subsequently to ,a new shell burning phase. The structure of the star 
becomes more , complex as layer develops upon layer. But there are two possibl^e 
ends to this sequence. 

As evolution proceeds the central temperature continues to rise, and 
when it reaches approximately a billion degrees, the star begins to lose energy 
in Bignificeuit amounts by the emission of neutrinos. Neutrinos are produced 
in all the earlier nuclear fusion processes; of course, and carry away some 
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ennrgy^ but energy loss by neutrino emission amounts to a f ev percent which 
is not enough to upset the steur's structure or to affect its energy bcdance 
appreciably. Hovever, above a billion degrees severed nev processes become 
possible that are copious producers of neutrinos, and these processed iacz|9a8e 
very rapidly with temperature. The significance of neutrino production is, 
that neutrinos, once, produced, hardly interact with matter at all. In fact 
a neutrino may traverse = 3000 light years of lead before interacting with a 
lead nucleus I Consequently they are able to leave the star without being 
stopped and are thus able to transport energy directly out of the central . 
regions of the star at the speed of light with no hindrance. A normal main 
sequence star suffers an energy leak of some 2 « 6^ by neutrino emission, 
but a massive star in the late stages of its evolution produces neutrinos 
at much higher rates than a main sequence star as a result of processes other 
thmi normal nuclear reactions that take place at very elevated temperatures 
and densities. The enhanced neutrino emission can quickly remove large am its 
of energy from the star triggering a collapse. This neutrino loss has been 
referred to as a sudden ''refrigeration of the core" that drops the pressure 
in the core which sustains etll the overlying stellar materieO.. With its 
foundations suddenly removed, the outer, matter collapses inwaurd in a stellar 
implosion. The refrigerated core rapidly shrinks to a' degenerate ball upon 
which the infalling stellar matter impinges, bouncing back in a cataclysmic 
shock wave probably enhanced by nuclear burning of infalling fuel (H, He, etc.) 
suddenly thrust into the billion degree interior of the star. The' resulting'^ 
explosion gives rise to a supernova whose peak luminosity 10'° L^) may well 
exceed that of Its entire parent galaxy. The explosion divests the core of 
most of the outer region of the star , usually several solar masses pf material , 



and produces an actlvei Irregularly .expanding nebula typified by the Crab 

Nebula. \ 

The extremely compact, degenerate core that remains ia a neutron star, 
probably rotating rapidly and obsenrable as a ptilsar. This neutroh star 
represents a different order of degeneracy than appears in the white dwarf. 
In the latter the electrons were so tightly packed that quantum effects aaso-' 
elated with the Pauli exclusion principle began to restrict the way in which 
they could be packed. In the neutron star few electrons remain, nearly all 
having combined with protons to form neutrons. But now the latter are so 
densely packed that their packing becomes restricted by the same Pauli prin- 
ciple applied to the neutrons rather than the electrons. The net resxat is ^ 
similar, however. The degenerate neutroiil matter behaves much as the deg^enerate 
electrons in a white dwarf in that it is virtually incompressible and hfUB high 
thermal conductivity. It also has some other interesting quantum propertieis 
such M superfluidity and superconductivity according to some theorists. 

Degenerate electrons are incapable of supporting greater gravitational 
pressure than that produced by about l.j^ so a degenerate core more massive 
than this cannot be a stable white dwarf and is destined to collapse to the 
next stage of degeneracy at the neutron star stage. A star whose initial mass 
lies between l*k and apprbxlmateiy 2 solar masses is believed to form a 
WD eventually by shedding all its mass in exceGjS of l.U M^. It does this 
late in its red giant stages mainly via radiation pressure blowing away its 
distended outer atmosphere in a kind of enhanced "solar (stellar) wind". 
Red giahts are believed capable of losing as much as 10-'* per year in this . 
way, and the giant phase lasts long enough for the star to shed more than 



1 to reduce itself below the l.k limit for white . dwarf s . This suggests 
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that the u«\mlL fat«> of stars Inltlaliy l«sa massive than « 2 M9 is to beco0O 
white dvarfs, and obserratlonal counts of WD' s tend to confirm this viav. 

The process of mass shedding mentioned above results In ejection of 
stellar matter with an outward radial veloc;lty of several tens of kilometers 
per second, and the ejected matter forms a spherical nebulosity around the 
central star. It Is very tenuous material apd Is optically excited by ultra- 
violet from the hot central star. As a result the ejected ms-tter appears to 
a distant observer not as a spherical ball of gloving gas but rather as s. . 
thick, faintly glowing, ring, resembling a smoke ring centered on the parent 
star. These planetary nebxilae, as they are called, are beautiful to behold 
and are sufficiently numerous to support the view that moat stars somewhat 
heavier than the sun up to « 2 form them late in their evolutionary life. 

The fates of stars e 2 has been described as explosive, ending ao 
supemovae that leave behind neutron stars and lose much mass in the explosion*. 
Observations of thn. rates of occurrence of supernovas have a large uncertainty, 
but rates at the upper end of the obsejrved rcmge are consistent with the 
view that all single stars I 2 become supernovas. 

Supernova explosions play an important role in the evolution of galaxies, 
for they provide a means, wherel^ the interstellar medium of gas and dust be- 
comes gradually enriched in heavier elements than hydrogen and heliuM tha^ 
we believed to be the initial material from which the galaxies and their first 
stars were formed. There is no doubt that the envelope ejected in a supernova 

explosion will be enriched in the normal products of nuclear fusion: He, C, 

■ ' ■ ' . » ■ ! 

0, N, Ne, F, Mg, Si, Fe and some intermediate elements . t There are also nuclear 
process that take place in the short time during which the inf ailing gas is 
heated to billions of degrees. These processes build up other elements lighter 
than Fe besides those listed and also build up elements heavier than iron all 



the way to uranium and a little beyond. In fact, SN explosions may well 
be the only way in which an element like gold can be formed. It is sobering 
to reflect that one's ring or the fillings in one's teeth were bom in the 
heaarb of one of nature 'BBwet awesome events I Stare subsequently formed 
from this material enriched with heavier elemi9nts evolve somewhat diTfferently 
in certain details than their earlier generation counterparts, and clues as 
to the ages and evolution of galaxies are contained in the relative numbers 
of metal-rich and metal -poor stars. 

In discussing neutron stars a hierarchy of degeneracies was mentioned: 
first electron degeneracy stabilizing the white dwarfs, then at greater den- 
sitiea there follows neutron degeneracy stabilizing neutron stars. Is ther# 
a limit to the mass that neutron degeneracy can support as there was for 
electron degeneracy? -The answer is a definite yes, but the exact limit of 
the mass is not known as precisely for the neutron star as for the white, 
dwarf because the properties of neutron matter are not well understood. The 
limit probably lies somewhere between 2 and 3 Mg, , and. a neutron star more 
nrassive than this appears to have no further brakes to its collapse. It 
shrinks until its gravitational pull will not even allow light to escape 
from its surface, and it beconies a "black hole", effectively cut off from 
our Universe by all save its gravitational pull on other bodies. This must; 
surely be the most extreme fate conceivable for a star. 

So far we have considered only the evolution of single stars unaffected 
by any neighbor. Stellar evolution takea on a new degree of complexity when 
we consider evolution of pairs of stars bound together by gravity into close 
binaries. Such binaries constitute over half of 'all the stars. 

Consider a close binary system of two stars of somewhat different masses 
on the main sequence. The more massive star evolves into the red giant stage 



f Irat , and in the^ proceis its swollen atmoBphere comas imder the grarltational 
influence of the other member of the binary. When this happens, a sizeable 
fraetion^and in some cases all^of the atmosphere of the giant may be pulled 
off by the less massive star. This material forms an accretion disk about 
the latter^ steur and leaves behind t^e hot core of the former giant, now in 
the form of a helium vhite dwarf. 

The star that has accreted the former giant's atmosphere is now more 
massive than before, and its evolution to the red giant phase is considerably 
speeded up. When it expands to the red giant configuration,' its envelope c<>men 
under the gravitational influence of the WD companion and may be captured by 
the latter, reversing the original exchange of matter. When this happens, the 

star that has recaptured its envelope duiy be able to proceed on its evolutionary 

f ... ■ ■ . .. 

course,. and if it Is massive enough to undergo a second stage of red giant 

3 ' • . ■ . « • 

expansion (when He burning ceases in the core), another transfer of envelope 
may take place. 

The extent of envelope transfer and the number of times it can take place 
depends upon the masses of the stars of the binary system and upon their opera- 
tion. "With these additional vauriables the variety of possible evolutionary 
scenarios becomes quite large Indeed ^ But there are still more ;>os8ibilitie8 
to complicate the picture. Suppose one of the stars becomes a supernova which 
can eject a considerable amount of matter* from the system. If the mass ejectied 
is less them half the total mass of the system, the binary stars will remain 
together as a binary, albeit with quite eccentric orbits. If the eject mass 
exceeds the above limit » the binkry system will be disrupted and each steu: 
will go its way as a relatively high velocity single star. ^. 



. Under proper conditions itcretion of oass by one member of the binary 
^ \ \ \. 

can lead to formation \f a black\hole. At the present time the mpat likely 

candidate for an obRerve\ black hdle, ia an x-ray emitting binary, Cygnua X. 



Such a star must originaliy\ have been too maoaive to form a dajgenerate core, 

and in a Bupemovst exploiionX It must h»ve fo^e*> co'^® *°° "^^^f^® *° 

■ ■'■ \ ' 

atop at the neutron star Btage\ Such stars vere probably initially above 

\' ^ ■ \ • ; 

3 M^, though this lower limit ia Vather\uncertain at our preaent at'ate of 

' • ■ ■ ■• .; ' \ . : • ■ - ■ ■ 

knowledge. A 

In jbhe process of accreting its partner's envelope a star pulls the 
acquired gas rapidly onto Its hot surface, the process is often energetic 
enough to produce copious x-ray emisBion. / Tlbus many binary stars shpuld be- 
come strong x-ray emitting sources during accretion phases of their evolution. 
If one member has evolved into a black hole; it\ should be a particularly 



vigorous x-ray emitter since the gas accreted.by; the black hole is accelerated 
to near the velocity of light as it approaches t^e black hole "surface", and 
collisions between such energetic gas atoms produce not only hard x-rays but 

gamma rays as well. " \ 

Accretion under less dramatic circumstances catn also produce rather 

- A -"' 

. spectacular results in the form of nova o^'^^):i\xrBtB: J^oyM are believed to 
occur only in close binary stars, and current opinion; ascribes these impres- 
sive mass ejections of novae to explosive burning 'df hydrogen that is rapidly 
atfcreted onto a very hot core such as the He white-dwKrf left after initial 
envelope transfer. Nova explosions are far less energetic than supernova 
events, and they eject relatively small amounts of; matter, of the order of 
0.1 M3 or less. This mass ejection can -be accp^ynt^d for by accretion theories 
but it should be emphasized that present th^jjoreticaT understanding of details, 
of binary evolution is in an early state, and much work remains before the com 
plications of the life histories of these stars are unravelled satisfactorily. 



STELLAR STRUCTURE AMD ITS DBTERMINATIOH 
I* Introduction 

' ■ . * * .** . 

Mtai's ability to deduce the structure end conditions of teopereture^ end 
pressure Inside a star that can only be observed from a gi'eat distance must 
surely rahk as one of the crowning achieveaeats of twentieth century science. 
Present knowledge of stellar structure and Its change with time rests upon a 
large base of astronomical observation and makes use of knowledge from many! 
branches of physics. " 

■ • • » • . 

The term stellar structure properly Includes both the study-^of the interior 
of a star and Its outer layers or atmosphere, and since these two topics. employ 
different methods of study; there has been a dlvlson of the subject into two 
special fields, interiors and atmoaphereB. Most of what we know about a star is 
learned from the light that it emits , and this light is generated in the outer 
reaches of the star, in its atmosphere. It is therefore of vital importance to 
understand the processes that go on in the stellar atmosphere in order to be 
able to relate the observed parameters such as luminosity, effective temperature, 
and chemical composition to these variables deeper Inside the star. The Stellar^ 
atmosphere is a boundary layer between the interior and the space into which the 

* . ■ . ^ 

star radiates such infomation as we are priveleged to receive, and the pro- 
perties of this boundary region must be understood if we are to be able to inter- 
pret observations in terms of interlc^- structure. 

By good fortune, it turns out that we don't need very detailed- knowledge 
about the properties of the stellar atmosphere to deduce rather good models of 
the internal structures of most stars. 3y most stars we mean those on or near 
the main sequence, the white dwarfs, and 'red giant s./fj^ ea^iy stages of their 
evolution away from the main sequence. The distended atmospheres of older red 
.glafats and special features of the ^atmospheres of some other types of stars' in 



late •t««ei of evolution play iaiportant rolea in determining the deeper atruo- 
turea of theae atara. Detaila of the atmoapherea of auch atara must be taken 
into account if their atructurea are to be reliably deduced. We ahall confiM 
ouraeivea to the atudy of the internal atructurnia of normal, main a«qu«ice atari 
for the moat part, and for theae the effect of the atmoaphere can be aimulitad 
rather aimply. We ahall alao be interested to aee how the interior of a atar , 
changea with time aa more and more of ita hydrogen fuel supply becomea converted 
-by "nuclear fusion into belium and heavier elements. Theae evolutionary effect8\ 
of exhaustion of the. nuclear fuel cem be followed aa the star leaves the main \ 
aequence. Reaulta of extenaive computer-aided calculationa of the changing , 
structure of- the star-give a picture of the various scenarios a star may be 
expected to" follow in its later life depending upon ita initial m£.Bs and com- . ~ 
position. This constitutea the basis of. oBtellar evolution. The main 

topic we will puraue, however, is an exposition of the general method by which 
details of the internal structure of a main aequence stiur can be found. ' The 
phyaical basis for these calculationa turn out to be diaarmingly simple, and 
once the baaic stellar structure equations have been derived, mainly through : 
atraightforward physical reasoning, they can be uaed to calculate aomeslmple4. 
a.tellar modela. - 

A atellar model is here taken to mean the radial dependence of aeveral. ' 
pertinent variablea, usually temperature, presaure, density, composition, and 
the net - energy flow outward through a apherical aurface aV radius- r within the 
star;": Thte mass contained within radius r is also generally determined as a . 
function of r. These variables presented either in tabular or graphical form 
constitute a stellar model. 



Thm object of stellar interior calculatlona is. to devise a laodel vhose 
total aass, M» Xtaainosity, L, smd surface temperature, T, corx'espond vlth a 
real obsr ;ved star. It might be suspected that many different aodels could, 
give the same M» L, and T^. If this were the case, hov could the correct model 
be identified » if indeed there vere a "correct*' modelt Under these circum- 
stances it vouliL^ hard, to do any meaningful i^tellar interior astrophyslcjt*^ ;^ 



The field of modetl calcxilations would lack the essential' reqiiisite of any good 
science » namely \mique contact with observation. Fortunately for the devlslop- 
msnt of our knowledge of stellar interiors and evolution » it turns out that 
with a few rare pathological exceptions, the structure of a Star is tihiquely 
determined by its total mass and chemical composition. Thus, we are assured 
that if we can calculate a model that giv<!id the same observable properties as ' 
.a real star, e.g., M, L, T-, and composition, y, this model will be \mique and 
therefore the correct structure of the corresponding real star. Thitf uniqueness 
theorem is due to Vogt and Russell and is a consequence of the mathematical • 
structure of the equations that determine the equilibrium internal structiu^e 
of a star. A more detailed discussion of the Vogt-Russell theorem and its 
limitations can be found in Cox and Guili.^^^ 

II, Obaex^vational Base 

Let US consider briefly the astronomical qbservations on stars that make 
u]p the data base upon which the theory of stellar structure rests and against 
which the restdts of such theoretical calculations must be< tested: These ob- 
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servations can\be grouped into several classes: 
• .., " , • *a. Mass . " . • ■ ■ 

^ Stellar masses are determined from jieasurements of, the dynamics 
of binary star systems. The measurements Include a varljety of 
of different types of .normal stare,. but the number of precise 
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meaaxirements 1^ not very large, nunAierlng perhapa a ftv liundradt i 
80 thla tiody of ^ observational data ia rather limited. 
Luminosity " 
The inyerae square lav ^applied to the apparent brightness of a 
star observed photometrically on fe^th gives its luminosity, L, 
provided its distance is knovn -^d provided that its total raAliant 
energy is measured by the a«<tronc»ier . Usually, the total radla- : 
xion cannot be measured, hnb only the radiation falling into a 
limited wavelength range. In this- case the total radiation can 
usually be Inferred from the spectreJ. distribution of the radia- ' 
tf on actually measured in the restricted region of the spectnan. 
Except for relatively nearby star^ vfaose distance can be measured 
accurately by triangulation (some 6V00 stars), stellar distances 
must be determined ^ various indirect methods subject to more 
or less error.. In the matter cf deducing the total radiation of 
a star liefer the whole spectrum froia photcnnetric mer^suremeiits made 
in the visible part of the spectrum, ther^j? is considerable room 
for error especially for very hot st^rs that radiates much \atra^- 
violet light and for cool, red stars that exalt much Infr^ed radia- 
tion. ; ■ ^. _ 

Temperature 

The surface temperat\u'e of a star (actually the ^tempOTature of the 
relatively thin shell, the photosphere, where the- visible radiation 
of the star originates) can be determined by observing the cc7or: 
temperature, or by deducing the temperature from relative intensi*- 
ties of several spectral lines, or from a knowledge of the star's 
radius (rarely known from independent mea^rurement , with the notable 



' ttxeeptlon of tte Sun). Th« surfaea tcaperaturt it aoat- oft«A 
found tigr mea«ui>lng th« oolor temperatura, i.e., the t«Bqp«rattnN&: ; 
of a blackbody that would hava th« same oolor (spectral dlatri* 
>mtlon) in the vlaible region of the epeetrum at the star unter 
observation. The temperature in the Interior of a star beloy - 
tho photosphere is inaccessible to observation, of couree* 

■ Radius ' ' • ..: 

Sxoept for the Sun vfaosd dlaaeter can >• nfsaaurad tyy:. trlangi£Latlon 
and a taiall number of nearby star a, nostly large red giants, iiihoae 
diameters have b^en measured interferometrically, and a number 
of eclipsing binaries whose light curves give the diameters of 
the components, ve have no other direct obserrational knowledge 
of stellar radii. Hovever, the luminosity ^ surfacia temperature! t 
and radius are related by the Stefan-Boltanann law, L ■ Uir a R^**; 
Thj>s assumes that the star radiates as a blackbody lAlch is a 
reasonably good assumption for the majority of main sequence stars. 
Departures from blackbody radiation can be taken into account by 
defining an effective blackbody temperature, T^^j, to be used in. 
place of Tg.. In either case the stellar radius is related to 
luminosity and temperature by a Stefan-Bolt zmann type of equation. 
Chemical CoEposltion ' ^ „ 

The Vogt-Russell theorem identifies composition and mass as the 
main factors that determine the structure of <, a star, so the role 
of composition is certainly of great Imiportance. It Islunfortunate 
that only the surface composition can be determined by observa- 
tion. Stellar spectroscopy is used to determine the surface com- 
position and it is applied In much the same way as it Is used in • 



tht laboratory to identify elements aad aiaqple'isoleeulei ft«^*: ll;? 
th« wavelength* of their cheraeterlitle .epectral lines or tofwadi, : 
These are seen in stellar spectra ialnly in absorption (or 
occasionally as enlssion lines in sboe rery hot stars). The 1 
relative strengths of the absorption lines in a stellar speetna > 
give quantitative information about the concentrations of various 
eleoente In the outer liyers of a star, but one oust rely Aipon 
educated guesswork to obtain fron this surface, coaqposition an idea 
of the Internal coi^position. The two- arie relatisd in different 
ways in different types of stars, depending to a lirge degree upon 
the amount of nixing of- the surface with undorlyli^ nateriai tirougl 
the action of turbulence and convection.- In some sVars this alzlng 
is imix>rtant, and in others. it is not. It is therefore a natter ; 
for ejcpert" Judgment to relate observed chemical abundances at the 
surface of a given star with the interior conqposition needed to 
determine unambiguously .the stars internal structure. The slovlest 
assumption is that the stellar material is well mixed in the early 
stages of formation of the star and that the star begins life on 
the main sequenpe with a homogeneous compositioa, uniform through- ^ 
out and the same as that observed in later life in its surface 
layers. This means simply that the observed surface composition 
is assumed to be the same as that of the interior of the star when 
it begins life on the main sequence and starts converting hydrogen 
into helium near its center. 

As time goes on, the internal composition wlll^ change, of courts, , 
because fusion converts lighter elements into heavier ones, but • 
•these proceaoes proceed at definite and predictfcale rates. If 



^ tie know thn initial- compogition, v« can keep track of hov it ^: 
changes with time at every point in the star. We rely upon the 
' absence of nuclear processes in the outer Isyers of the stsr to 

preserre the priatordial composition there where the astronomr 

■ * • . *• • , ' ' . ' ' • ■ ' ■ * • . 

can oeasure it spectroscopieally* This scenario (l^oes not bold , 

universally, hut it appears well justified in the iaajbrity of 

stars and serves as a useful voi^king hypothesl;^. On the basis 

of spectral analysis of stellar light ve expect the mass concen- 

trations of nost stars to consist priaarily of hydrogen {^70%) . 

and helium (a25-30J{) with relatively little U3^) of heavier 

elements, chiefly neon, oxygen, nitrogen, with still nailer 

lunotrnts of carix>n, magnesium, silicon, sulfur, and iron (^0.1% 

each). . Other elements than these are usually present only in 

trace amounts in the outer layers of the star, but deep inside 

where nuclear trimsformations take place, the heavy eloients may 

build up appreciable concentrations in time. 

Besides the observed features listed above, ther^ are others of a generally 

diffez;ent character whose bearing upon the subject of stellar structure is 

either less direct or less well understood at present. Among these are the sur 

face magnetic fields and rotational rates of stars. Most stars have small 

magnetic fields similar to the field of the Sun and amounting to less than ^ ' 

about 100 gauss at the poles. Some stars, however, have strong fields in the 

kilogauss range, while fields at the surface of white dwarfs are thought to 

range from.^lO^ to as high as 10® gauss. ^ 

Stellax* magnetic fields are measured spectroscopically by observing the 

splitting of si>ectral lines through the Zeeman effect. Technically the mea~- 

surements are very difficult except in the stars with v>&ry strong fields 
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b«0ftui« in nott weak-field stars, the spectral lines are not rery sharp* ^^ley 
are broadened by" effects of collisions and thermal relocities , by turbulent 
MDtlons of the photosphere / and by I>oppl« ■hifti originating frw stellar ro- 
tation. The Zeeaan splitting is Igenerally only a minor contribution to o^rerall 
line shape in sueh stars and is obtainable vith considerable uncertainty wlgr 
after all the 6ther broadening effects have been sorted out and quantitatively 
accounted fore It is fortunate for progress in the field of stellar interiors 
that the magnetic field of a star seems to play only a small role in determlnlagj 
the internal structure, of a nbrmad star. ^Present knowledge about the origin of 
the field in a star like the Sun indicates that the field originates from con- 
nective motions of ionised material fairly far put near the photosphere, at 
most involving less than the Outer 20% of the Sim's radius. In this region the 
interaction of the field with the turbulent gas should affect the scale and 
detailed dynamics of the turbulence • Considering how little ve know at present 
about turbulence in general, it is safe to say that magnetic interactions do | 
not seem to play a major role in determining the structxire of the outer convec« 
tive layers of a solar type star. Whether this holds trtxc for stars with au^^ 
stronger fields As open to debate. 

Rotation has been grouped with magnetic fields for a good reason « The twcr 

(2) 

appear to be intimately connected:, and in the case of the Sun, Parker ^ ' has 
shown in detail how the solar field arises out of dynamo effects driven in part 
by differential rotation of the various layers of the Sun. It is also strongly 
suggestive that the strong-field stars are all rapid ix>tators as wei:^ e In 
fact, observation reveals thj^t the younger, more maatsive hot blue stars (those 
"earlier" than spectral type A, having surfai^e temperatures >SOOo\ and masses 
greater than about i.ffsolar laaases) are the rapid rotatojt*a and include all the 
known strong-field stars. Rapid Irotatlon, besides giving rise to strong fields 
^ which should influence convec tive » turbulent regions of the star, also makes the 
star non-spherical.* Rapid rotators are oblatei because of centrifugal forces 

■ 8 7,^ . - 



and can bulge ^ry considerably at the e<iuator . Thl. aeana that the apherlcal 
8y«natry that" 1« so helpful" in reducing complexity of the general edullibrlia 
•tructure equations no longer holds, and the structure calculatioasSjecoBe bucH 
Bore difficult. Bren so. in the cases in iAieh rotational distortion hai been 
taken into account, the structures do not differ greatly froa the non-rotatln« 
stars/ Rapid rotation night be. expected to influence the rate of bms loss 
froa stars via the "stellar vlnd". a counterpart of the solar wind of aatter 
ejected from the Sun. Mass loss is particularly pronounced in the hotter, aore 
aasslTe stars where the radiation pressure is so great due to the high radiation 
density at the surface of the star that significsnt amounts of matter are 
blown away. Rapid rotation should enhance this mass loss from the e«iuatorlal 
regions of the star where centrifugal" force aids radiation pressure. 

Observationelly, stellar rotation la measured by the symnetric broadening 
of spectral lines due to the D6ppler ehlft. Every surface element of the star 
is either approaching or receeding from the observer at a velocity which demands 
upon its position on the stellar aurface, upon the rotational velocity tm. 
star, and upon the angle between the observer and the axis of rotation of the 
star. It is possible, by sep^atlng ^hia rotational part of the apectrid line 
shapes from other apurcea of line broadening, to R£-rive at a atiur's rotational 
velocity. The data actually give n sin e where n la the rotatione.3 v«,locity 

* 

and e is the inclination of the ^ucis of rotation. Statiatical- studies of large 
groups of stars of similar type serve to 8«p<»rate from 6 since the Ifettsr may 
be assumed distributed randomly ov?r all angles. MaSa loss can sometimes be 
mffis^ured directly, the best case being plan».tary nebulae that have ejected mass 
i& the form of a visible shell >«^5io8e extent can be seen and whose density can 
be inferred from spectroscopic data on the shell radiation. .Maes loss from , 
certain other types of stars wuch aa the very luminous , young blue O-atars and 



th* older, post Baln-aequenco K-glanta ean often be deduced froa speetr^d ' 
Masiiremnts of gae ejection relocltlee, egaln Baking nee of the Doppler effeot 

There reoaina to be Mentioned one unique observational aeans for "seeing" 
ixlto the center of a star, neutrino astronosqr. leutrinos at>e ejected In the 
nuclear fusion proceseas that provide « star's enorgjr, and ths r&ts and mmrgj 
of neutrino OBlsslon are calcvlable. Since neutrinos are hardly stopped at all 
hy »atter> they travel outwards froa the center of the star quite uniapeded 
and off . into space. The Sun is. a close enough source » end the only source doss 
enough » to generate a aeasureable neutrino flux on Earth, MeaauresMhts of 
this flux are very difficult because of its snail interaction with matter, but 
Davis and his coworkers have succeeded In establishing the fact that the 
Molmr neutrino flux Is at least a factor k less than that predicted by current 
■odels of the Sun, and therein lies ^ an important current problem in astrophysics 
to be solved. 

Finally, all of the observational information la5)llcitly contained in the 
H-R diagram constitutes a primary pool of data against which to test theoretical^ 
stellar models and their evolution with time. The H-R diagram Is really quit%: 
a restrictive obserrational constraint. Any valid stellar model must fit some-"^ 
where on this diagram in regions where real stars are founds and' it must con- 
form to the other known properties of real stars in that partilof the diagram. 
It must, for example, also fit the obserjred mass-lumlnosity relation for real 
stars, another restriction on possible models. 



III. Equations of Stellar Structure 

As in most other branches of physical science the aatrophysicol determln- 
ation of stellar structiure and evolutio;i proceeds by a regular Interplay between 
theory and observation. Starting with initial, reasonable aosumptions about 
chemical coasposition, a, model. star is calculated, and its mass, luminosity. 
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•urfttec t«ai>ttr»ture and any other obsanrable properties are calculated and 
checked against observed stars to see vhether the nodel star Batches any real 
star. If not, the assuaptlons are changed, a nev aodel Is calculated, and it 
is again, cellared vlth the ganit of real stars. Ifhen the aodel . aatches a rami 
star, the Vogt-Russell theorem assures Its correctness provided conditions for' 
the applicability of the theoron are net. These are quite general, And the 
theorem almost always holds. ^ 

We now proceed to develop the principle equations that govern the interior 
structure of a spherically symmetric star in equilibrium. By equilibrium w 
mean that structural changes may taie place, but only slowly on a time scale 
of at least the order of 10^ years which is the time it takes r star like the 
Sun to readjust itself by the action of gravity to small structural disturbancei 
This time scale is called the Helmholts-Kelvln contraction tlme>, T„. It is 
one of several , important stellar time scales and the one which &>easures the 
relaacation time of a star for departures from energy balance (i.<4., thermal 
equilibrium). It la in order of magnitude, 

y , . ■ ■ . ' . . _ ■ . , 

Thk « 5th ^ (1) 

■ L " • 

where E^jj is the total internal thermal energy of the star. If the star has no 

strong magnetic fields, and if Its internal energy is made up mainly of non- 

relativist ic monatomic particles, EJ^j^ can be evaluated in terms of M and t for 

thie star. The resulting expression is approximately 



Tjjj, «.2 X 107 (^j2 (^) (J.) years 



(2) 



Where M,, L^, are respectively the solar mass, luminosity, and radius. The 
significance of Tjjj^ la that the radius of a star adjuets to non^equlllbrlum 
conditions at the rate ■ 

R(t) - R(o) e"*/'^HK (3) 
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Another tlae scale of Interest to us is the nucleeur time scale » Tj|, that 
■easures the tlM it takes for the properties of a star to change significantly 
as a result of nuclear reactions in the interior. This tine> toOt vill hare 
diissnslons of an energy dirlded hgr ;t;he luainoslty:, in this case 'E||/L ^ere 
Xj^ is* the available nuclear energy. For a star that derives its energy froa 
conversion of to He^y a process vith an efficiency of 0.1% ot the total 
aass-Honergy Mc^t one would expect to be proportional to .007 \Mc^. The pro- 
XK>rtionality factor is about O.l because it has been found that normaJ. iiain- 
sequence stars begin to change Into red giants vfaen about 10% of their Kass 
has been converted into helium. Thus ve find 



Thus Tjj and are of the order of lO^'^ years and 2 x 10^ years respective]^ 

for a star like the Sunt though for very massive stars both time scales become 

much shorter. With the help of the observed mass-luminosity relation « L can be 

eliminated from these time scales. The M-L relation Is usually given in the form 

where n depends some^^t upon mass and ranges from about 2,73 for low 

' mass main-sequence stars to about h for steLrs more massive than about 0.6 M . 

e . 

Observational data on stellar radii also l^ad to the empirical relation R « 
where a also varies but Is approximately 1 for main-sequence (MS) stars more 
massive than «0.6 M^. With these M-L, M-R relations the time scales can be ex- 
pressed in terms of M alone: 

Tjjjj. « 2 X lo"^ I^^Q^ years (5) 

Tjj a 10^0 years » 

^ \M/ 



or 



T ■ 

*HK . „5 . 

r — ■ 10 for all MS ttars nor* BuisslTe than «0.6 M ." 



A noraal star can therefore adjust its atrueture to achl«T« enargy balahi)* 
and hjdroBtatlc equlllbrlwi on a tlaM seal* that la rary short eootparad »lth 
thf seala of aajor awlutlonaiy changes, Tj^. This conclusion establlshas 
tha Talldlty of treating normal atars' as "equlllbrlwa" objects hound togethar 
by gravity on one hand and prevented ftron collapsing on the other hand by tha 
outvax^ flov of energy fro« nuclear sources in the Interior and by the isensa 
pressures generated In the star by the therail energy It contains. . 

This quasl-equlllbrlua property of the star greatly simplifies calcula- 
tlon of Its structure, for It allows us to apply relatively simple static aqua^ 
tlons to the Interior. Consider first the requirement that the star ba In 
hydrostatic equilibrium irtdch means that a given blob of gaseous material In 
the star Is at rest, balanced by the Inward pull of gravity on one hand and by 
the outward buoyant forces of the pressure gradient on the other. This buoyant 
force Is exactly analogeos to the familiar buoyant force in any fluid. Archlaedei 
Principal says that this force is equal to the weight of the displaced fluid 
which ia simply the mass of the fluid times the gravitation acceleration at Its 
rest position. 

Suppose we take as our aample of gas in equilibrium a thin spherical shell 
in the star. Let it have radiua r and thlckneas dr as shown in Figure 1: 
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Figure 1» Thin spherical shell inside a star for vfaich conditions of aqulli- 
hriuB are applied. ^ 

Let the material vlthln this shell have density » p(r)» tenqperature » T(r), and 

pressure^ P(r). The ma8s» dM(r)» of material in the shell is given siaply by 

dM(r) - Uir r^ p(r) dr (6) 
and the vcight (i.e., gravitational force acting on the shell) Is Just g(r)dM 
vhere g(r) Is the gravitational acceleration acting at radius r^ 

In a spherically symmetric distribution of matter the gravitationail pull 
at any point r from the center is as if all the matter inside a sphere of 
radius r vere concentrated at the center • If G is the gravitational constant 
(6.67 X 10"^ dyne cm^ gn""^), it follows that the gravitational force on (dM) 
at r .is» making use of (6),. 

F » GM(r) dM(r) « Utt GM(r) p(r) dr (7) 

. r*^ 

vhere M(r) is all of the mass lying vlthin radixis r; that is» 

M(r) - /dI4(r) » p{r) dr « (8) 

o o ^ \^ 

The btioyant force on dM(r) is supplied by the difference In pressure between 
the inside surface of the shell vhere the pressure Is P(r) and the outside 



surface of the shell vhere it is P(r * dr) - P(r) . Pressure, being 

foi-ce per unit area, exerts a total outward force of W P(r) on the InBlde 
surface of the shell and a net force, due to the pressure gradient. dP(r). of 
Fp - r2 dP(r) 

It i8 this net outvard force on the shell. that is balanced by the Inward pull 
of gravity given by (7). At equilibrium, therefore (7) and (9) must be equal. 



so 



dP(r) » -G f (r) dr 



(10) 



The Hdnus sign has been added to signify that the pressure decreases as r in- 
ereaaea. It must eventually go to zero at the "surface" of the star. 

This, equation (10) together with (6) or (Q) defining M(r) make up two of 
the four principal equations that determine a normal stellar st^i^cture. 

The third structure equation deals with the energy generaiion in the star 
and its connection with the total outward energy flow. L(r). This quantity 
will be called the luminosity at r by analogy with the total energy radiated 
by the star which we have termed ita total luminosity » L. 

Let e(r) be the rate at which each gram of matter In tli* »hell of radius 
r in Figure 1 is generating energy. Then, since there are W p dr grama 
of matter In the shell-.-the-net -contribution of the-shellTto- the-l^^^ 

dL(r) ■ r^p(r) e{r) dr . ^^^j 

This is the third principal structure equation and the first so far to introduce 
explicitly a property^of the material of the star." e. 

The remaining structure equation deals with the processes by which energy 
flows outward -in the atar and. is therefore expected to contain additional pro- . 
parties of the atellar matter. ' The three principal modes of energy transport 
are radiation, convection, and conduction. We rule out the latter as 

■ . ■ 15 ' , ■ 



playing only a-^exy olnor role in most start (oxoeptiont will ba diteuased 
later) bjr noting that even at the Iflnense pressures and densities present in. - 
the central regions of ordinary stars («S.x 10^^ atmospheres and «l6o ffas/cm*: 
in the Sun) the teaperatiares are sufficiently high («15 Billion degrees in the 
Sun) to n^intain nuitter in very nearly an ideal sm state. Likmise conditions 
outvard frba-the center also confirm the validity of the ideal gaseous state = 
throughout the star, again with some exceptions to be discussed later. Oases j 
are Tsmim to he poor, conductors "^of heat so ve expect that energy vill be trans- 
ported through a gaseous star either by radiative transfer (stepvise/ absorption 
and reemiaslon of radiation) or by convection (turbulent vertical motions of 
the gas vith hotter gas from belov rising to mix^and transfer its heat to the 
cooler overlying gas, the cooler gas sinking in the process to become further 
heated belov) . Usxxitlly one or the other of these proeeslses will predominate 
in a given region of a star, though both may be going on vlthin the same star./ • 
The predominance of one form of energy transport over the other in a given stelli 
shell allows each process to be created separately and grea^fcly simplifies mattbra 

Radiative transfer wiil be considered first. Asomeutioned^srli^r^this 
mode of energy transfer involyss absorption of radiant^* energy by stellar matter* 
followed by reemiasion of the energy as radiation again a short time later. 

The process is not very efficient if only a single step is considered because 

■ ^. ' ' . ■ ■ 

very nearly as much energy is reradlated inward as outward from any small region, 
but the absorption and reemission is rapid, so many steps can take place in 
unit time, ftvnd the slight predcniilnance'of outvard over inward reemission, a 
consequ^^ri/ e of the spherical geometry and temperature gradient in the star , 
suffices to make the overoir process an effective means for getting energy to 
the stellar surface, \ * 

The radiation need not be visible light; indeed, it la not, as soon as 
one gees a short distance into the star • There tempei^aturea rise rapidly imtil 
most of the radiation is u.Vr. iet light, x-rays, and eventually g^ona rays 



Mir the center. The priaary radiant energy giv«n off In the nuclear fusion 
. reactioni" is in the fona of ganma radiation, anf it is only after a large 
number of absorptions and reemlssions that the 'energy of this radiation is 
degraded into the riknge we fsniliarly refer to as "light". It is s^aetiMs 
remarked that the center of a star would appear dark to a human observer. 
Thir 5c, ftot tirue because other ridiative processes take place involYlng charge 

collisions that produce visible light (bremsstrahlung and Compton 
ttcattaring. for exam^>le) but it is true that the central re*ioa of a star woul< 
appear far less bright than its teiiperature would suggest. 

The luminosity^ I.(r) . gives the total energy flowing outward through a 
spherical surface of radius r inside the star. The outward energy flux at r 
is therefore Just L(r)/J^7r r^. The temperature gradient dT/dr across a shell 
of thickness dr must be sufficient to maintain t.his radiant flux, We can get 
the required dT/dr by Asimple argument . 

Consider the radiant energy flux L(r)/J»ir as being propelled outward by 
radiation press^ P^ad. In-auch the scune way as hydrostatic pre88ure_-forces 
water through a pipe. If the radiation pressure is P^^^ at r and (P-ap),.,^^ 
at r + dr. then it- is the difference. -dP^^' *ha* gives rise to.,the outward 
flov of radiation, and this differential pressure must be exactly balanced by 
the resistance met by the radiation along its outward p^^h. .The resistance 
is provided .by absorption and reradiation processes iii the shell dr. • • 

When an atom absorbs radiation the process conserves momentum, so the 
JBomentum of the absorbed photon is gi^en W the absorbing atom. If. befpre 
it "later reradlates a photon, it collides with a nearby atbta. some of the mo- 
mentum it picked up from the ab8orption_of_fche .^>hotGn-ls-transfered-to i;he-~— 
second atom. This net transfer of momentum from photons to gaa atoms con- 
stitutes a force acting upon the gas'arising out of the radiation prese.ure. 
A photon . Of energy E carries momentum^ E/c which it transfers to'k gas atom \ 



upon b«lng absorbed. We •xptct, theraforo, th»t ono factor contributing to 
the i j«i«tanc« of matter 4e the flow of radiation irould t>e proportional to 
L(r)/Uw r^c (of the fpr« K/c) since this is the net boacntua carried by the 
radiation flux, L(r)/i»w r^. A fraction of this aoMntus flov Is stopped by 
the atoos in dr by absorption. 

Let the opacity of the gas to radiation be designated by k (oWk letter 
kappa, the conventional notation in the stellar interior literature). This, 
is Juet the usual oaas absorptipn coefficient in physics. If a been of light 
passes throuirh a thickness x of natter of density p, the bean intensity de- 
creases according to the faniliar equation 

^(X) " ho)^'^^^ ' (12) 

which is equivalent to saying that dl/I - -Kpdx. Applying this definition of 
•c to our spherical shell dr. a fraction fcp of aKaentua flux L(r)/Uir r^c is trans- 
fered to the gas by absorption in dr. and it follows that for this to be balance* 
by the differential radiation pressure , we mist have 

- <P^r). dr J 
hit r^c 

Electromagnetic theory provides a sisple relation potween radiation pres- 

> ■ - ■ ■ . 

sure, Py^jj, and the energy density of the associated radiation field/ E 
this is , . 

^rad • ^ Erad * J aT** . (lU) 

where a - Uo/c for radiation In local thenaal equilibrlua at temperature T 
(a is the Stefan-Boltamann constant ), The differential pressure dP is 
therefore simply ii aT^ d'l'^and ^(13) may be written 

dT ^ 3pKL(r) \ 

dr " i6if r.2acT3' (15) 




This la the final interior eqiuition needed if the iter'e- ener^ transport is 
entirely- liy radiation with convection uaiiqportaat Hote that both p and ic 
rmry with r sinc« they ar« functioM of tht Yariables T(r) /P(r) t and oospoti- 
tion. The opacity is another constitutlr* paraaeter specifically referring 
to stellar aatter that has to be Introduced*. 

If the thermal gradient, dT/dr, is large enough, convection vill develop, 
because the underlying gas vlll be hot enough and therefore have low enough 
dcneity cosspnrod vith the gaa in the upper part of our shelly dr, for^the V 
buoyance of the lover afeterial to overcome gravity • . JPiiia lighter underlying 
gam will riae through the^ overlying gas forming' Jbiurbulent convective "cella" 
that tranaport thermal energy,. Frcm another point of viev/if (dT/drV radiative 

ia.too nail (because pK is too small) to allow the energy L(r) to escape by 

** * 

radlatj.ve transfer, the star will find an alternative way for the energy to 
get out » and convection is the alternative mode of energy tra^isfer. 

It is possible to devise a simple criterion to test whether convection can 
'take place at 'any given place in the star. Thiols done by- seeing whether a 
blob of gas of density and pressure p(r) and P(r)""SM the 1^^^^ 
dr would rise or sink if moved outwards by the shell thickness dr and allowed " 
to, expand adiabatically to the dehsity and pressure of the surrounding gas at 
r ♦ dr.,. If ItvsinJcs under these conditions, convection will not take place 
and energy vill, be trtmof erred by radiation. The criterion that results from 
these consi^eratlol^is Is . 

I i (dl) •< I (l - I) idP I for stability . (16) 

^ re^d V Y' P dr 

If the inequality is not obeyed, convective transport" will predbmina^^ In 

this expression y is tiie usual ratio of specific heats at constant pressure an<l 

constant volume for the stellar gas (y - c- /c ). .This ratio enters into ex* 

. • ■ P V 

presslon (l6) because during adlaba1;la^expansion of the test sample of gas 



. aoved fro« r to r ■♦• dr Do'haat is Interehuiged with the •urroundlngg , «n<l 
th« •diab*tic expaniion lav, PVY ■ Pp-'' . eonit ant, mwt hold. - 

Onca conrectiYe traniport has been oatabliihed, a different expresnlon 
than (15) for dT/dr Buat .be applied because the energy 'transport rate is 
different. Schwarxschlld^'*^ glres a simile approzinate theory for conyectlve 
energy floir in terms of the excess temperature that an adlsbatically e-vpcmded^ 
ToluM of gas BOTod to r ♦ dr has over the temperature of the . surrouhdiiii- gas - 
at, r > dr. ■ . . . . ' 

If the blob of gas trayele on the average a distance 1^^^ the "nixing 
length", before becoming mixed vith the surrounding gas and losing both its 
identity and the energy it carries, the final expression for the convoctivj 
heat flux can. be written^**) 

W - Cp' (Mrl,'^ (AVT)^/2 ^2 „ ^^^^ 

where ATT signifies the excess temperature of the rising gas over its sui^round- 
ings above its starting level. This flux w be compared yi^h the radiative 
flux, . ■ ■ . i_ ; •- 

Uir r^ 3 Kp . dr ^ ' 

obtained by revriting (15). The total flux H - H^^ + H^J^^ can be solved for 
the actual temperature gradient dT/dr ujing (17)' and (l8) . When this is done, 
it is fount that ATT, tfie excess over the adiabatic gradient is very small,' 
of the order of one millionth of the 'adiabatic iradient itself, or of the order 
of only about lO K/cm.. A typical convective velocity v « 30 cm/ae'c which is 
several orders. of magnitude smaller than thermal velocities, insuring that con- 
'vective flow does not upset overall hydrostatic equilibrium. The convective 
heat, flov still contains a parameter, the mixing length, which is difficult to 
evaluate. Schwarrschild and others have taken I to be of order 0.1 R, but 



this, it 4 problem w« need not purnue here. Before leaYing thle brief dlecu*- 
Bion of con7«ctlon, it should be reaarked that the convection 1b turbulent and 
rerres to mix thoroughly the contents of any convectlve region in a itor so 
the region ii praticolly unifora in coapoiitlon. 

Since AVT ii so onall cdipared vith the adlabatic gradient. 



dr ad P dr 



(19) 



the latter con be used in place of (dT/dr)^^^ fron (15) wherever test (l6) show 
convection to be present. .■ 

The four (or five^ if we include (19)) structure equations, hav • written 
in terns of P, pi and T, but these ore related by the equation of stava on 
ideal goo. 



which is the way astrophysicists cuatoaiarily write the nore fomilior fom 
P7 ■ nRT. In (20) k is the usuol Boltnnann constont, 1.38 x 10" ergs/^l^'j V 
is the neon moleculor (or atomic) weight of the gas, ond H is the moss of a 
proton, 1.6T x 10~2'*,gm. Solving (20) for p and substltutiog into equations ('8), 
(10), (11), arid (15) serves to eliminate p giving 



(i) ,QHU M(r) p 

^ k r2 T 



(ii) . dM(r) ' U-n H 5 P 
dr. "T-^'^'t 



(iii) dL(r) . _ P 

dr 



(iv) 



(v) 



k .^^T 



(21) 
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7.2 mh 



■^'conv ^^Pdr 



Th««« oquatlona. tog«th«r with » ?K»ov3.«d«« of »c wid c and their depcudtnoe 
on P. T. and with sultabla boundary conditions anabla u« to conatruct atollar 

. ■ " , 

■odols. 

. The boundary condition, will be dlaouaied flr«t. They can generally be 
caet In a very aimple fora. While It is true that the stellar surface Is not 
sharply defined, the density does fall to sero rapidly beyond the photosF^ere 
.In a «»H fraction of tbt star's radius R. Blallarly the tesperature and 

pressure fall to aero rapidly beyond R. so for cur purposes It Is pantlsslUe, 
to take as boundary-conditions: . . . ^ 
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at r ■ R 



(22) 



at r ■ 0 



where >:,h<s .subacrlpts c refer to values at the center of the star. The condi- 
tions M():) - L(r) -.0 follow from the definitions "of M(r) and L(r). Since.M(r) ' 
is the mass vilchln a spherical 'surface of radius r. it clearly -/anishea at 
r - 0. Likewise at r - 0 there^can be nothing left to generate energy so . 
1.(0) « 0. 

The. Bean atonic weight y is often more conveniently expreaaed in tenaa of 
the compoaition.- The conrection In stellar interior work la to expreas compoal- 
tion.ln teraa of three variables X, Y. Z whicli are respectively the maas • 
fractlona of H. He, and all heavier atoaa in a aaaple of atellar matter. It 
follows that X + Y + Z - 1 so only two of theeeVy^riaU^^^^ are needed to specify 
tha coDpofiitlon. v , 



Tho reason vtay ftpeclfyios only H, Ho, and "heevlor^' atoan ii adaqiaatii < 
in ttoat stellar interior vork is tvofold: (a) in normal, main sequenca stars 
the nuclear energ/ generatins p:!2^efises represented by e inrolre conversion 
of H into He and oco»i»ionaXly conrersion of He into C^^ (i,^,, x Y and 
T Z)« The detailti^ oaturi of Z is needed only if one wishes to study the 
precesses of nucleosyBtii^tdis by ^Adch all the heavier elements are built up in 
stellar nuclear reactions; and (b) as a*^ consequence of the high temperatures 
in sili^ellar interiors the atcns are fully ionised, for the most part, and the 
heavier atoms present in appreciable amounts (generally less than 3% by veight) 
ma;,; be coaaidered completely ionised into bare nuclei plus free electrons equal 
in nuMtber to the atomic number of the heavy nucleus, namely « A/2 electrons 
for elesents C^'^ fIS^^f. up vhich are the ones of interest. 

Th® moleculwr weight of a gas is the weight of « 6 x 10?^ particles* 

If 6 X 10^^^ ^tcwBs of (1 gram) are fully ionised, we have 12 x 10^3 particles 

and 1 gram atcgn of these (6 x 10^^ pax^icles) weighs 0*5 gram since on average 

haJLf will be protons and half electrons. The mean molecular weight of ionised 

^^TOffFa is the:refore Similarly a gram-atom of He^ {h grama) when fully 

ionized gives 3 x 6 x lO^^ particles (each He** gives 1 alpha particle plus 2 

electrons ); thus the mean molecular weight of fully ionized He^ is ^. Finally, 

heavy litoms of normal atomic weight A give rise to A/2 electrons per nucleus, 

and the mean molecular weight of the fully . ionized matter is approximately 

2A ' ,^ 2 gram of a mixture consisting of X gma of hydrogen, Y grams of 

A+1 

X ■ 

helium, cuttd Z grams of heavy "atoms", fully ionized, produces 2 — particles 

from plus 3Y particles from Hel* plus ^AZ^^ z particles from heavier atoms 

. ^{j ~ -V AH"'v [ ---- •■ 

vfaere H is the mass in grams of a single hydrogen atom. The mean molecular 
weight of thiS" graun of matter is by definition the weight of 6 x 10^3 particles. 



8ine« I - 6 X 1023. it follow that M !• 1/H dlTld«d by th« su* of partiulo 
jSlv«n 4boT«, or 

u • 1 ■ 
^ 2 

Tn tens of X, 

5X..-Z ♦ 3- "x; ■ 

The range of U 1« from h for pure Iqrdrogen up to 2 for heavy Batter (Z ^ !)• 
An He'* in the deep interior. X decreatei and y increases, raising each 

gradient dP/dr, dM(r)/dr, dT/dr. and dL(r)/dr in the basic structure eauations* 
It is not surprising that the increase in y that results from conversion of 
H-"- into He'* has large effects upon the structure of the star. 

Consider the next constitutive variable, the opacity, K. It is clear fro» 
the structure equations, (15) or (21,iv), for dT/dr that K controls the t«(per- 
ature gradient where radiative transfex^ predominates. In regions of large- 
opacity, the energy cannot escape easily. As a resvilt, the local pressure and 
temperature rise, causing the star to eacpand until the energy can get through 
the high opacity region at Just the rate at vhich it is produced in the interior 
It ir clear « therefore, that the structure of the star vill depend in a sensi- 
tive way upon the opacity, for if k changes, the star has to readjust itself 
throughout in 'such a way as to allow the energy generated in the deep interior 
to get out without being blocked' anywhere. The opacity is« expected to depend 
upon P, T, and composition and will vary widely between the center and surface* 
of a star. 

Accurate values of ic as functions of the variables too cooiplex to be 
represented a simple analytical formula. For detailed structure work, 
values of ic are tabulated as functions of P»'T, X, Z, and used In that form in 



ciMqmttr calculated stellar KOdels. Hovever, analytical representatloiia oac 
be obtained vhlch are sufficiently accurate for approximate structure deter-* 
■inations and vhich hare a further advantase in shoving clearly the functional 
dependence of k upon and cdsposition* 

Except near the surface, tenperatures are so high in nomal stars that 
the atoois are noptly fully ionised, and the stellar aaterial consists nainly 
of bare nuclei, chiefly protons and helium nuclei, plus enough electrons to 
give electrical neutrality. These electrons My be free (not bound to any 
nucleus) or bound temporarily toe nucleus until a sufficiently energetic photon 
comes along to free it by photoionization. 8\ach bound-free transitions consti- 
tute one of the DOBt important processes contributing to opacity. In fact, 
bound- fr^e transitions ordinarily contribute the dominant term in k. 

Hext in Importance are free-free transitions. An electron in an unbound 
state in the ricinity of a positively charged nucleus travels in a hyperbolic 
orbit. If it encounters a photon, it can be. shifted <>to a different \mboimd 
(hyperbolic) orbit vlth absorption of energy from the photon « Such a proeess 
can only take place near another interacting particle (the positively charged 
-nucleus- in- 1 his -case) "!!! -^order " f or^moteenti^^ free 
electrons can scatter electromagnetic radiation by the Thomson scatterixig pro- 
cess. In the cooler stellar envelope, the more familiar; bound-bound transitions 
can take place, giving rise to more or less sharp absorption lines vhose wave- 
lengths are characteristic of the absorbing atom. But this process is of 
negligible fmportance in the interior. . 

The bound-free and free-free transitiqns are much more important than 
electron scattering in moat stars, and it^ is possible to express the b-f and 
f-f contributions to < analytically by a formula first obtained by Kramers. 
We will not derive the Kramers opacity formula but will simply present it 
and mention the physical origin of the various factors that appear in it. 



(25).: 

ic^ » U X 1022 g fig 5 (i-z) (1 ♦ X) ca"^ 

The nvMrioal factors aad* up of phyiloal oonttaati -swiJi at «, h» e, 
•te. Th« opacity vould be •scpeetttd to b« proportional to tha density of the 
gas,, p, since absorption Is certainly proportional to the nuaber of absorbers 
,xu the light path. The temperature factor enters In a different vay. One aay 
say that at higher tenperature the radiation pressure Is higher so the radia- 
tion flovt mor« Mslly through the matter leading to reduotd opaolty. Or one 
aay argue that as T increases |b*-f absorption nust decrease because the chaaoe 
is snaller for a given atom to have an electron temporarily in a bound state 
at any instant. The f-f absorption would also be expected to be smaller at 
higher temperature vhere the electron kinetic energies are higher, the hyper- 
bolic paths are straighter , and the electrons are farther from the nuclei on 
the average, giving a reduced interaction necessary to conserve momentum in 
the process. All of these vays of looking at the effect of T lead to the ex- 
pectation of nw^rnve^^ >c and T. T^The calculations of Kramers 
show the dependence to vary approximately as T'^'^. The Gaunt factor, g, 
arises from the wave nature of the electron and is nearly unity. The ^^guillotine 
factor", t, vas Introduced by Eddington to take into account the fact that once 
an atom has undergone a bound-free transition, it is for a time unable to repeat - 
this process until it later acquires a replacement bound electron again. This 
effect acts to reduce the absorption below that which would exist if every 
atom were available for b-f absorption at all times. The guillotine factor 
depends upon the amount and relative abimdances^^of the heavy elements, upon 
the density, end upon the temperature in a complicated manner, so the valuiss 



ot t ulually t*l)ul*ted. Sl^pltr approximate M«lytlc»l for«ulM for t 
•re •oMtlBtts used In less detailed eeloulatlons . Por our purposes It is 
•ttfflclent to take s and t both of order unity. 

At low densities and high teaperatures scattering of radiation hy free 
electrons (Thompson scattering) can doainate ov»r b-f and >-f transitions 
vhich require close proxlalty between electrons sind .lons. Tboaqpson scattering 
glTes a slMple contribution to the opacity. 

•c, « 0.19 (1 ♦ X) (26) 
This is proportional to the nuaber of free electrons. % (1 ♦ X). as would be 
expected. The number of electrons is slaply the sub of contributions , 

+ IJBL ♦ 1 A Z£ - 1 (1 ♦ X) £. of electrons from X gma of hydrogen. Y gM of 

He, and Z gM Of heavy ^leaents . Remember that H in the denominator ^of the 
suMMtlon is the mass of a proton or the reciprocal of Aragadro's nuitber . \ 

Figure 2 below shows the ranges of T and p in which the various kinds of 
opacity dominate. 




Figure 2. Temperature density diagram for opacity. The run of p, T for the 
Sun is shown. 



This figure Includes a form of conduction that can b« laportant in regions 
of Tsry high densities vfaere the electrons are packed so closely that they 
becQM degenerate as they are in ordinary Mtals at Ipv tenperatxires* Under 
these conditions heat can be transferred ty conduction via the degenerate 
electrons in auch the saae vay as heat is conducted in^ a metal. This prooees 
is Tsry efficient and leads to neaily isothemal conditions throughout the^ 
.degenerate region. Electron degeneracy is cosnaonly found in vfalte dvarf atari 
In the cores of low mass stars late in life, and in masslTe red giant cores 
after eonrersion of to He** has ceased in the core. • 

In ordinary naln-sequence ■tars there are tvo main processss by vlileh 
is fused into He**. In stars Ie«s massive than about 1.2 central teaperatui 
and densities favor the proton-proton chain ^ose principal reaction steps 
are given belov with the energy released in each reaction: 

♦ h1 d2 ♦ ♦ V l.UU Mev (l.U x lO^O yr) 

♦ — *" He3 -1. V 5.U? Mev (6 sec) (27)' 
y..^ He3 ^ He3 — >- He'* * 2a^ * v 12.85 Mev (10^ yr) L 

6* is a jpositron and v is a neutrino. The times in parenthesis are the mean 
times hctvoen reactions for conditions in the center of the Sun. For eacample, 
a proton will live. on the average for l.U x 10^° yrs before it undergoes a< 
collision with another to f ona a deuteron in the first reaction. The D?. has a 
mean life of only 6 sec before it is converted to He^, eaA so on. There are 
other branch reactions in the set, but the three shown above are the main ones. 
These generate a total of 26.7 Mev of energy per He** romed (2 each of the firs 
two reactions must take place to provide the two He ^ atostts for the final react! 
Each neutrino carries away .26 Mev, leGving 26.2 Mev per He** to contribute to 
the star* li luminosity. 
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In wov BMBlVtt •tars in which ctntr*! teotperaturci «xe«ed ■ 15 x 10^ ®K 
th« following cycle of reaction* becooes the dotrlnaat energy generating prooeet 
This let of reactions couprise the carbon cycle, soaetlMS Jloiovn as the ClfO 
bi-ejrcltt since C, I, and 0 are all InTolved. 

Cl2 ♦ Hi — h13 ♦ V 1.9l» Mev (1.3 x lo7 yr) 

— *C13 ♦ e* ♦ V 2.22 Mev (T ain) 

Cl3 ♦"Hi ■^Hl'* ♦ V 7.55 Ker (2.7 x lO^^r) ° 

Hi** ♦ Hi— --0I5 ♦ V 7.29 Mev (3.2 x loQ yr) ^^^^ 

~* »15 ♦ s+ ♦ V 2.76 Mey (82 sec) 

h15 ♦ Hi — -Cl2 ♦ Ho'* ♦ V U.97 Mer (i.i x 105 yr) 

* Abput one in 2500 of the last reactions produces an qI^ and giyea rise to 
the following set of reactions irtilch are not rery Important froa an energy pro- 
duction standpoint but do illustrate the first steps toward building heavier 
elenents : 

jfl5 ♦ fll — *>-0l6 + V 

qI^ ♦Jll . rr*- f11. -»- V - -: 



f17 — ^ 0I7 4 o* ♦ V 

0I7 + Hi — V h1^ <' He"* 
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The carbon cycle uses C^^ as a catalyst since it eventually regenerates the 
Cl2 except on the rare oecaBiona when the second set of reactions takes place. 
The main reaction is again fusion of four protons to form a He^ nucleus. Energy 
released in this reaction is 25.0 Mev, a little less than in the p-p chain be- 
cause more energy is lost through neutrino eaission in the carbon cycle whose • 
neutrinos are acre energetic. The aean reaction times in (28) are again for 
solar cejitral conditions and vary widely with temperature and density. But they 
give a rough indication of the slov, overall rate determining steps. 
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A third kind of nuclear energy generation procean vlll be Mntloned briefly 

li 12 
the triple-idpha process by vhich 3 He^ nuclei fuse to foni C ^ in a.tvo-step 

process vith Be^ as the intemediate product. This process only goes at rery 

high teaqperatures , above 10^ ^» because the strong coulonb repulsion force be«> 

tveen tvo He nuclei requires a large thermal energy to orercoae it during a 

collision in or^er that the nuclei can approach close enoiigh for the strong, 

short range, attractive nuclear forces to operate. 

The rate of energy generation, by any of the above processes flniat con-- 

, * if . 

tain several factors vhich can be Justified by simple physical arginients. 

v.. 

Pir«t, ve note that collisions must occiir between tvo nuclei, such as tvo prp«- v 
tons, if a nuclear reaction is to occur. In fact, these collisions sust bring 
the nuclei vithin the range of the nuclear forces. lOT^^cn) for a reaction to 
occur. It follcivr^ that one of the factors in e must be the number of times 
per second, on the average « that one kind of nucleus collides vith another. 

Let the distance betveen nuclei for interaction 'via nuclear forces be 
denoted by d. This may be thought of . as the radius of the "sphere of influence*^ 
of "^a given" nucleus; if it comes vithin a distance d of ainother vith vhich it is 
capable of reacting, a re::.ctl6n may place. As a nucleus of one kind, sa^ . 
typejl, moves about at velocity v, it sveeps out a "cylinder of influence" of 
voliu&e each second. - 



If thmrm art nuoloi of typ« 2 per unit voIum (vith vhlch typf 1 ou rtMt)t 

• • ■ 2 ■ 

our noTlng type 1 vlll encounter Td V^v type 2 nuclei per eecond, and vlll past 
theie^enoush to cause a reaction. If the particle density of type 1 nuclei Is 
Hi per unit toluMr, then the total reaction rate Bust be proportiqnal to 
wd^j^ffg. This can be conrerted into bums densities fay noting that £c pro* 
portional to pX, and Hg is proportional to pX^ vhere X^ and are t he re spectlin 
abxmdances by weight, llius the reaction rate vlll be proportional to td^p^vXj^Xg 
This ean be further sioplified by teking isto account the vare properties of 
natter. The deBroglle varelength of a nucleon is h/nv, and the collision dls«* 

tance, d, -^i^t be proportional to the deBroglle vavelength. Thus d Is proper- 

- 2 
tional to T"*-'' and the collision factor becomes sitoply'fi. X^ X^. 

T 1 2 

So far only nuclei soTlng vith a single relative Telocity, r, hare been 
Considered. Actually , a distribution of velocities Is present, and at a temper* 
atuiTG T, the Maxvell-Boltuann formiaa for the distribution of relocities in a 
gas ssys that the number of nuclei that hare relative velocities, v, is pro- 
portional to v^T*^/^^ vhere y is the reduced mass of \the colliding nuclei 
(m"^ ■ m^"^ ♦ "2"^^ and k is the Boltzmann constant. This factor must also 
appear in the final formula for e. 

The colliding nuclei both have positive charges and repel each other 
electrostatically. In order to get vlthin range of nuclear forces they must 
penetrate this Coulomb barrier. The Coulomb repulsion^ is very strong at dis- 
tance of the order 10 -'cm, and it would require very high temperatxirea to pro- 
vide enough thermal energy to drive two nuclei this close together. These 
temperatures would have to bo far above those in star centers were it not for- 
the fact that the wave nature of the nuclei permits them to tunnel throi;ifeh a 
Coulomb barrier that they- could not otherwise energetically surmount. Figure. 
k, illustrates this tunnelling proceao Bchematically, 



Coulomb potential barrier 



da Broglla wavt uaoclatad vltli nuclton 



Small probability 
of nuclaon wava 
tunnalllng thru 
barrlar* . . p 




^Nuolaon of tnargy E Insufflclant 

. to go over Coulomb barriar 



Range of nuclear forca 



NUclaar binding anargy 



Figure Schematic rapreaentation of nucleon tunnelling through CoxdiOBb barrier* 



The tunnelling probability, first calculated by Oaaov, reduces the number 
of effectire . collisions by a factor proportional to 



ezp 



vhere Z^^ and Z2 are the charges on the colliding nuclei, and h is Planck* s 
constant. 

Finally, even if tvo nuclei collide and penetrate the Coulomb barrier, a 
nuclear reaction will not cdvays take place. To account for this, another factor 
P, the **nuclear factor" must be included. It is the prbbability that a reaction t> 
vill take place betveen tvo nuclei once they are vlthin the range of nuclear 
forces. It depends veakly upon the relative velocities of the nuclei and r" 
strongly updn the particular nuclear species that &re colliding. 

When all the factoii's are put together, the energy production rate, e, is 
proportional to ; ^ 



colll.lon »ta»«ill-Bolt„am» 0«o» Hucl.ar 
ffcctor factor, factor fwtor 

Thli •xpre.Hlon liiToinni •xpon.ntlal t«nu th«t d«i>«nd upon T oan b« ajrprox- 
l«at«d by Buch norm eonTanlent fonu, 

,.v epp - Apx2t» - (31) 

^afO*^ BP^Cl'" (32) 

vher. A tnd B •P. duaal-conatant. . (that 1.. they chang. .lovly with T). . and 
n are exponents ^that al.o change .lovly vlth T. and X^^ l. the »a.. fraction of 

heavier;i.ateriarthat con.lat. of carbon and nitrogen. X^^^, l^ roughly i Z 
for aolar type atari. ^ 

..>■■■• 

i . . ■ . _ 

In (31) the exponent n 1. -6 at T «!, x lO^ ^k. dropping to -S.? at high 
temperature, around 2 x lo7 ok. The value of ■ in (32) also varies from -20 ' 
at 13 X 10^ °K down to «13 at 50 x 10^ ^K. 

In the range of temperatures found In the energy producing core of the^Sun. 
(l.eV. T «12 to 15 million degrees) n «1» and «»20 while A^IO"^? erg./.ec-g. ' 
«»d B«1.6 X 10-1'*2 ergiT/sec-gm. ' TaWng - i z the energy generation e^ua- 
tlons for a solar type star become, approximately, 

^ epp - 10-29 px2t^ erg^ec-gm (.33 J ^ 

^ciro" 5 X 10-l't3 pzxt20 (3l») " 

The much stronger temperature dependence of e^^ reflects the higher Coulomb ° 
barrier Of the heavy nuclei whose charges are .A/2,. A being their atomic number, 
6 to 8. For completeness ve add the approximate expression for for the 

triple alpha process as given by Schvarzschlld^^) 

■ . ^3€ ' 10-2'*8 p2y3t30 ^^.^ ' 



•of course the relatlYe 



vhloh Is Talld for T «lUO Billion S/^ST^^i, 
abundance of He'^. . ^ 

Figure 5 belov gires the dependences or tm< t^hree c»s upon T for the 
pirticular choices, px^ » lOO. Z - .OIJX, and ^ ^V^ - 10^ for the triple-alpha: 
process: 

♦ 10 



4* 



lo«€ 



-r - 





7.0 



Figure 5. Nuclear energy generatioD rates as functions of temperature. 

The energiec released per graa of tranaauted matter (i.e., per graa of He** br 



formed) for the three processes are ^ . 

Kpp « 6.3 X 10^8 erg/ga Ht'*- formed 

^CKO ■ 6.0 X 10l8 arg/gm H*'^ formed 

830 ■ 6 X erg/gm formed 

With these energy release values and the e's, it is straightforward to calculate 
the rate at which. is converted into He**: 



(36) 



dx . epp 
dt - - ^ 



^CNO 



(37) 



dY 
dt 



These equations, as Schwarsschlld points out, contain the core of stellar 
evolution: the changes in chemical composition produced by. nuclear reactions 



3k 



vlthin ft .t*r. The structure mist chang. In r.0pon.. to chimge. In X, Y 
•ecortlng t6 the Yogt-Ru...ii theorem Wch ..y. that the etructur. 1. dettr- 
ined .olely hy (fixed In • given .tar) and c(»po.ition which change. 



vlth tiM. 



now that e ha. been found in term, of p. T. X. Z. it can be expr»..ed 
vith th. help Of the ga. lav p - ^ ^ ^ 

the .tructure equation (21 iii),. All of the .tructure equation, then depend 
.^pon juat 6 unknown., p. t. L(r). X, Z. If X and Z are known or inferred f«» 
-pectroacopic »ea.ure«ent."of the outer layer, of the .tar the number of un, 

known, i. reduced to 1* and equal, the number of equation. (21). 

Since the four equation, only determine cha^e. in P. T. L(r). M(r) with 

radlua. the procedure for aolvlng the equations involve, aaeumlng value, of P. 
. T. L(r). M(r) .omewhere in the .tar to start with, usually the center or the 

surface, and using the eq^t ions (21) to calculate values at slightly different 

r . These values are then used tl5=>progressively calculate values at another 

nearby r. and so on throughout the star. 

^ — Calculation of-a Simple Stellar Mni!i>i 
- To Illustrate the method by a simple example, assume m Is constant through- 
out the star. Such a.star would have to-be very young, at an age before any 
appreciable fraction of Its hydrogen has been converted Into helium. If ve 
take X - 1. Y - Z . 0. ve have the Bimpiest possible case. 

K.rt. Choc, .ulfbl. ,iau« of P„ T. rt th. cMor vh.r. Ur) ...d M(r) ' 
■ f.... valu.. « will lnt.8r.t. th. .tructur." .quatlon. outvard. di- 

viding th. r.dlu. of th. rt„ mo . „™b.r of thin .h.U. of thlcta... ir. 
A «o.r dlrl.lon «,uld glv. 8r..t.r .ocur.cy. Th.r. 1. .o r.«on vh, m th. 
.h.ll. Should h. Of th. thlcta..., hut dlvHoo into .,u»l 4r »m h. u.,d 
in thl. 0... for lllustratlv. purpo.... W. „hltr«-lly pick .tartlog valu.. of 



■ 8 X 10^ °K and « 8 X 10^^ dyn«s/ea^ b«oause th«s« li« In th* g«a«r«l 
rang* of raluas for a star soMvhat less than on* lolar juuia. Tha gaa lav. 

Pc - - .5 X i.6T X lO-^** X 6 X ipl^ 60.5 g«/c«3. 

1. 38 X 10-16 X 8 X 10^ ' ^> 

1 AB}r of tha four rarlablaa P, T, L(r),M(r) at r can b« oalculatad at 

r Ar by Mans of a Taylor's aerias expansion about r. Vor azaaple, 

P(.r ♦ Ar) - P(r) ♦ g I^Ar ♦ 10 I^Ar ♦ ..... 



Ve will take only the first derivatives (froa equations 21) except irtiere tbej 
Taniah at r«0. In the present exaaple Ar'vill be chosen to be lO^cn, snail 
enough to give a fairly snooth run of the variables P, T, etc. 
In the first region, at r ■ lO^ca, designated r^^ 

ll(ri) - iijr ri^ Pc - 2 X 1030 g». 

' ' 3 ' . . . 

Since at the center /— ) »^ 0» it is necestaxy to take the second order tem 

Vdr/e 

In the series expansion: 

10)/ ■> [^'' ^ - ^ e • - P;- °- 

Taking « 6.67 x 10"°, its value in cgs units, and « 6O.5, ve find 

Htj^) - Pc - I (|^c^) '1^ ■ ^ * ~ 2 X 10I5 . 7.8 X 10I6 dynes/cm^. 

The luminosity at r. becomes 



By (11) (dL/dr)^ ■ 0 since it is proportional to r^. Again the 2nd derlvatlTe 

" . ■ u . • 

must be used: 



Ihorefon a third tmrm in thm •xpaniibn is nMd«d. It is 



\,{Ty) ■ ^lij ^^)3 • Si jgpg ^ proportional to r, r^, atoijAl 



vhiar« m bare used Ar ■ r, - 10^c», as befors. 

This star has a relativsly low Tc and tha p-p chain viU he the only im- 
portant mde of energy production. Therefore, ve can take for e at the center 

ftroM Kqn. (33), 

Taking ■ 60.5 gm/ca^ and T^ ■ 8 x 10^ ve obtain €^ • 2.U8 ergs/sec-gn, 
and L(ri) ■ x lO^O erg/sec. \ 

Since T has a maxinnxn at the center » (dT/dr)^ ■ 0, and a second term in 
the expansion is required, it is 

Therefore T(r- ) ■ T^ • "^ePc 'c^^l ^ 

1 \ r— 

Sac Tc^ 

The opacity ic^^ must be eval\iated before TCr^^) can be calculat|5d. Taking the 

Kramers value for fx*ee-free transitions (b-f transitions may nsglectsd whsn 

„ • • .... . • 

•Z ■ 0) »nd adding the electron sci&tterlng opacity (equation 26) gives 

K - 0.39 + 8 X 1022 Pc . 0.39 + 3.3»» - 3.73 cm"-^ 

Tc3.5 ^ 

With this opacity T(r3^) « 7.86 X lO^ ^ic. 

37 ■ '' 



In Table 1 •!! of tho aboT* Y»lu«i »t r • 0, Jr^ ■* 2 x ao^e* and at 
■ueeesslTe Ar - 2 x lO^cm Intarrali ara lJ.it ad: Tha appllcabla «truotura 
fbrmtilaa ara eollcietad balov in nuaerlcal ' form 
Alfi(p) - la.SfMf^Ar 

; •> 

AP(r) - -6»67 X lO^p'^^ Ar 

AT - -263 Ar (radlatiTa tpanaport) 

• ' -1.62 X 10-16 M(r) /s^j. (coDTactlve transport) 

r2 

- ,12:6 r2peAr 08) 

- 10-29 pT** « 

- .39 + 8 X 1022 £3,5 



AT 
AL 



" 6.05 X 10"^ I 



T 



< |0.U AP_| condition for stability against convection 
TAr PAr 

p ■ 8.95 X 10-9 p3/5 ,1.30 X 10^^ T"^*^ (convectiire traaspoi 
(HOTE: The formula for AT (cqnvective) coaes from ■ (1 - y) ^ vi*^ 

. ... ... ■ . ... - ■ - . - ^ ■ ■■ - ■ 

AP ■ ' 

:f > /a^^' irtiich is appropriate for a monatonlc gaa. Replacing ^ by 
^ QMCr)p ^4 noting that p ■ ^ gives the formula shown. The relations 
between p and P for conrectire (adiabatic) transjpiort is the familiar 
adiabatic gas law p - ftP - i^T^^'' with constants which are 

evaluated frrai particular valuos of (p, P) and (p, T) calculated In the 
convective region of Table 1 belo\f. ) " 



(lOlOci) 



0 



.2. 



■ .6 

I 

.8 
1.0 

convection .1.2 
l)eginB, , + ' 

• ' 1.1* 
1.6 
1,8 



TABI£1: Apurelvdrogen«odel»tar'iategritedfroir.. 



0 outwd., 



M(r) 
• (10^3 ga) 




.002 



.010I» 
.063 

. .109 
.255 ■ 
.^15' 
.82a 

1.13 



P(r) 

(10^5 dyneB/cs2) 



TO 

6?.8 

■ 51.9 
l»1.5 
1T.9 
-2.3 



T(r 



T.86 
T.5T. 
T.30 
T.OO 

6.30' 

3.51 
3.00 

2.31 
1.26 



.L(r) • 
(1033 erg/BM) 



M 
.060 
.ll»T 
.288 
.i(20 
.l»65 



Mi 



60.5 
60.0 
59.1 
58.0 
56.9 
5T.9 



83.T 



'(erg/sec-j 



2.i»8 
2.29 
1.93 
1.65 
1.5^^ 
■.,91. 
.Ill 

.0T5 
.013 




3.73 
3.92 

'>.35 : 

l».8l 

1.78 

82.3 
ll>3.6 
201 
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J Thl«. calculation has be^a^ carrlad only aa.far aa r - .l.e x l6^Pc» 
irtiara P(r) becomea negatlva. It Bust go to la^o Juat a bit short of this 
radius vhsrs^T chouid also becoms ssro as rsqulrtd by the surfacs bound^ 

. conditions. T&« fact that both P'4nd T do not slBultaneou*!^ bacQ«i>s«r9 

f ■ ■ 

' t ■■ ■ ^ 

for SOBS r which woiad then bs tbs star's radius, R, is a sign that tha . 

correct initial conditions vers not chosen at r ■ 0. This is no surpris'e 

sinee Pg an^^g vere chosen arbitrarily fox- the 'first run through the cal- 

culatioa.' A trial and error approach could be adopted, going back to r ■ 0 

and changing the central paraaeters Tg slightly and repeating until 

T(R) « P(R) « 0 slnultaneouBly. . There are more sophlitlcated methods for 

insuring a fit at the aurface^^^, but theae will not be deacribed here. 

Our model Is unsucceasftd in the sense that Our first guess at initial 
conditions did not produce a possible star. It did serve to ill\istrate the 
numerical integration of the structure equations considered as finite dif- 
ference equationa (38) rather than aa differential equationa (21). The 
Bodel alao illuatratea-the development of a convective envelope on a radia- 
tive core.^ Initially at- about r ■ 1.2 to l.U a Buapicioualy l?a*ge AT/Ar - 
waa found uaing the radiative tranaport equation for AT/Ar. When thia unas 
tested by the stability criterion (Eq. l6), itvaa found that convection 
takea place in this region. So AT/Ar vas recalculated uaing the convective 
for* of AT/Ar froB Equation (38). At every aucceeding atep thla teat had to 
be repeated to bie aure convection peraiated. 

Stellar modela are generally preaented in tabular form aa in Table 1 
atwve or in graphical form in which varioula parametera P(r), T(r), L(r), 
M(r) are plotted againat r. In Figurea 6 and 7 below, thia is done for our 
first trial aodel. 



^ .* /.a IX ' M /.i M 2.0 

Mass and lualnoslty yarlatlon vith radius for trial m^a^l star, 
The arrov oa'tho abscisaa shova tha position of the l)oundary 
between the radiative core and the convective envelope • 




O »2 .4 > 

Radial dependence of P» T,'i;i a trial Bodel atar^th radiative 
core and convtpctive envelope. The arrow at r « 1,3 x lOlOcn 
mBTkn the transition between' tliecie regions. 



'^In Pigttre J* the^bbund&^ betvttjt nldiatlvt aad qoj(iT»otlt!!r r«glont is ^ 

\ * ■ . ' • . ■ ^ , » i'. • . . * .• . • ' * . 

cleaae'ly indicated l>y; the marked levelling of f of the J9!Cr) funetion. Thie 

is to be ejqpected bebtuse eonTectioh iir sn, effioient Mans of ' energjr tnoii- 

port* and Qoiisequently co'nTective regibps haTS smaller temperature gradients 

than adjtfbent radiative*^ regiDns. ^n Figure 6 the oonirectiTe regioji is seen 

to lie at about tl^e point vhere the Ixminosity L(r) ceases to increase be« 

cau^e nuclear energy generation has ceased ('c « 0} vfaere the teoperature has 

beeone too lov to sustain a^ppreciable hydrogen burning. ; 

J- ' '/ " ■■ . ■ - \ ' ■ - . " ^ ■ ^ ■ 

The ■ass'of the trial nodel cannot be determined \mtil the model has 
been adjusted to give proper surface conditions, P(r) » T(R) « &• Figure 6 
shows that M(r) is still rising at r 1,8 x lO^Ov At this point ;M.» 1.13 
X lo33 « 0.57 so it appears that this model ^properly readjusted, vould^ 
yield a star of « 0,6- - O.tS" golar masses. Its luminoAityl\irouid be^« * T 
:o.U8 x\l033.erg/sec.or about' 1/iB L^,t-and its iradlus^vould be rbuijhly 2 x l(P-0 
cm, bfVabout Q.3 R^* With ';^8e raxxgh 'estimates of L and R, thc^ffective ^ 
surface temperature becomes « dO,OOO^K from the Stiifan<*-BoltuBann equation. 



/ This vould be a very blue, hot, 8ub«*dWfLrf star, veil beloiv the main; * 

■ If • • . " . _ 

sequence. an<l far to the left , in the region generally associated vljbh the^ 

. ; \ " ' ■ ■ ■ ^. " 

central stars of plcmetary 'nebulae. / ' 

^ * The general' features of the curves in Figures 6 and 7 &re repeated in 
the folloving figtures that illustrate real, detailed nibdel calqiiLations for 
■ain sequence stars of various masses^ 



l42 




8: braphlcai retrtientatioh of . • aod«l itar of 2.5 M^, R ■ 1.59 R ' 

• • « - . . . ■ ' * 

of "cpiflpoiltion X ■ O-.pO, Y ■ 0.09,. Z ■ 0.01^^),. L ■ 21 L^.for - . 

^ this upper main sequeneo itar . It liai eonv«etlv« core iiadc 

radiative enTolopa aad^ energy i4 geMra.ted oalnly; l^"^ the carbon 

cycle. The boundary betviaen^core and^<n^reldi>e is ait r « 2 z lO^'^c 




9: Model star of 0.6 M^; R ■ i^375 of conpoaition X ■ 0.72, 

Y "^0.26^ Z ■ 0.02<T), L ■ 0.1» for thil atar vhich haa a radia- 
tive core and^ convectlve envelope Joined at r =r .O.iBj x lO^^cma 
The fltiar generates its energy by the p«-p process. Its surface 
tenperature is 7300^, somewhat hotter and bluer than the Sun. 



ni« BOiUil star shovn in Figure 9 !■ rathtr J>lk« our tarlltr trial 
Bodal In si so and maaa. It la a good» tquili^rlua atar vhllt tha trial 
■odal vaa not . The higher central preaaure and tasperatxire auggeat tbi^t 
the trial nodal ought to be Modified In thla direction. ^ 

The aatrophyalca llteratur^contalna aany exaaplea of atellar sodela 
calculated for a vide range of naaaea and conqpoaitiona. Theae, plotted on 
the H-R diagraa and maaa-luBlnoaity diagraa reproduce quite aatiafyingly 
the dlatrlbutlona of real atara on theae dlagraaa. 

The field of atellar erolution irtilch la a genefallaatlon of the field; 
of atellar ^atructure determination actually inrolTea the calculation of 
vhole aerlea of nodela allovihga small coispoaition change betVeen auccea- 
ali'e Bodela of the kind that vould take place in a atar converting Ita 
internal atorea of into He and heavier elementa. Theae aucceaaive nodela 
involye a calculable time lapae'^^betveen them during vhich the compoaltion. 
Change occura. When the aeries is connected by a line on the H-*R diagram^ 
auch a^llne represents, an eyplutionaxy track on the diagram along vhich a 
atar mores aa it ages. The various ^forms auch tracks tsike depend upon both . 
initial mass and composition of the star, and they form the basis for a 
diacusalon^of atellar evolution. 
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PULSARS 

I« History 

P»l.«. vr. flr.t dl.c.v.„d 1„ 196T. «,r. or 1... by „cid«t. 

A gro^ Of r^o, „tro„oa,r. .t Ombrlag., s„gi„4, (H.ri.h. feu. PlmngtcB. 
8co.*f W conm.) (« W JuBt put into operation . n« radio t.l..cop. 
optlM«.d to a.t«t .cmtluatlona In r«Uo .our«. ..„..d by .^i 
wftlon. in the InterpW^ g„ or ".ol^ vlnd". Thl. ». p.rt of „ 

l.^rt«rt progr«, to m.«ur« th.-.u.. Of «trono«lc,a radio liourc.^ ' 
fro, .Ue «a3ur.m.nt. to d.t.rMlne tholr dl.t«>«.. to l^rt«t future 
Of th. r^o d.t,=tlon .yrt^ Of th. t,l,.cop. «. It. capabUlty to raspond 
to, r^Jd fluctuation, m. .Ignal atrangth on a luiaaeond time acal*. Thl. : 
r..p«,., had bean dealgnad Into th. ayat-.to p»nrlt « to record .olntUl.-: , ' 
tlona. but It alao gave If the capability to r..pond to puiai- radio ,*1.- 



8ion. 



Shortly after thl. nev teleacope >.nt Into operation, the 0«.brldg. ' 
group bega. to notice occaalonai .poradlc' lnt.rf.r.n.e. the .trongeet appear- 
ing to come from a location m the conatoUatlon Vul,«=ul.. Thla .ourc. ™ 
periodic and conaLted of .hort pulae. bf 81.5 MK, energy (th. f«,uency to ' 
vhlch the teleacope wa. tuned) repeated p^laely „e^ 1.33T .econd,; ai.' 
flrat Piaaar «a reported the following Pebruaxy, 196B. In Nature and va. ' 
conflr«d promptly ^ other group, at Jodr.il Ban*. Areclbo, OSIBO In Au.- 
tralla. and at JPL In California. MeanvhUe «yla and Balley(2) m«,e a more 
precl.. determination of the location of thl. pul.lng ^adlo aourc. u.lng th. 
Cambridge 1-mU. interferometer, ftey found tie radio aourc e to- be located 
near enough to a faint, aome.hat reddl.h .tar.ti .ugg..t .trongly that th. 
.tar and radio aource were one and the same. 



Pilklngton -jid his eoworkera^^) rechecked other locations where inter- 
■ - ' ■ • ■ • , ' 

ferehce.had been found among the aeintlllation measuraments and quickly - 

found three more pulsars . Meanwhile, worktri at Kitt Peak began searching 
for op^cal pulsations, and early in lynds. Maran. and Trumbo^**) found 
What seeae^ to be a regular but very small variation tn light from the star 
believed to be associated with the first pulsar, but these light variations • 
had * period twice that of the radio pulses and a different shape. . 

The search for more of these Wsterious objects was on In earnest, and 
to date more than 100 have, been found and studied. Probably the most valu- 
able of these from the standpoint of the help it has provided In understanding 
the nature of pulsars la the Crab pulsar; NP0532. This object was first dls- ; 
cowred in the Crab Nebula in mid-l968i and it has. the shortest pulse period i ' 
33 mllllsecdnds, of any p^sar yet found. The following January a group, of. 
aatronomer8t5) in Arisona (Cocke, Disniiy. and Taylor) made a momentous dis- 
covery when they fouid regular pulsations in the light from a central star 
in the Crab Nebula long believed to be. the remnant of the supernova that pro- 
duced the nebula in AD 103k. Thia optical pulsar had exactly the same period 
as the radio pulsar (.033095 seconds) and was within the region from which . 
the radio signals, came. Identification of the Crab pulsar with the Crab 
supernova remnant was important in two respects r it connected pulsars and# 
supernova residues positively, for the. first time, and it allowed theorists 
to draw upon a large body of information that had accumulated about the Crab 
Nebula which; has been carefully studied because it is the closest recent 
supernova and the one that [can be most fully resolved by our telescopes. . 



II > Obs erved PropertleB of Pulaars ^ ' i 

The highly unusual character of the pulsar Implied by its rapid asd ijp- 
credibly precise pulse period generated enomous interest imediately fol- 
lowing the discovery of the first of these objects. As a result, most of 
the world's radio telescopes vere. quickly put to searching for more pulsar, 
and studying those already found. The results of this observational activity 
have led to che following picture, of pulsar characteristics. : 

bursts of radio' pnergy that constitute the observed pulsds contain 
radl^uency energy covering a' very wide frequency band, in some cases i^to 
the opti<5al, ultraviolet and even x-ray regions of the spectrum. Each pulsar 
has its individual pulse shape but only when averagk over a large number of 
pulses. There 1b considerable .variation in shape from pulse to pulse frim a 
given pulsar. Figure 1 shows characteristic pulse shapes froo four pulsar.. 
These shapes are averaged' bver a large number of individual pul.es in each 
case.. The shapes also depend somewhat upon the radiofrequency at which the 
observations are pade. 
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Plkure 1. ' Average pulse shapes of four pulsars, The, ordinate is energy flux 
received on Earth but normalized to the same peak. value for each 
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_ . .The most striking and 8igniflcai.t feature of pulsar emission is the 
ieeularity of the pulse repetition rate. These periods range from 33 mllli- 

,aeconda up to k.B second, and the pulae period averaged over a large number 
of pulses is precise to a fey parts per billion/making them among the most 
precise ^clocks" known. Even a laboratory Vality quartz crystal oscillator 
can scarcely match this precision which is in the range of atomic frequency 
standards. 

A second important observation on the pulae 'rate is that all pulsars: 
are slowing down at measureable rates as detemihed by a lengthening of the 
pulse period. The rate of slowing down can be used to. deduce- a "lifetime" 
for the pulsar. The observed rate >of change of the^pui^Iriod has been 
found proportional to the period. P.Ja fact which can be expressed . 



dP , P 



(1) 



T»,i8 leads to a time dependence of the fom P(t) - P(0) e"*^^ which identifies 
T as. a lifetime, the: ^ime it. talces for the pulse rate .to decrease by e"*. In 
. this way it has been found-^at pulsar lifetl^^ 

10' years, rather shorytimes by cosmic standards. 

Two brief exce^ions tee the sUady runni^^ 
been observed. Several pulsars have abruptly exhibited a small increase in 
pulse rate, and^^ollowlng these "glitches" they have again settled down to 
reW their^ormal Slowing dow,. rates. We shall later see , that this phen- " 
omenon has.^revealed an intr 

'.Detailed analysis of individual pulse structures shows that several T 
things are apparently going on simultaneously in a pulsar, for there are 
three distinct tim% scales involved. There is' first the long time scale of 
the pulse rate, namely, .033 to k.B seconds. Then there is a shorter time 



scale, on the order of milliseconds, that accounts for much of the vari- 
ability |n shape from one pulse to another. Finally, the shortest observed 
time scale is on the order of 100 microseconds, and this time scale sets 
an upper limit to the size of the region that radiates the pulsar pulses/ 

The crgument goes like this: if a radiator emits periodic electro- 
magnetic waves at a more or less definite frequency, all parts of the 
radiator that contribute to this radiation must act in phase. If this were 
-not the case, different parts of the body would radiate out ^f phase with 
each other and would on the aVera^e c;ncel each other's radiation. Now if 
all parts of a radiator cooperate to radiate-in phase, each part must some- 
how receive information about the oihers , otherwise they could never ^in 
synchronize/in phase. : Relativity theorjr teaches us that no information can 
be transn/tted faster than the speed of light, c. Therefore, parts of a " 
radiato/ that are farther apart than an electromagnet ic; wave can travel in 
about a period of the radiation frequency cannot stay synchronized because 
they, are effectively out of touch with each other on the time scale of the 
radiation endtted. Thik means that a radiator emitting, regxaar pui 
intervals of l£l'0-y sec must not be larger than light can travel in 100 w sec. 
This a^ts an upper limit on the size of the radiator in a.pulsar , and this 
size limit is .about 30 km (lO- sec x 3 x l'0» km/sec). This does noV necesi 
sarily mean thkt the" entire object that constitutes the pulsar cannot be , • 
larger than this; the radiator might be a "spot" of this size on the surface 
of a larger body, and additional information is needed before a choice can 
be made between these possibiiities. 

Another important observation is that the radiofrequency energy in the! j: 
pulses, is polarized. Polarized radio waves have the samt-g-enerarproperties 



as polarized light, and „radio polarization is recognized by ising receiving 
antennas that respond to polarized incident waves. Polarization in pulsar 
radiation was found to be very complex. It"depends/upon freq4»ncy and Is 
highly variable, changing fro. pulse to' pulse and even within ^ single pulse! 
But the Significance of the finding that the radiation is Btronily polarized 
18 that it vexy likely originates in a region in which tkere is'^strong 
nuignetic field because we kpow that radiation' frop regions permeatSd by 
n«.gnetic fields is characterized V strohfe polarization. This. conciLion 
is supported by the observed spectral distribution of the radiation an^- by 
the fact that the crab Nebula contains strong fields, as detennined ^ ibde- 
pendent Observations Of other features of the Nebula.; Tl»us it se.ms-^^r 
^.^est^lished that strong nu^etic fields are associated with pulsars.#^ ' 
In order to deduce the amount of energy, radiat^cf ; by a pulsa*. its^ dis- 
tance must be known. One of the important consetiuences of identification 
Of the Crab Nebula with NP0532 is that the distance of the Crab is accurately 
known, and this permits a goo.d determination of the energy radiated by the 
Crab pulsar. 

■iir on. .m.thoa-or «,oth.r all th. >mo«. ,ul.«.s hav. .^.n-.-fpSd to U. 
-Ithln our Oolaxy. Pulsars undoubtedly «i.t 1„ other ^^^s, y^.^-' 
ore too to u obs-irvabl,. The distribution of pulsays^-„<^e aalaxy 
is consistent .1th. their Identification as supemova>4ants. ku.nu^a?' 
..ISO .upports.thls vlev. ^i'ce U^y^rijel^ti^^ .. 
(< 10' years). ther^uld-not-he-W^r^^ to observe at w given tl.1; 
if,they-«.e-«,e consequences of , supernova explosions . their number should be ^ 
correlated with the observed fr;,u^ncy, of supernova occurrence, an. this turns 
out to be the case, allovln«. for'a number of obse^atlonal uncertainties. ' 



^ The enerar loas r.t. frc. a ^ ^k. th. Crab ha. been foind 

to Of th. order Of lo", erg, ^second, l„ei„,es radiation over- 

th. enttr. electro^agneti, .pecWind Is approximately lo' tl«. g^eit., 

than -the. r«.latlon .W tj-i Sun, If an object radiates at this rat. for a- 

lifetime of = lo' years. It .mist lose to. that time a total of => lo" e^g. 

W. Shall us, thl. flgurc later to W what pul^^r models 'o^, lontaln'thl. ' 
much energy. - - 

lit. Basic Pulsar MnH#>1c» 

we too, only three general types of massive- objects In nature that are 
good timekeepers and thereVe possible .candidate; for models of pul.«.. * ■ 
.hose Chief Characteristic Is Its remarkable regularity, aese classes of 
Objects that are capable of relatively slow periodic motions are vibrating, 
.tars, blna^ stars,:™d rotating stars. L^t us .„mln. each brlefly ln ' 
light of observed pulsar characteristics. 

. Certain types of stars are knov^to pulsate a. a result of radlar^^^T 
latlons Of their atmo-spheres. The P.rlod. r.ng^m-U^r;r;:;^ths and the 
periods do not exhibit th^^rd^f,p„«-ron of , pulsar pulses. The genera 
theonr Of vlbratl5g,bod«ri«ds to a dependence of the period . P upon the . 
d^BWy^Of the form P - p% or upon mas, of the form P . n,ese two 
Simple expressions pemit u. to rule out vibrating stars as pulssr models ' 
because the second predicts the ^ng variation of pulse period .1th time. 
As time goes on the mass of a star ™st ordinarily deerease.^^ It loses ' 

by matter ejeetl™ in the form of a stellar .Ind. ,L it radiates a«y 
.»as m the form of energy (E - Mc'). Consequently If M decreases, the vi- 
bration period Should decrease and the pulse rate ought\to increase rather 
than decrease, as observed. This ardent is not by it.Wlf conclusive. 



toweTer. because stars may also accrete matter under certain circumstances 
by gravitational attractlon^f the gaseous enrelope from a close binary red 
giant cojnpanion, for example. The accretion rate could exceed the mass loss 
by radiation «id steUar wind and could lead to a pulse rate slowing as 
observed, But ^en 6ne calculates the densities required to give pulse rates 
of a second or less, a problem crops up/ The most dense* white dwarf stars^ 
are not able to vibrate rapidly enough, th fact, unless one EPStariteTrather 
unlikely complex modes of vibration . likely to be^bon-^^niel out . it does 
not appear j^ossible to devise a plausme^vib^^^Sng white dwarf ^m^^^^ with 
a period shorter than = 5 secpnds'r^a is close enough to the \ . 8 sec . 
period of^«^e^lo5n»sr^ 

JCrab)^to .be explained away by uncertainties in the calculation. Only a 
-star i^lth density in the ran^sejof nuclear densities (> 10>» ^/cm>) could 
vibrate ^aa rapidly as. the Crab pulsar. 

A second general group of candidates for pulsar models is the clwis of 
Closely spaced binary stars revolving about their com«,n center of mass much 
like a rotating dumbbell. Let us see what kind of binary would be required 
to rotate 3p times per second. The laws of orbital motion {Kepler and Newton)' 
tell us thai-the masses, periods, and sizes of the orbit of 'the -members of - 
a binary pair are related by Kepler's 3rd law 



(Jt^ + Mg) P» - hif Ga» 



(2) 



vhere G is the gravitational constant and a is the semi major axis of the 
elilptlcal trajectory^^the-members of a binary are comparable in mass 



to the Sun. then using P.. .033 seconds we find a . 2 x io' cm. about ^ 

• \v 300 
tb, raaiuB ofthe Earth, Thus a M„a^ „od.l Wd 'require appSoxWeay.' 

solar i^ass objacts to ravolv. In orbits Much smaller/than the radius of the ;, 
Earth. If theee condensed star, were In contact their Individual radU could , 
not exceed 2 X. 10' c and their densities would -be of the order of p ^ M /i,.' 
_ or about 6 x 10" g./c„=. Again ve find that a stellar object vlth a denslV 
Of the o^er of lo" ^ c.- or higher Is required. This Is the lover UMt 
for the density Of a .eutron star, a body In ,Hlch nearly all of the Matter - 
normally m the fon. of protons and electrons has been converted Into neutron, 
by the enor,ous pressures and ten^ratures Involved. Th^ neutrons are fomed 
by the reaction . e-V n, and the protons ind. electrons come froM the 
matter of the star. >alnly hydrogen and helium, that ha, been c6,*let.ly lonl.e< 

«e ,ee that whatever -th. 'defects of the vibrational or binary star ■«,a.l. 
"ay be. both call -foj object, of denMtle, In the neutron star range, w. ■ 
^shall soon see th^ the rotating object n,odel also, de»^d, a neutron' , tar . 
so.the decision is „nanl„„s. But before the rotating «.del: li con,ldered. ' 

°"" '"-^'^'f^ts of"the binary Model of a pulaar. The 
%4^isf:ihe,e is very slnaple. If the orbital kinetic energy of the ^in^ ' 
-is What .feed, the pulse energy, as the binary loses energy.lt MU,t ,peed up.- 
not Slow down as observed. ^1, follow, froM Kepler's law above and tht fact 
that the orbital- sue. a, ^st decrease as energy 1, extracted froM W systifc. 
in this respect the blna^ «odel has the- sa»e drawback as the vibrational- '/ ' 
»odel: both Should speed up rather than slow down a, tl„ goe, on. There 
is another ^re esoteric defect of the binary Model. Oeneral .relatlvlty -theory' 
ssys that two ve^., «assive bodies 'Orbiting close together,,, the Model demand. • ' 
are ,tr.ng radiators or gravitational waves. These wave, cal; energy away 



from the binaiy. system, and calculation reveals that a ;33- milliBecoi{r'i,inkry" 
would be. such a strong gravitational; wave source that' it would lose' all its 
energ^ in a rerWveiy short time/ far leas than 10' years. " > V " ■■■ 
Finally, we tuni to. the third possibility for a pulsar model, a Totiting 
body with „a radiating source soaehow attached to its surface so the radiation 
aw^ ki^'und in space like a lighthouse beam!. The pulses would be received 
•on Earth each time the beam sweeps by. and the pulse rate would give the.ro- 
tational period of the object. 

.^A simple argument can be used to set an uppei^ limi^ on t^e sizT^such 
a pulsar model. . We require that its surface gravity at least equal the centri 
fugal. force at its equator; otherwise it would lose matter rApidly from the 
equator. If the tot^. mass and rkdius of the object are M .and R respectively, 
^the gravitational acceleration at the surface is simply 

" GM ' . ' ■ " ■ 

■ . ^ ". ■ ' , ■ " ■ (3)\ 

^ vhere G is the gravitational constant (6.67 x lO"-' cgs units). . At the equatbr 
the centripetal acceleration is Ra,^ where a. is the angular frequency <,f rota- 
tio9. For ssero mass loss we require • . . 



R 
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The ^ss'niay be written in terms^ of the average density as M » 1 7tR» jr. '^Also 
vritlng a, = 27r/P allows (V) to be put into th<r more-simple form^ 

" . /'^^^"^ '(5) / 

This, establishes the minimum average density for the slowest pulsar at . 

= 10' gm cm-', in the vhite dwarf range, and for £he Crab pulsar p := 10»V g«cm- 



:!:\T " " - - - »-n-.w .... 

wext, the rotation*! ■ ^mge. 

«ar o, . ^^^^^^ ^ ^^^^ , 

E ■ 1 lu^ '■■ - 

2 ■ ■ ■ ■ ■ ■• ■ -• ' ■ . 

(6) 

"ipated whenju qhangea la given 

■" ■ ■ -■ _ ^ .: m] 

ror the noment of inertia I « mr^ v. " ^ 

.bove involving aght tr.v.l „ . «" ?on.l4.r.ti„n. 

, «u rotational energy E a 2 x 10*" P-a - 

E * 2 X 10" ergs The ^ • ^ " ' ^33 sec. 

Of. the nrtula (= 3 , 10" . 7" ' " H th. 

, , ,0" .,g3 in r " ""^ *° ^' yP-xi- 

- 1° tM» partloulM .pulsar TM. . " T . 

atar case. •> ; ^ ^ - iiretlme as In the binary 

- * neutron star' la to ^ 

' * small Integer factor n - ^ 

■ ; , • , ' - ■ » n»ybe there 



». tvo or rour pul... p., ^^ ^^^ ^^^^ 

.t„ conaerv.. ^^^^^^^ coU.p.e rotation 

»t. „eo„„t.d for . copanar . .t„ u*. the 8u» vlti, ^ . T , 10» c, 

ani rotatlona period » 21 days (. 2 , io« ,„ , „ , . 

■V V £ I iu BOO ). If Buck a body were to 

contract to R . ^ . =,=o„Bervl„« .^ar ^„t« . r,. 

t-tlooal period «,^d decre... ^ a factor (B./B^,-. . „ , 

5 X10 In tMB^c^^ie. noal period would th«, be . 0.5 ^XUaecond. " 

~r. than abort e„ou«b for a ,.w p„..„. „ ' 

uncertainties in this crude calculation 

e caicmation. the agreement with observation is 

satisfactory. 

. It va. ^^^^ '" Vfew pulBar. there have beeo'abrupt 

■ chan^ea i„ p,,iod.U.din« to the concXuBioo that the pu^aar BpeedB up 
«d then .ooo aettie. down to it. fo„e. rate of .W„, down. The., .peed- 

<iP» or rtitchM are quit, .mall/ tut ^. 

waaurau.. Nhat can cau.. a 

r^tatins ohj.ct to .peed up if no ..t.^ tor,u. i. .ppu.d, Xh. .n„.r 
i. .i^l.; it. radiu. d.cr...e. cauain, a decree.. i„ it. ^ : 
I; « . oon.e,uence. th. rotational fre,u.ncy,„ «oe. up and the p.riod de- - 
crea.... ^. ^^ .^ ^^^^^ .hould .udd.nly decreaae 

U. radiu. by about 1 c. the a^unt needed to produce the obaerved .peed-up . 
vao .olved by Rudeinan'^' who Bhowed tw 

iwer Of a neutron .t«. 
™. probabi, .oxid and in fact, cryataiiine. strea.e. in thi. c™t fro. ' 

centrifu^ „d Coriolia force. occionaX^ cauae it to crac. or to rec^. ' 
.telli^e. producins a ..ata.,ua.e... .he reauXtin. re„r„«e„e,t of the .t.Uar 
•urface l.ada to a te^ra^ ^ ^^^^ 

lncr.a.e in the,at..of rotation. Thi, ,atch oniy Xaata a f^ wT... until ^ 
the «.r can rcu™ it. former .tructur.. The =^.talUn. natur. of a neut„n 



exp}.alnlng the peculiar , sudden apeed-up of some piilsars. 

If we accept the factiimplied by: the Blmilarities betveen pulsars of 
videly different behavior^ that all pulsars are similar obJectJs and that the 
long period ones are simply older and have slovtd down more, we pust conclude 
that initially they form as neutron starp with densities > 10** gm cm""* and 
with radii ^ 10 km rotating at periods as short as a .020 sec, the extrapolated 
period of the Crab piulsar at its formation. Moreover, these neutron stars 
must possess enomously strong magnetic fields, of the orde^ of 10*^ gauss 
or more at the poles. , This est iniate fallows from the reasonable assumption 
that magnetic flux is conserved as thie.Btar collapses in the supernova event 
lihat forms the neutron star. A star collapsing from the size of the Sun 
(R^ a 7 X 10*^ cm) to the size of a neutron star (R »lp* cm) undergoes a change 
of radius by a factor 7 x 10* and a change of cross section area equal to \ 
the square of this factor,.or « 5 x 10*. If the stiur had an initial polar i 
field of the order ofxonly 2 gauss, the neutron star field would increase to 
a 10*" gauss. Many stars have much stronger initial fields, and if they are 
appreciably more massive than the Sun, as theoiy says th^y must be to become 
supernovae, then the final field may well exceed 10*^ gaussw The presence 
6f a strong field has several consequences of importance for the pulsar. It 
provides for polarization of the pulsar radiation in agreement with observa- 
tion. It also provides a possible connection between the. neutron star sur- 
face and a region some distance away from \rtiich the observed radiation may- 
actually be emitted as required dn some theories of pulsar emission. The 
field may also provide a sweeping beam or sheet of radiation to explain the 
reception of precisely timed pulses, y 



^ When.veget down to details of poealble mechaniama that might prodMpe 
a beam or aheet of eiectronagaetlc energy rotating with the atar. it i. no 
longer ao aimple to aingle out one of a number of competing theoriea. ' We 
, ahall l^ve to be content to outline the baaic requirementa of a aatiaf«.tory , 
theory and then to examine three different type, of theoriea. brief ly pointing 
out the advantagea;and weatoeaaea Of eac^ It ahould be .aid at the outaet ■ 
that no fully aatiafactory theory of pulaar electrodynamlca haa yet been 
developed, and it la not even certain that any of the three to be diacuaaed ^ 
here will prove in the long run to be c6rrect. though each aeema capable of . 
explaining 8(Mne pulsar feaiures . 

Th. -UTo, fl,u. eltted by . pul.„ cm b. cremated i„ , .tr,agWfor- / 
«rd fuihlon. «,a It give. „ i^,^ clu. to th. kind „, ^Siatlon m.ch«i« 
to be sought. If a pulsar Ilk. th. Crab «lt. pn th. or*.r of 10» ,rg./..c 
from « Object or radlu. = 10 to>. 1,h. flux -ritted b. at l..,t lo" .rg. 
o«-' and probabljr «r, like 10» erg „-» ..c-. .u.. the pulear do..„.,, 

radiate uniformly over It. *ple .urface. . A black body radiating thli flux 
Would hay. to have a t«.perati^ =■ 2 xlo' <* (.Inc. Plux . o T»). Thl. lo,^ 
po..lbl. until on. note, that the .h^- of the .peCrum i. definitely not > 
that Of a black body at 10' "K. If on. =on.ld.r. th. total .nerar in the ' 
r^lo .p.=tru«^ (. 6 X 10" .rg./.ec for the Oreb pul.ar) and =al=ulat.s the 
black body te^eratur. r.,ulred to produce this flux at ridlo fr.,„.ncleB. 
the. re.ultlng Wrature 1, eno»o„.. about 10» Ok. Only a coh.r«it radia- 
tion sourc. cs„ possibly produce such radio flux... Thu. on. „.st s..k . 
^chanls. that ^Ul Produc. coh.r.nt osdilatlon. Ilk. thos. of a tuned riilo 
tr«,a»ltt.r or las.r or «..r, n,ls cohsrenc. re,ulr».nt puts a strong " . ' 



reBtrlctlon on possible radiation mechanisms » but tlxere appear to be at 
least tvo ways to satisfy this requirement , both pf which depend upon the 
presence of a strong magnetic field. 

One of these mechanisms is synchrotron radiation that results when rela- ' 
tlvliaitlc electrons are injected into a strong magnetic field. The electrons . 
follow helical paths about tho^magnetic lines of force and radiate at the 
freqiiency with which they spiral about the field lines. This is classical 

radiation in the sense that classical electromagnetic theory requires an 

■ j-{ 

accelerated charge to radiate. Naturally, the revolving charge suffers 
centripetal acceleration in its helical motion about a field line and radiates^ 
accordingly. The radiation is moderately coherent at the circular frequency 
of the electron. The latter slows down as the electron radiates away energy, 
leading to emission of a broad band oiT frequencies with the highest frequencies 
in the spectrum being proportional to the magnetic field strength. A field 
of loV gauss will produce frequencies up into the far infrared while fields 
of 10** gauas will give frequencies into the hard x-ray range (X a 1 A). 
Synchrotron radiation is polarized in the plane in which the electrons spiral , 
and relativistic effects cause it to be emitted anisotropically In a narrow- 
angle disk about the orbital plane. These features are shown in Figure 2 below. 
(See next pa^e) 
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Figure 2. (a) Helipal notion of an electron In a magnetic field. 

(b) Synchrotron radiation of an orbiting electron. 

(c) Spectral distribution of aynchrotrbn radiation. V lathe 

^<iuency, the orbiting frequency of the°electron 
- about the field lines. . - 

The other type of coherent radiation to be considered involves a oaaer ' 
or laser type of emission but of a most unusual kind.^''^^ Consider the 
radius of the circular or helical motion of an electron In t magnetic field. 
This radius is inversely proportional to the field strength, and if the field . 
gets strong enough, the orbit rkdius can shrinictp atomic size. When this 
happens, the orbit must obey quantum conditions Just as It doas in an atom. 
Once this quantum regime is reached, the electron no longer behaves classi- ■ 
cally. and in, particular It no longer radiates synchrotron radiation. Now. 
as in the hydrogen atom, radiation can- be eaitted only when an electron Jub^s 
from one energy level to 'another. The levels can be calculated/ as In the 
hydrogen atom, by requiring that a stable orbif have a circumference that ia ' 
an integral Cnwnber. of dfeBroglle wavelength, of the electron. Detailed analyai. 



along these lines shows that the energy gap between neighboring orbital states 
Of an electron is ll6 H/H KeV where H ■ k.k x 10^' gauss, the value of field 
at yhich the magnetic enerc^ of the orating electron equalp its rest energy 
{0.31 MeV). At low fields <<H^ the energy levels of the electron in the field 
are virtually continuous since Ihey are very closely spaced. But, at high 
enough fields, H « il , the levels are widely spaced, and the spacing aay be 
large compared with kT, In^his case, the. excited levels are only sparsely 
occupied ty thermally t it ed electrons. The system is much like an atom at 
low temperatures. If thefe is an efficient non-^themal means of populating 
the excited jata^eia of thls s^^ then the possibility for laser action ezlits, 
and this 'Wgneto-laserV- can in principle produce a strong, coherent bean of 
radiation aa required, 7; The necessary population inversion to cause laser 



action is postTilated to be caused /k^ electron collisions among energetic elec* 
trons accelerated by electrodynamic processes in the surface regions of the 

o ■ . ' ■ 

neutron star where there is a thin layer of "normal", non-degenerate matter* 
Ope such acceleration mechanism will be discussed below. 

The rather exotic magnet;ic laser mechanism described above meets the 
, requirements of pulsar emission in respect to brightness temperature, cotaer** 
ence, and ani sot ropy (beaming). T^e magnetic field also provides polarisation • 
Electromagnetic radiation from such a luer is able to escape along the field 
lines, so it beams out from the magnetic poles. If the magnetic poles of 
the neutron star are not aligned with the axis of rotation, the radiation 
beam from the poles will sweep around in space like the beam from a lighthouse 
to give the characteristic pulse of the pulsar as seen by a distant observer. 
The figure below makes this clear .^JLt also shows a case in which there are 
two p\LLses per revolution of the sVar (Figure 2b). 




^•9^tiit. Fitiibnts 



(a.) 



Figure 2. "Liihthouse" model of pulaar with t.\ i i 

PU1.„. lnf.ct "lnt.rpul.„" Which „. w..k.r p„i... 

th^ »«in imlaes. The.. mi^b« 4 „„„ireBtatlon of ng^ 

-__^On.-of- th. .arll,.t thXjj^of pul'.„ r.dl.ilon «. .dv«,c.a by aoid'» 

itat. with th. atar a. In Tlgar* 



1*0 proposea that the magnrtlc field 



? Oh«-g.d partlcl.. fr„„ th. .urfac. of tb. .tar tl.* to the 

field lln.,, 1,.., they rotate tightly about the«. , a. th... particle. «>« 
»«y fro. tfc .tar along a given field l^„. t^^,^, .j^, ' 

tangential velocity of a fl.ld lln. l^crea... with dl.tab„ the of 

rotation according to th. f«ailar r.latlon v^ 4 ». (s.. Pig„. 3 
« a sufficient dlatauc. thl. t^ientlal velocity approache. tO. .peed of 
light and rf,.n thl. happ«.. th. partial, radlafi «.t of It. aicuMulatcd ' 
Mnetlc Wgy .l.ctro.«gnetlc radiation. The entire proc... 1. a bit 
Ilk. the gaae of "cr.ck-th.-.hlp.... .h.n the .Lctron .Ing. around fiat 
..nough, it let. go ln.bur.t of;.l.ctro«ig„.tic relation. The apeed of 
light 1, attained on a cyllbder of radlua r. - co,- which 1. about 1.700 
Ion for the Crab pulsar. 



r. 
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^Sp««d of light 
cylindtr 



Figure 3. Gold's Model 



The radiation at the speed of light cylinder will be somevhat beamed, as ' 
shown, by re. tivistic eff«etfl.,^ The "besia" will actually be a continuous"' 
dlsfc of radiation' radiating in all directions from the equator of the cyllhder 
V aa the situation Is sketched in Figure 3. This would not produce pulses at 
a distant receivei^. To do that, it must be postulated that electrons are 
emitted mainly either from localized "hot spots" on the surface of the neutron' 
star.or from its polar regions. These then foliowr the local field lines out ■ 
to the speed of li^t cylinder where they radiate from a localised region 
on the cylinder. This regioiT gives in tum'a localized beam;of radiation 
which races around, the cylinder once each' turn of the star. la this way the ■ 
model pan give pulses . /The existence of a fairly strong field at the cylinder. 
(«10» gauss) provides polarization of the radiition. The polarization, 
frequency distribution, bean width, etc. . are all characteristic of synbhrotron 
radiation as would be ej^ected since the electrons are being spun about in - | 
a gi«mt circle by the j-evolving field lines and the electron orbits are siaUiarl 
to those in a synchrotron./ • r 
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One difficulty with the Gold model whijjh ia conmon to all models that 
propose a source of radiation located some distance away from the neutron 
star surface is that it is hard to see how the relation between source lo- . 
cktion and star surface can be maintained so accurately over a long period 
of time as to give the observed pulse rate precision. Any drift b«iw^en the 
surface and the distant rmdlating region ^would spoil the preclBe repetition 
rate^of the pulsus. 

Gold's model calls for relativistic electrons (« 3 MeV) to be delivered 
to the speed of light cylinder. This requires cdnaiderable acceleration of 
the electrons between the neutron star surface and the light speed cylinder. 
Two mechanisms that have b««n proposed to accomplish this acceleration will 
be described next. / 

Goldreich and Julian^^O) suggested that the familiar mechanism of uni- 
polar induction will operate to accelerate electrons. This can best be d«- 
scribed In terms of the ele^ntary laboratory demonstration of the Faraday 
disk dynamo. CbnBlder a conducting disk rotating in a magnetic field as 
shown in Figure U iwjlow with slip rings on the axis and the edge of the disk 
to make electrical ^onnectldns. 




Figure 1». Slmjile^ Unipolar Generator 
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An electron of the rotating aiakls carried acroaa magl^tlc field lines 
vi.th veloQlty V - m and experiences a radial force F,\ - evB - erB... 
This force Is equivalent to a radial potential gradient U dy .or 



jiv;. 



flu rdr 



^.'-'l BUR* 



(8) 



^ere V Is the totia v|ltage developed across the disk of radius R. 

The important 4a|ure Of this device f or the pulsar application 1. that 
if electrons in the lo|zed, gaseous plasma surrounding the rotating neutron" 
^star are not rigido; att^hed to the «,tatlng magnetic, field initially, then " 
the field lines do pass i the electrons with some relative velocity and 
|hereby Induce an electri,. field, ^e 1^ x that can accelerate the electrons. 
Detailed consideration of a |lpole field shows that this mechanism is capable 
of acceleratlng0jc^;.pn^ an|l6ns away from the surface of a neutron star ^ 
to energies W9li;f.:lnt6^^ range, f^ aore than enou^ for the . 

Gold mod^. The electric |llii4nerated in this manner is strong enough 
to separate the electrons aj^d^i^is in the plasma near the star.^and tiL bharged 
particles mbVe alon g magnetic field lines as the ' diagram below sho^. : 
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«... corctatln, ^^,tcsp>,.„ is .ad. up of ohar..d pa«l.l„ atta.... to 
• Clc... ^.uo u„.,, ^a «au una. t^t „oss th. spae. of u«.t 

. Clindar ara opa„ at t>» a^faoa of o.Und.r tt... =lo.. fa«..v out , . ■ 
«a parties.. t.at ..aiata at t... U«^- aurfao. foUov t..a. opan fial, 
a„aa. ...ctrons an. proton, folio.. «ff..e„t seta of Una. ,aa a coo- : 

..quanoe of thalrdlffewnc. Inmaaa." 7 ' ■ 

It can be seen from Plffure S +ho+ -t-u^ j 
, , V ■ ^ ^^'^^ * region on the light cylinder 

Uhere the particles intersect It 

V ^ - compriseB a source of radiation which.5 

sweeps arou„d ^ththen.a«netic:pole as the star . otat;s i^ie ^ 

^oTia-ces, ^ This provides the >^ 

«,ui...,aarc,ai«.t affaot to p„.uc. o.a.^. pul.,. . . an. fa radiation " 1 
is Of t>» a:™o,.„tron t„a „ or «aar f. li«,t cninder aa in Oold-, .oaai. t 
TH. Ooldr.io..^ian ,t..o^ ia an i.p„v..nt . ..cauaa it a.o.a t^ ^ 

. Plaa^ and ^atoap.a„ ^t ..iat naar t.. aurfac. of a nautron .tar in- ^ 
.tead Of a vacuum. „d tHa- tHaor. fun,!.,., ^ ^^^^ ^^^^^^^ 

_-is-Pl^^.^ .e acoelaratad ava, fro. tH. atar to t.a li^Hapaad o.liW.. . 

/ diffarent. aooaW.tion .achanls. ™a propOaad ^.^unn and Oatri.ar(U) ' 
vHapomtad out t.at a-rotatin^ .a«natio dipol. . aourc. of lo. fra- 

.^r r.diatio, t.at could acoalarata .l.=tron...to,.i«. W.r«i.a. . rotating 
magnetic dipole will generate in its vlcinitv . 

^ . ^''^""^r a time-varying magnetic field': " 

In itself such a field can aeeAl#»T.o+- * 

hetatron. Basic electromagnetic theorr shows th^t -'^^^^ 

^ a varying magnetic field 

generates an electric f lei h p « «iTi/jx ' r 

At dist - . - ^ t^^^'^ield accelerates the particles. 

At distMces beyond the lightspeed cylinder th<o " 

an electromagnetic waye for relative 

. relatlvlstic reasons. . o„e might make the simplified 

argument that the field lines wh4ov, «- - 
. . . . - ^ «n«rgy cannot move faster t^^^ 

speed Of, light, and when they reach that velocity thev n \ 

^ velocity they no longer co-rotate -but 

-veaway.at the speed Of light as>electromagnetic waves; 



J. 



Cl«,lc«Uy th« >bt.ting akpole,raaiat«rat tli^ • ■ ' ' 

lie'- . ««"'•«■ , / ■ (9) 

Bi i. th. Meld .tr.»8th .t th. poX. of the .t.^-, H l'. the r«Uu.. 
«a th. «t.tio„.l fr.,u«cy. Taiing b„ . lo" gau... R^. io ta. ' 
^ ^ . 186 corriapondl-g to . p«,iod of . 033 .« . th. total pp«,r radl.t«l 1. 
3 X 10" .rg./„c. approxtBatoUr the total po«r onput oftli* Orab'pul..,. ^ 
th. rotating dipo^. 1. capahli of th. «,ulr.d .n.rgy o«tpat. ■ 

Chargad partlcl.. r , acc.l.rat.d to the «av. valocltjr th. itir 
f«. «.d th« rid. th. .l.ctro«*5„.tlc vav. out to th. light cyllnd.« plcklbg 
up additional ,n.rgy frp» th. vav. m th. **oi.„;; Th. "««„ riding" «.h>u.l.. 
1. .i«ll«- to that lo a llnaar accl.rato^ and 1. th. .I.,tr6i«gn.tlc analog' ' 
Of Burflng on an ocean vave. / 

•n..r. 1, a prohl„ ^ich plagu., th«,rl.. that po.tm.te acceleration 
of partlcl.. near th. .tar foUo«d ^ transport to th. light cyllnd.r h.fore 
?~llatlon take. plac. 01e.»lcally/ ,ynchrotron radiation «lll tak.; place tl 
th. accleratlon r.glo„, and a ckculatlbn ahow. thl. proce.. .convert, the 
energy gained hs- an electron Into radiation ao ■.fflclently that the electron 
will not be ahie " r.ach r.latl,rlstlc .p..da, ne,d.d to xeach th. llgh; cyllhd.r 
hy either .,f th. .ech«,l,™ outlmei ahov,. At flrat gl^c. this W ..e« t 
to he Of no l^rtance: Inatead of valtlng ™tll It. reache. the light cylinder 
an electron «l^ly radiate, prematurely; a„d ™ net converalon of ^tatlonal- ' 
energy Intp relation 1, pretty much t^ ^. Howerer. In a den.e cutting ' ' 
medium .elf aW-^Uon (i.e. . reahsorptlon of emitted photon.) ta..e. place ' V 
he-fore «e , addition c«, e»cipe,from„.the:me41um, and the resultant fluxicannot 
«cc,-: that Of a blaclchody ^t the te»p,ra^n^ of the^ medium ibout bne phoW 



mean. frtfejt>ath down into the .medium fiHDm its "aurface". Since the plasma ^ 

a neutron ator cannot exceed a temperature of the oyder of lO'oK.A 
- the escaping aynchrotpon radiation flux cannot exceed the blackbcdy flux from. ; 
ar-souree at that' tempe rat xure. ^But the meaaurad radiofrequoncy flux from a 
^jpulsar haa a brigfetheaa temperature at least 10** Or as "pointed out earlier, 
^ hence the discrepancy Introduced by self-absorption. \ 
. One way- out of the self-absorption dilenma is offered by the Basnet io 

\ laser mechanism of 'Chiii, Canute, and Occhionero, mentioned earlier. Two 
^factors act to eliminate the problem in this theotyV When there is a p6pu- 
, lation inversion with the upper energy levels of the system heavily pppulated,' 
V t^e -^system is no longer, in thermodynamic equilibrium and the blackbody liasit 

oi^Y radiated flux no longer, applies. It is well known that the brlgJitnsas 
■ '''ix \ ■ * " " ■ ' ^ -. ■ ' 

-tenii^^^ture of laser emission is extremely high, or to put' it another way, 

by emitting its radiation coherently a laijer is able to exceed by many orders 

■ ■ • \i ■ • ■■- ^' ■■ ■ .■.■■"„ 

of- magnitude the radiation output from a black body operating at the , temperatun 

of the lasing medium. The other factor mitigating the synchrotron self-absorp- 

tion problem is that, at the high fields under consideration the. electron orbits; 

. about the field lines are quantized as mentioned earlier, and they do not 

I* radiate ayhcKrotron radiation while in one of these quantized state a. More*- 

i over. If the field is high enough, the magnetic quantum states are far apart^ ^ 

:\ in ene3rgy,and only synchrotron radiation of very short wavelengths could be. 

I absorbed by such a quantized system. The longer wavelength radiation, ir 

"I present, would travel out of the magnetosphere without absorption. - 

\ Only a few of many theoretical proposals for a mechanism of pulsar radia- 

•^ion have been discussed and these only very briefly. The subject is both 
' \ ■ . . ■ ' 

complex and difficult, and to date no fully satisfactory theory has emersed. ' ' 

\ • ;• ■ ■ • ■ . . 



Hum ttqm its "eurface". Since the plasma 

Dt exceed a temperature of the order of lO^OK.^V 

Lon flux cannot exceed the blackbcdy fltix from. 

■„'.■■•■■'. •■ « 
3ut the measured radlofrequency flux from a 

;ure at least 10^*^ ^ as pointed out earlier, ' 

1 by self^abeorption. \ 

sorption dilemma is offered by the magnetic 

> and Occhionero, mentioned earlier • Two 

)blem in this theory^ When there is a pbpu- 

energy levels of the system heavily populated/' 

lodynamic equilibrium and thta blackbody lladt 

ea. It is well knovn that the brightness 

I extremely high, or to put^ it another way, 

fntly a laiier is able to exceed by many orders 

Lt from a black body operating at the , temperature 

factor mitigating the synchrotron self-absorp- 

" • ' -1 ■ 

fields under consideration the electron orbits; 

' ' ' . ■ - ■ • ' . ' ■ /' 

zed as mentioned earlier, and they do not 

• **.*,■ * . ' 

lie in one of these quantized states. More*- 

' <^ 

h, the magnetic quantum states are far apart^ 
adiation of very short wavelengths could be. 
tern. The longer wavelength radiation, iff 
magnetosphere without absorption, 
cal proposals for a mechanism of piilsar radla- 
ae only very briefly. The subject is both 
te no fully satisfactory theory has emerged, ' 



Those outlined above have varloiii features in f-^iKu^m 'theory , but 
^each has its weak points as well .X: They serve to iliu^trke in so^ie degree 
the wide variety of phehonena that aiust be considered in dbaUng^with the» 
phenomena and p»int out several directions one might expect, a successful 
theoiy eventually to take. L 
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